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COASTAL  MODELING  SYSTEM  (CMS)  USER’S  MANUAL 

CHAPTER  1 

COASTAL  MODELING  SYSTEM  (CMS)  OVERVIEW 

PART  I :  INTRODUCTION 

Backaround 

1.  Numerical  modeling  technology  is  increasingly  being  employed  as  a 
tool  to  study  complex  problems  in  many  engineering  disciplines.  In  the  field 
of  coastal  engineering,  processes  studied  with  numerical  models  include: 
hydrodynamics  and  transport  associated  with  astronomical  tides,  storm  surges, 
and  tsunamis;  refraction,  diffraction,  shoaling  and  breaking  of  waves;  wave- 
induced  currents;  shoreline  change  induced  by  littoral  movement  of  sand; 
erosion  caused  by  short-term  severe  storms;  and  the  fate  and  stability  of 
dredged  material  placed  in  open  water. 

2.  Numerical  models  use  sophisticated  methodology  to  solve  governing 
equations  describing  the  physical  processes  of  interest.  When  a  physical 
process  is  well  described  by  the  specified  equations,  the  numerical  modeling 
technique  is  capable  of  providing  accurate  solutions  to  engineering  problems. 
However,  the  model  developer  must  take  extreme  care  in  implementing  the 
equations  in  a  numerical  scheme,  and  the  model  must  be  thoroughly  tested  for  a 
variety  of  conditions  and  compared  with  prototype  measurements  to  ensure 
correct  model  operation.  Of  equal  importance  is  the  care  the  engineer  must 
take  in  applying  the  model  to  a  field  problem  by  ensuring  that  the  model's 
governing  equations  are  appropriate  for  representing  the  dominant  physical 
processes  at  the  site.  Also,  errors  in  input  specification  could  dramatically 
alter  model  results,  and  without  diligent  inspection  of  model  output,  the 
errors  could  go  undetected. 

3.  The  Coastal  Engineering  Research  Center  (CERC)  has  developed  a 
number  of  numerical  models  that  are  used  for  studying  a  variety  of  coastal 
processes  under  existing  and  proposed  configurations.  Some  of  these  models 
are  large  and  complex  and  may  require  a  substantial  amount  of  input  data 
and/or  computing  effort.  The  models  usually  have  a  battery  of  supporting 
software  for  creating  computational  grids,  running  the  model  on  the  host 
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computer,  and  evaluating  large  amounts  of  output  from  the  simulations. 

Because  of  the  complexity  of  the  model  software  and  the  diversity  of 
modifications  to  input  data  associated  with  individual  applications,  these 
models  have  required  considerable  effort  to  apply. 

4.  Traditionally,  most  of  the  models  were  developed  by  researchers  as 
individual  efforts  to  solve  a  specific  problem.  Consequently,  models  have 
used  different  input/output  structures,  have  been  sparsely  documented,  and 
have  been  frequently  altered  to  accomplish  application- specific  tasks.  In 
addition,  each  model  had  its  own  set  of  supporting  utilities  to  plot  results, 
reformat  data,  etc.  Often  code  alterations  resulted  in  several  versions  of 
the  same  model,  each  with  specific  or  enhanced  capabilities.  This  situation 
had  the  potential  for  introducing  errors,  either  by  engineers  applying  the 
wrong  version  of  the  model  or  by  programmers  inadvertently  altering  the 
operation  of  a  previously  tested  portion  of  the  model. 

5.  Because  of  these  factors,  some  of  the  most  useful  models  at  CERC 
have  also  been  the  most  difficult  to  apply;  and  in  some  cases,  it  was  neces¬ 
sary  for  the  model  developer  to  oversee  the  specific  application.  This 
dependence  on  the  model  developer  meant  that  he  or  she  had  less  time  to  devote 
to  providing  improved  models  to  solve  Corps  problems. 

6.  The  condition  of  the  models  described  above  implied  a  need  for: 

a.  Maintenance  of  a  single  version  of  each  model. 

b.  Standardization  of  the  model  input/output  processes. 

c.  A  central  location  for  using  the  models. 

d.  Versatile  models  capable  of  covering  a  wide  range  of  problems 
at  different  sites  by  using  model  options. 

£.  Comprehensive  model  documentation  in  the  form  of  a  user's 
manual . 

Objective 

7.  The  Coastal  Modeling  System  (CMS)  is  a  software  package  aimed  at 
organizing  CERC's  larger  numerical  models  and  their  supporting  software  into  a 
user-friendly  system  that  is  available  to  all  Corps  elements  having  a  need  to 
apply  the  supported  modeling  technology. 

8.  Several  objectives  are  followed  in  developing  and  expanding  CMS. 
Since  some  of  the  models  share  similar  input  requirements,  output  capability, 
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and  procedural  implementation,  efforts  are  made  to  standardize  these  portions 
of  the  models  as  much  as  possible.  This  standardization  promotes  efficiency 
because  coding  effort  is  reduced,  new  users  learn  the  models  in  the  system 
more  rapidly,  and  chances  for  errors  in  entering  input  or  interpreting  output 
are  reduced  because  of  user  familiarity  with  the  system  structure. 

9.  Most  numerical  models  are  applied  to  specific  areas  by  representing 
the  spatial  bathymetric  and  topographic  features  as  depths  or  elevations  on  a 
matrix  called  a  computational  grid.  In  some  instances,  it  is  possible  to 
share  a  common  numerical  grid  between  models.  This  commonality  allows 
efficient  application  of  several  different  models  to  the  same  site  without 
additional  effort  in  building  new  grids.  Similarly,  processing  of  model 
output  data  can  proceed  in  the  same  manner  for  several  different  models.  To 
the  extent  possible,  models  in  CMS  share  the  same  numerical  grids,  utilities, 
plotting  programs,  post-processing  routines,  and  job  control  files. 

10.  FORTRAN  77  programming  language  is  used  exclusively  in  the  system 
software  to  ensure  portability  of  the  models  and  supporting  programs  to  other 
computer  systems .  Graphics  programs  also  make  use  of  DISSPLA1*1  software . 

When  appropriate,  attempts  are  made  to  produce  a  code  that  is  capable  of  being 
vectorized  for  efficient  and  economical  use  on  vector  array  processing 
computers  such  as  the  CRAY  Y-MP,  where  CMS  resides. 

11.  Models  selected  for  inclusion  in  CMS  must  be  well  advanced  in  their 
development,  and  they  must  have  been  rigorously  tested  over  a  wide  range  of 
conditions .  In  most  cases ,  a  selected  model  has  already  been  applied  to 
numerous  field  problems  and,  thus,  has  reached  some  level  of  maturity  through 
careful  application  to  the  types  of  problems  encountered  by  the  Corps. 
Incorporation  of  these  time -tested  models  into  CMS  involves  modification  to 
the  input/output  structure  and,  possibly,  re-coding  of  sections  to  produce 
more  efficient  use  of  computer  resources .  Once  included  in  CMS ,  the  models 
are  not  expected  to  be  modified  unless  errors  are  found  or  new  features  are 
added.  Therefore,  the  models  in  CMS  can  be  considered  tested,  reliable,  and 
mature . 
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PART  II:  IMPLEMENTATION  OF  CMS 
Organization 

12.  CMS  software  Is  organized  into  three  major  groups:  models, 

supporting  utilities,  and  procedure  files  that  draw  components  together  for 
execution.  Several  software  elements  use  common  algorithms,  and  efforts  are 
made  to  place  these  software  elements  into  shared  libraries .  This  placement 
eliminates  redundant  software  and  reduces  associated  development  costs. 

13.  Another  level  of  organization  concerns  user  interfaces  to  the 
software,  which  are  arranged  according  to  various  activities  encountered 
during  a  modeling  endeavor.  These  major  procedure  files  provide  users  access 
to  the  various  CMS  software  elements  on  the  available  computers . 

14.  Although  the  models  and  supporting  software  are  written  in  standard 
FORTRAN  77,  library  structures  and  procedure  files  are  specific  to  the 
computer  system  hosting  CMS.  In  addition,  some  plotting  utilities  rely  on 
specific  graphics  software  that  resides  on  the  host  computer,  and  certain 
terminal/plotter  configurations  may  be  required  to  produce  plots  of  the  model 
output  at  the  user's  local  site. 

Host  Computer  System 

15.  Most  models  included  in  CMS,  or  targeted  for  addition  to  CMS,  are 
both  memory- intensive  and  computationally  intensive,  requiring  use  of  large 
supercomputers  for  efficient  operation.  It  is  possible  to  run  some  of  the 
models  on  smaller  minicomputers,  but  double  precision  would  probably  be 
necessary  to  avoid  accumulated  round-off  errors,  and  applications  could  take 
from  several  hours  to  several  days  of  CPU  time  on  the  smaller  machines. 

16.  To  satisfy  the  objective  of  Corps -wide  access  to  the  models, 
initial  installation  of  CMS  was  on  a  mainframe  computer  operated  by  Scientific 
Information  Services  (SIS),  formerly  CYBERNET,  which  provided  Corps-wide 
mainframe  and  supercomputer  service.  Now,  with  the  installation  of  a  CRAY  Y- 
MP  at  the  US  Army  Engineer  Waterways  Experiment  Station  (WES),  Corps-wide 
access  is  available,  and  CMS  has  been  transferred  to  that  system  to  take 
advantage  of  reduced  costs  and  the  likelihood  of  a  permanent  home  for  CMS. 
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17.  Corps  personnel  have  access  to  the  CRAY  Y-MP  through  several 
communication  networks.  Presently,  the  supercomputer  can  be  accessed  through 
INTERNET,  MILNET,  ARPANET,  ASNET ,  SURANET,  BITNET,  NSFNET,  through  a  1200, 

2400,  or  9600  baud  moden  or  a  dedicated  line. 

18.  Generally,  CMS  users  are  not  required  to  learn  the  operating  system 
associated  with  the  supercomputer  (UNICOS)  because  aost  of  the  job  control 
commands  normally  required  to  submit  models  and  data  files  to  the  computer  for 
execution  are  accomplished  by  the  CMS  procedure  files.  This  setup  reduces 
learning  time  appreciably  and  minimizes  errors  caused  by  improper  commands. 
However,  users  must  be  able  to  manipulate  files,  create  and  edit  ASCII  files, 
and  download  output  files  to  a  printer  or  plotter.  These  functions  are  easily 
mastered,  and  manuals  are  available  to  all  Corps  users  of  the  CRAY  Y-MP  at  a 
nominal  fee  by  contacting  the  Information  Technology  Laboratory  (ITL)  Research 
Library,  Ms.  Susan  Hicks  (601)  634-2296.  Presently  there  are  26  CRAY  manuals 
covering  such  topics  as  UNICOS  User  Commands,  CF77  Reference,  UNICOS  Support 
Tools,  and  UNICOS  Symbolic  Debugging. 

Model  Support 

19.  Including  a  model  in  CMS  represents  a  technology  transfer  from  CERC 
to  the  field.  CERC  will  maintain  the  CMS  on  the  WES  CRAY  Y-MP  and  will 
provide  support  services  to  Corps  users  of  the  system.  Support  includes 
correcting  recognized  flaws  in  the  codes,  updating  the  models  with  new 
capabilities  and  technology,  improving  the  user  interface  to  the  models, 
improving  graphics  and  visualization  capabilities,  updating  the  User's  Manual 
to  reflect  changes  to  the  models  and/or  CMS,  conducting  periodic  workshops  for 
Corps  personnel,  and  providing  telephone  support  services. 

20.  Additionally,  CERC  staff  can  assist  Corps  personnel  in  applications 
of  the  CMS  via  "one-stop  services"  or  by  direct  participation  in  site-specific 
studies.  One-stop  service  is  intended  to  address  questions  or  problems  that 
arise  during  field  application  of  a  model  in  CMS.  Usually,  these  questions 
can  be  satisfactorily  resolved  in  a  short  time  over  the  telephone.  More 
involved  questions  requiring  a  substantial  effort  by  CERC  engineers  or 
scientists  may  require  reimbursement.  Experience  at  CERC  indicates  that  field 
application  of  these  models  usually  requires  a  significant  initial  consulting 
effort  by  CERC  engineers  until  experience  has  been  gained  by  the  field  user. 
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Coastal  Modeling  System  (CMS)  User's  Manual 


21.  The  Coastal  Modeling  System  (CMS)  User's  Manual  is  intended  to  be 
an  evolving  document,  and  it  is  structured  in  a  modular  fashion,  much  like  the 
modeling  system  itself.  Individual  numerical  models  and  major  supporting 
utility  software  are  documented  in  separate  chapters.  Attempts  are  made  to 
structure  all  chapters  in  a  similar  format  to  facilitate  learning  the  system 
models . 

22.  The  unbound  format  of  the  user's  manual  allows  efficient  and  cost- 
effective  updating  of  the  manual  as  models  are  added  to  CMS,  and  it  allows 
users  to  remove  chapters  for  convenient  reference  during  model  applications. 
The  documentation  for  each  new  model  will  be  an  added  chapter  to  the  CMS 
User's  Manual.  Updates  will  also  include  additions  of  (or  alterations  to) 
utilities  and  procedures. 

23.  Initial  distribution  of  the  CMS  User's  Manual  will  be  to  all  Corps 
Divisions  and  Districts  with  coastal  interests.  A  register  of  all  manuals 
distributed  within  the  Corps  of  Engineers  will  be  maintained  by  CERC,  and 
updates  will  be  provided  for  all  Corps-registered  copies  of  the  manual. 

24.  Training  on  the  usage  of  CMS  and  on  application  of  specific  models 
within  CMS  is  accomplished  during  periodic  workshops.  Workshop  participants, 
especially  new  users,  will  be  introduced  to  CMS.  Each  workshop  will  demon¬ 
strate  sign- on  procedures,  building  input  data  files,  file  transfer  methods, 
accessing  CMS,  running  workshop-specific  models,  and  post-processing  model 
results.  Technical  presentations  of  workshop  specific  models  will  also  be 
given.  More  intensive  training  can  be  provided  at  CERC  as  part  of  joint  field 
applications  between  CERC  and  personnel  from  a  field  activity. 

Point  of  Contact 

25.  Each  model  residing  in  CMS  has  a  CERC  point  of  contact  (POC).  Most 
often  that  person  is  the  model  developer  or  someone  with  extensive  experience 
in  applying  the  model.  Table  1-1  provides  POC's  for  models  presently  included 
in  CMS.  This  table  will  be  updated  periodically  to  assure  that  it  continues 
to  be  a  useful  reference  for  CMS  users.  The  modules  listed  in  column  1  of 
Table  1-1  are  briefly  described  in  Part  III  of  this  chapter  and  are  extensive¬ 
ly  documented  in  later  chapters. 
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Table  1-1 


CERC  Points  of  Contact 


Sublect 

Point  of  Contact 

Office  Symbol 

Phone  Number 

CMS  general 

Inquiries 

Ms.  Mary  A.  Cialone 

CEWES-CR-P 

601-634-2139 

Using  CMS 

Ms .  Lucy  W .  Chou 

CEWES-CR-P 

601-634-2843 

SPH 

Mr .  Dave  J .  Mark 

CEWES-CR-0 

601-634-2094 

WIFM 

Ms.  Mary  A.  Cialone 

CEWES-CR-P 

601-634-2139 

RCPWAVE 

Mr.  Steven  M.  Bratos 

CEWES-CR-0 

601-634-4230 

CLHYD 

Ms.  Mary  A.  Cialone 

CEWES-CR-P 

601-634-2139 

SHALWV 

Dr.  Robert  E.  Jensen 

CEWES-CR 

601-634-2101 

STWAVE 

Mr.  Steven  M.  Bratos 

CEWES-CR-0 

601-634-4230 

HARBD 

Dr.  Edward  F.  Thompson 

CEWES-CR-0 

601-634-2027 

GENESIS 

Mr.  Mark  B.  Gravens 

CEWES-CR-P 

601-634-3809 

SBEACH 

Mr.  Randall  A.  Wise 

CEWES-CR-P 

601-634-3085 

PBLWIND 

Dr.  Edward  F.  Thompson 

CEWES-CR-0 

601-634-2027 

HURWIN 

Dr.  Edward  F.  Thompson 

CEWES-CR-0 

601-634-2027 

WICM 

Ms.  Jane  M.  Smith 

CEWES-CR-P 

601-634-2079 

CMSGRID 

Ms.  Mary  A.  Cialone 

CEWES-CR-P 

601-634-2139 

CMSUTIL 

Ms.  Lucy  W.  Chou 

CEWES-CR-P 

601-634-2843 

CMSPOST 

Ms.  Lucy  W.  Chou 

CEWES-CR-P 

601-634-2843 

CMS  SAMP 

Ms.  Lucy  W.  Chou 

CEWES-CR-P 

601-634-2843 

PART  III:  PRESENT  CHS  COMPONENTS 


26.  The  following  briefly  describes  the  numerical  models,  major  utility 
software,  and  major  procedures  currently  in  the  Coastal  Modeling  System. 

Cbapttr  .2.; _ Using  the  Coastal  Modeling  System 

27.  This  chapter  provides  information  on  execution  of  the  CMS  on  the 
WES  CRAY  Y-MP.  The  new  user  should  refer  to  this  chapter  to  learn  to  compile 
or  run  a  model  for  a  specific  application.  Once  the  user  becomes  familiar 
with  the  system,  this  chapter  can  be  used  as  a  quick  reference. 

Chanter  3:  Standard  Project  Hurricane  (SPH) 

28.  The  numerical  model  SPH  is  a  parametric  model  for  representing  wind 
and  atmospheric  pressure  fields  generated  by  hurricanes.  It  is  based  on  the 
Standard  Project  Hurricane  criteria  developed  by  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA  1979),  and  the  model's  primary  output  is 
hurricane -generated  wind  fields  that  can  be  used  in  storm  surge  modeling.  It 
can  be  run  separately,  or  it  can  be  invoked  from  within  the  WES  Implicit 
Flooding  Model  (WIFM) . 

Chapter  4; _ WES  Implicit  Flooding  Model  (WIFM) 

29.  The  numerical  model  WIFM  solves  the  vertically  integrated  Navier- 
Stokes  equations  in  stretched  Cartesian  coordinates.  The  model  simulates 
shallow-water ,  long-wave  hydrodynamics  such  as  tidal  circulation,  storm 
surges,  and  tsunami  propagation.  WIFM  contains  many  useful  features  for 
studying  these  phenomena,  such  as  moving  boundaries  to  simulate  flooding/ 
drying  of  low-lying  areas  and  subgrid  flow  boundaries  to  simulate  small 
barrier  islands,  jetties,  dunes,  or  other  structural  features.  The  model  may 
be  driven  at  the  outer  boundary  by  tide  elevations,  flow  velocities,  specifi¬ 
cation  of  uniform  flux,  or  inverted  barometer  effects.  WIFM  also  accepts  wind 
fields  for  including  the  effects  of  wind  stress  during  hurricanes  or  other 
strong  storm  systems. 
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Chapter  5: _ Regional  Coastal  Processes  Wave  Propagation  Model  (RCPWAVE) 

30.  The  numerical  model  RCPWAVE  is  a  short-wave  model  used  to  predict 
linear,  plane  wave  propagation  over  an  open  coast  region  of  arbitrary  bathyme 
try.  RCPWAVE  uses  linear  wave  theory  because  it  has  been  shown  to  yield 
fairly  accurate  first-order  solutions  to  wave  propagation  problems  at  a 
relatively  low  cost.  Refractive  and  bottom- induced  diffractive  effects  are 
included  in  the  model;  however,  the  model  cannot  treat  diffraction  caused  by 
surface-piercing  structures.  This  model  does  not  include  nonlinear  wave 
effects  or  a  spectral  representation  of  irregular  waves . 


Chapter  6:  Curvilinear  Lonp-Wave  Hydrodynamic  Model  (CLHYD) 


31.  The  numerical  model  CLHYD  simulates  shallow-water ,  long-wave 
hydrodynamics  such  as  tidal  circulation  and  storm  surge  propagation.  CLHYD 
can  simulate  flow  fields  induced  by  wind  fields,  river  inflows/outflows,  and 
tidal  forcing.  CLHYD  is  similar  to  WIFM,  with  the  added  feature  of  operating 
on  a  boundary- fitted  (curvilinear)  grid  system.  However,  CLHYD  cannot 
simulate  flooding/drying  of  low-lying  areas  as  WIFM  can.  This  feature  will  be 
incorporated  in  a  later  release  of  CLHYD. 


Chapter  7:  Spectral  Wave  Modeling  Module:  Model  SHALWV 


32.  The  numerical  model  SHALWV  is  a  time -dependent  spectral  wind-wave 
model  for  computing  a  time -history  of  wind- generated  waves.  The  model  solves 
the  inhomogeneous  energy  balance  equation  using  finite-difference  methods.  It 
simulates  the  growth,  decay,  and  transformation  of  a  wave  field  over  a  spatial 
area  (i.e.,  an  ocean  basin,  bay,  or  lake)  for  a  given  time  period.  SHALWV  can 
simulate  the  wave  climate  for  a  specific  storm  or  idealized  events,  such  as  a 
standard  project  hurricane. 

Chanter  8:  Spectral  Wave  Modeling  Module:  Model  STWAVE 


33.  The  numerical  model  STWAVE  is  a  near-coast,  time- independent 
spectral  wave  energy  propagation  model.  The  model  solves  the  spectral  energy 
balance  equation  (including  refraction,  shoaling,  and  wave  breaking)  using 
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finite -difference  methods.  This  steady-state  model  simulates  wave  propagation 
over  a  spatial  area  assuming  wave  conditions  vary  sufficiently  slowly.  The 
variation  of  waves  at  a  given  point  may  be  neglected  relative  to  the  time 
required  for  waves  to  pass  across  the  computational  grid  if  the  model  is 
limited  to  near-coast  applications  in  which  waves  move  quickly  across  the  grid 
(within  30  minutes). 

Chanter  9:  Harbor  Wave  Oscillation  Model  (HARBD 

34.  The  numerical  model  HARBD  is  a  finite  element  model  for  predicting 
wave  oscillations  in  harbors.  HARBD  is  a  steady-state,  linear  monochromatic 
wave  model  that  assumes  bathymetric  gradients  are  small  and  neglects  wave -wave 
interaction,  wave-current  interaction,  wave  breaking,  and  wave  transmission 
and  overtopping  of  structures .  The  model  has  been  used  in  the  design  and 
modification  of  numerous  harbors,  the  study  of  dredging  effects  on  wave 
propagation,  and  in  the  design  and  planning  of  wave  protection  structures  for 
existing  harbors. 

Chapter  IQ; _ Generalized  Model  for  Simulating  Shoreline  Change  (GENESIS) 

35.  The  numerical  model  GENESIS  is  a  one -dimensional  ("one-line") 
model  for  simulating  long-term  shoreline  change.  It  can  predict  shoreline 
change  and  longshore  transport  rates  under  a  wide  range  of  beach,  coastal 
structure,  wave,  initial  and  boundary  conditions,  which  may  vary  in  space  and 
time.  Input  data  include  the  initial  shoreline  position,  measured  shoreline 
position  for  calibration  purposes,  structure  positions,  depths  along  the 
nearshore  reference  line  (where  input  from  a  wave  model  are  saved) ,  and  the 
wave  height,  period,  and  direction  for  every  time-step.  Model  output  includes 
the  shoreline  position  and  longshore  transport  rates  at  user-specified  time- 
steps  . 


Chapter  11:  Storm- Induced  Beach  Change  Model  (SBEACH) 

36.  The  numerical  model  SBEACH  is  an  empirically-based  two-dimensional 
model  for  predicting  short-term,  storm- induced  beach  erosion  and  post-storm 
recovery.  Short-term  erosion  is  treated  as  a  process  dominated  by  cross-shore 
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sand  transport  processes.  Input  requirements  include  a  time- series  of  wave 
height  and  period,  a  time-series  of  water  level,  median  beach  grain  size,  and 
initial  profile  shape.  Optionally,  model  input  may  include  the  presence  of  a 
seawall  or  revetment,  and  SBEACH  simulates  the  profile  response  to  such  a 
structure.  Model  output  includes  the  beach  profile  at  all  user-specified 
time-steps,  as  well  as  the  cross-shore  distribution  of  various  process  parame¬ 
ters  (maximum  wave  height  and  water  level,  and  wave  height  and  water  level  at 
user-specified  intervals). 

Chapter  12:  Extratropical  Storm  Planetary  Boundary  Laver  Wind  Model  (PBLWIND) 

37.  The  numerical  model  PBLWIND  is  a  generalized  numerical  model  used 
to  predict  winds  near  the  water  surface  based  on  atmospheric  pressure  gradi¬ 
ents  and  temperature  differences  between  the  air  and  water.  PBLWIND  is  a 
steady  state  model  which  includes  computational  modules  for  geostrophic, 
gradient,  and  surface  boundary  layer  winds.  The  model  should  not  be  used  for 
modeling  of  hurricanes  and  tropical  storms  because  their  compact  size,  intense 
pressure  gradients,  and  rapid  changes  in  time  require  special  treatment. 

Output  from  PBLWIND  consists  of  wind  speed  and  direction  at  a  desired  eleva¬ 
tion  above  the  water  surface  which  can  be  used  as  input  to  a  hydrodynamic 
model  (SHALWV,  WIFM,  CLHYD) . 

Chapter  13:  Tropical  Storm  Planetary  Boundary  Laver  Wind  Model  (HURWIN) 

38.  The  numerical  model  HURWIN  is  a  two-dimensional,  time -dependent 
model  for  predicting  surface  stress  and  wind  speed  and  direction  in  the 
planetary  boundary  layer  of  a  tropical  cyclone.  Wind  information  is  calculat¬ 
ed  from  meteorological  storm  parameters  available  for  historical  hurricanes 
and  provided  at  a  user-specified  elevation.  The  model  is  based  on  the 
momentum  equations  that  are  vertically-averaged  through  the  depth  of  the 
planetary  boundary  layer.  Options  are  also  provided  for  estimating  surface 
wind  over  terrain  of  specified  roughness  including  lakes,  marshes,  plains, 
woods,  and  cities.  It  has  been  used  extensively  in  the  Wave  Information 
Studies  (WIS)  to  hindcast  historical  hurricanes  along  U.S.  coasts. 
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Chapter  14:  Wave -Induced  Current  Model  (WICM) 

39.  The  numerical  model  UICM  is  a  two-dimensional,  depth- averaged 
model  for  predicting  wave- induced  currents  and  water  surface  setup.  UICM  can 
simulate  flow  fields  induced  by  waves,  winds,  river  inf lows /outf lows ,  and 
tidal  forcing.  Similar  to  CLHYD,  this  finite  difference  model  is  developed  in 
boundary- fitted  (curvilinear)  coordinates  and  cannot  simulate  flooding  and 
drying  of  low- lying  areas . 

Appendix  A:  CMSGRID 

40.  Grid  development  is  a  major  part  of  successfully  applying  a 
numerical  model  to  a  specific  site.  Module  CMSGRID  contains  software  used  in 
the  generation  of  stretched  coordinate,  rectilinear  computational  grids  for 
several  models  in  CMS.  The  software  employs  sophisticated  techniques  that 
allow  concentration  of  grid  cells  in  regions  of  interest,  or  where  geographic 
features  are  irregular,  and  wider  spacing  of  grid  cells  in  areas  where 
conditions  are  not  expected  to  change  rapidly.  The  ability  to  generate 
variably  spaced  grids  provides  economy  in  computational  time  and  costs. 
Generated  grids  can  be  plotted  to  scale  on  Mylar  for  overlaying  bathymetric 
charts  to  obtain  depths  and  elevations. 

Appendix  B:  CMSUTIL 

41.  Module  CMSUTIL  contains  useful  programs  that  supplement  numerical 
models  in  CMS.  Presently,  there  are  two  programs  in  this  module:  a  program 
to  determine  tidal  constituents  from  a  time -history  of  tidal  elevations,  and  a 
program  to  generate  a  time  series  of  water  elevations  from  tidal  constituent 
input . 


Appendix  C:  CMSPOST 

42.  Most  numerical  models  generate  large  output  files  containing 
results  saved  at  user-specified  grid  cells  and  time-steps  during  the  simula¬ 
tion.  Module  CMSPOST  contains  post-processing  plotting  packages  that  plot  the 
model  output  for  comparison  and  analysis  purposes.  Four  basic  types  of 
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plotting  are  available:  (a)  time-histories  of  field  arrays  such  as  water 
surface  elevation  or  velocity  at  selected  grid  points,  (b)  "snapshots"  of  the 
entire  flow  field  at  a  given  instant  in  time,  (c)  wave  ray  plots,  and 
(d)  profile  plots  that  show  the  spatial  variation  of  a  model  variable  at  an 
instant  in  time. 


43.  Module  CMSSAMP  is  used  to  access  sample  input  files  for  each  of  the 
models  in  the  CMS.  This  provides  the  user  with  the  required  format  of  input 
data,  and  can  also  be  used  as  a  template  for  generating  the  user's  own  input 
file(s).  In  addition,  sample  files  can  be  used  to  run  a  particular  model  in 
the  CMS  to  help  the  user  gain  familiarity  with  the  model.  Appendix  D  contains 
a  table  of  sample  input  file  names  and  the  procedure  for  accessing  these  files 
from  within  the  CMS. 
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HURW1N:  TROPICAL  STORM  PLANETARY  BOUNDARY  LAYER  WIND  MODEL 

PART  I :  INTRODUCTION 

1.  The  HURWIN  model  in  the  Spectral  Wave  Modeling  Module  of  the  Coastal 
Modeling  System  (CMS)  computes  surface  stress  and  wind  speed  and  direction  In 
the  planetary  boundary  layer  of  a  tropical  cyclone.  Wind  information  is  cal¬ 
culated  from  meteorological  storm  parameters  available  for  historical  hurri¬ 
canes  and  provided  at  a  user-specified  elevation.  The  full  time  history  of 
the  surface  wind  field  is  described  by  linear  interpolation  in  time  using  a 
series  of  characteristic  wind  field  "snapshots."  Snapshots  are  calculated  at 
discrete  times  in  a  storm's  history,  assuming  the  storm  can  be  represented  as 
a  series  of  steady-state  configurations. 

2.  The  model,  developed  by  Cardone,  Greenwood,  and  Greenwood  (1992)  and 
upgraded  by  Cardone  et  al.  (1994),  is  based  on  a  numerical  primitive  equation 
model  of  the  planetary  boundary  layer  in  a  moving  tropical  cyclone.  A  surface 
drag  formulation,  based  on  Arya's  (1977)  similarity  model,  is  coupled  with  a 
roughness  parameter  specification  for  a  water  surface  consistent  with 
Cardone' s  (1969)  law.  Options  are  also  provided  for  estimating  surface  wind 
over  terrain  of  specified  roughness  including  lakes,  marshes,  plains,  woods, 
and  cities. 

3.  The  CMS  includes  another  hurricane  wind  model,  the  Standard  Project 
Hurricane  (SPH,  Chapter  3).  Although  the  SPH  model  and  the  HURWIN  model 
provide  similar  information,  the  SPH  model  has  traditionally  been  used  in  the 
Corps  of  Engineers  (CE)  as  input  to  storm  surge  models  (e.g.  CLHYD,  Chapter  6) 
and  the  HURWIN  model  with  wave  models  (SHALWV,  Chapter  7).  The  HURWIN  model 
is  more  flexible  and  more  nearly  represents  the  physical  processes  in  a 
hurricane.  It  has  been  used  extensively  in  the  Wave  Information  Studies  (WIS) 
to  hindcast  historical  hurricanes  along  U.S.  coasts  (Abel  et  al.  1989).  The 
CMS  will  allow  the  two  models  to  be  more  easily  evaluated  and  compared  in 
various  applications.  A  broad  overview  of  wind  processes  and  CE  modeling 
applications  is  given  by  Thompson  and  Leenknecht  (1994). 

4.  Theoretical  background  of  HURWIN  is  described  in  the  following 
section.  Operation  of  the  HURWIN  portion  of  the  CMS  Spectral  Wave  Modeling 
module  is  described  in  the  remaining  sections.  The  structure  of  the  module  is 


described  in  Part  III.  Organizational  procedures  for  a  HURWIN  s inula t ion  are 
described  in  Part  IV.  Part  V  covers  the  required  data  for  a  successful  HURWIN 
simulation.  An  example  problem,  simulation  of  winds  in  Typhoon  Russ,  is 
presented  in  Part  VI  to  illustrate  application  of  the  model. 


PART  II:  THEORETICAL  FOUNDATION 


Overview 


5.  The  basic  model  is  based  on  a  vortex  model  developed  by  Chow  (1971) 
and  modified  by  Cardone,  Greenwood,  and  Greenwood  (1992).  The  basic  model  was 
further  modified  by  Cardone  et  al.  (1994)  to  include  options  for  more  detailed 
spatial  resolution  and  more  general  specification  of  radial  pressure 
variation.  Chow's  model  concerns  the  planetary  boundary  layer  (PBL)  only  and 
solves  for  the  wind  field  and  horizontal  convergence  in  the  PBL  of  a  moving 
tropical  cyclone  from  the  equations  of  motion.  The  pressure  field  in  the 
boundary  layer  is  prescribed  and  fixed,  so  that  there  are  no  atmospheric 
gravity  waves  excited  in  the  numerical  solution.  This  facilitates  the  use  of 
a  nested  grid  system,  which  allows  grid  spacings  as  small  as  2  km  near  the 
hurricane  inner  region  without  sacrifice  of  overall  computational  efficiency. 

6.  The  model  is  based  upon  the  equation  of  horizontal  motion,  vertical¬ 
ly  averaged  through  the  depth  of  the  PBL,  written  in  coordinates  fixed  to  the 
earth  as: 
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p  -  pressure 

JCH  -  horizontal  eddy  viscosity  coefficient 


CD  -  drag  coefficient 
h  -  depth  of  the  planetary  boundary  layer 

It  is  assumed  that  the  vertical  advection  of  momentum  is  small  compared  to  the 
horizontal  advection  and  can  be  neglected  and  that  the  shearing  stress 
vanishes  at  the  top  of  the  FBL. 

7.  Pressure  is  prescribed  as  the  sum  of  pe  and  p 


P  m  Pc  +  P 


(13-2) 


where 


pc  -  pressure  field  representing  the  tropical  cyclone;  not 

necessarily  axisymmetric;  assumed  to  translate  with  the  storm  at 
a  specified  velocity 

p  —  large-scale  pressure  field,  which  may  be  specified  by  the 
corresponding  constant  geos trophic  flow  as 

fj?x  •  -  j  V  p  (13-3) 

With  this  pressure  specification,  Equation  13-1  may  be  rewritten  as  follows: 


^  *  fit  x  (?-  Vg)  .  -  i  vPc  ♦  v.(k„v?)  -  (13'4) 

A  moving  Cartesian  coordinate  system  (x,y)  is  now  defined  such  that  its  origin 
always  coincides  with  the  moving  low  center  of  pc  .  In  terms  of  the  moving 
system,  Equation  13-4  is  transformed  into 
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horizontal  wind  velocity  relative  to  the  low  center;  «  P  -  Pt 
effective  geostrophic  flow  relative  to  the  low  center;  »  P  t  -  Pt 
velocity  of  the  moving  reference  system  relative  to  the  fixed  earth 
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time  change  local  to  the  moving  coordinate  system; 


8.  Equation  13-5  can  be  expanded  from  vector  form  into  equations 
involving  the  scalar  x-  and  y-components  of  velocity.  After  some  rearranging 
of  terms,  the  equations  to  be  solved  are 
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fv  -  P  +  H-  F. 


(13-6) 
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9.  The  general  formulation  is  completed  with  the  specification  of  the 
fora  of  CD,  Kh  and  the  boundary  condition  at  the  outermost  boundary  of  the 
grid.  Following  Saagorinsky  (1963)  the  horizontal  eddy  viscosity  coefficient 
is  given  by 


(13-8) 


where 

a  -  mesh  size 

k  -  non- dimensional  constant  (*  -  0.4  is  assumed) 

The  term  in  brackets  in  Equation  13-8  may  be  viewed  as  a  total  deformation  due 
to  tension  and  shearing  strains.  Drag  coefficient  is  computed  based  on  a 
revised  PBL  parameterization  discussed  later  in  this  chapter.  The  boundary 
condition  at  the  outermost  boundary  of  the  grid  is  obtained  by  neglecting 
acceleration  and  horizontal  diffusion  of  momentum,  implying  a  balance  between 
Coriolis,  pressure  gradient,  and  surface  friction  forces.  With  this 
assumption,  Equation  13-5  can  be  simplified  along  the  outer  boundary  to 

fit  x  (V-?g)  *  -  lvPc-  !g\V+Vcl(  P  +Vc)  (13-9) 

10.  Several  options  are  available  for  specifying  the  pressure  field. 

The  simplest  and  most  used  is  an  axisymmetric  pressure  field  defined  by  the 
well-known  exponential  pressure  law.  This  expression,  which  is  also  used  in 
the  SPH  (Chapter  3) ,  is  given  by 

pc(r )  *  p,  +  a pe~*S'  (13-10) 

where 

p0  -  storm  central  pressure 
ap  -  storm  pressure  anomaly  defined  as 
R f  -  scale  radius 
r  -  radial  distance  from  the  eye 

Ap  -  p„  -  p0  (13-11) 

where  p.  -  axisymmetric  ambient  (far- field)  pressure 
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11.  A  more  general  option  is  available  in  which  the  parameters  ap  and 
R p  in  Equation  13-10  can  be  specified  by  stora  quadrant.  The  final  pressure 
field  under  this  option  is  devised  by  smoothing  between  the  pressure  variation 
specified  for  each  quadrant. 

12.  Some  storms  are  not  well  represented  by  Equation  13-10,  even  when 
the  parameters  are  varied  by  quadrant.  Options  for  more  general 
specifications  of  radially  symmetric  pressure  variation  were  introduced  by 
Cardone  et  al.  (1994).  A  more  general  form  of  Equation  13-10,  which  has 
proved  successful  in  describing  some  storms  (Holland  1980) ,  is  given  by 

pc(r)  (13-12) 

where  B  is  a  constant  in  the  general  range  0.5 -2. 5. 

13.  Field  evidence  has  accumulated  in  recent  years  that  the  radial 
pressure  profile  in  the  inner  core  of  some  storms  is  so  irregular  that  even 
Equation  13-12  provides  a  poor  approximation.  Willoughby  (1990)  and  Black  and 
Willoughby  (1992)  describe  the  tendency  for  concentric  rings  in  the  radial 
wind  distribution  to  be  a  fairly  typical  characteristic  of  intense  tropical 
cyclones.  The  rings  appear  to  be  related  intimately  to  storm  intensity 
evolution.  A  capability  for  modeling  radial  wind  distributions  with  broad  or 
double  maxima  can  be  added  by  extending  the  generalization  of  Equation  13-12 
to  a  double  exponential  form  (Cardone  et  al.  1994) 

pc(r)  =  p0  +  ^CdPi  e~lgpi/l) "  (13-13) 

where 

dp,  -  pressure  anomaly  for  the  i ' th  component 

R,,  -  scale  radius  for  the  i '  th  component 

B,  -  Holland's  B  coefficient  for  the  i'th  component 
Since  the  values  of  dp,  are  constrained  by  the  relationship 


(13-14) 


only  the  values  of  Ap  and  dp,  need  be  specified  to  determine  the  value  of  dp}. 

14.  The  boundary  layer  depth  h  was  investigated  by  Cardone, 
Greenwood,  and  Greenwood  (1992)  for  two  historical  hurricanes  (Camille  and 
Carla) .  They  concluded  that  an  appropriate  boundary  layer  depth  is 
approximately  500  m. 
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15.  In  sunaary,  the  main  Independent  variables  are: 
p,  -  central  pressure 

ap  -  storm  pressure  anomaly 
Rj.R,,  -  scale  radius  values 
B,B,  -  Holland's  (1980)  8  coefficients 
dp,  -  pressure  anomaly  for  inner  component 
u„,  vc  -  components  of  storm  forward  speed 
h  —  boundary  layer  depth 
u(l  vt  -  components  of  geostrophic  steering  flow 
Either  one  or  four  values  (by  quadrant)  may  be  provided  for  each  of  the  first 
three  variables  if  the  simple  pressure  profile  is  used  (Equation  13-10). 

16.  The  computational  grid  is  a  system  of  rectangular  nested  grids 
within  each  of  which  the  mesh  is  constant.  Figure  13-1  shows  the  inner  three 
nests  in  one  quadrant  of  the  grid  system.  Mesh  size  doubles  with  each 
successive  nest,  and  each  mesh  extends  10  mesh  spacings  from  the  center.  The 
number  of  nests  is  at  most  seven  and  generally  at  least  five.  If  five  nests 
are  used  and  the  mesh  size  of  the  innermost  nest  is  5  km,  the  second  through 
fifth  mesh  sizes  are  10,  20,  40,  and  80  km,  respectively,  and  the  entire  grid 
covers  a  square  1,600  km  on  a  side.  If  seven  nests  are  used  and  the 
resolution  of  the  innermost  nest  is  refined  to  2  km,  the  other  mesh  sizes 
(nests  two  through  seven)  are  4,  8,  16,  32,  64,  and  128  km,  and  the  grid 
covers  a  2,560-km  square. 


EXAMPLE  OF  NESTED  GflD  SYSTEM  (1ST  QUADRANT) 
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17.  Details  of  the  finite  difference  formulation  and  of  the  compu¬ 
tational  scheme  are  given  by  Chow  (1971)  and  are  not  repeated  in  detail  here. 
Basically,  a  combination  of  diagonal  and  ordinary  upstream  differencing  is 
used  for  spatial  derivatives  in  order  to  reduce  computational  errors  in  the 
calculation  of  the  advection  terms  and  at  the  interior  mesh  boundaries.  The 
computation  starts  with  an  initial  guess  field  consisting  of  the  gradient  wind 
components  (computed  from  pc) .  At  each  grid  point,  Equations  13-6  and  13-7 
are  integrated  forward  in  time  until  the  acceleration  (3?/3t)c  is  tolerably 
small.  The  computational  time -step,  in  seconds,  is  taken  as  ten  times  the 
finest  computational  mesh  size,  in  kilometers.  For  example,  if  the  innermost 
grid  mesh  size  is  5  km,  the  time-step  is  50  sec.  Chow  found  that  800  itera¬ 
tions  (equivalent  to  13  hr  20  min  if  the  time- step  is  60  sec)  were  sufficient 
to  achieve  a  steady- state  solution.  Cardone  et  al.  (1994)  recommended  3,200 
iterations  when  using  the  generalized  pressure  specification. 

18.  The  drag  law  used  by  Chow  (1971)  was  revised  by  Cardone,  Greenwood, 
and  Greenwood  (1992)  to  provide  an  accurate  and  consistent  PBL  wind  and  stress 
representation  in  the  vortex  model. 

19.  The  framework  for  parameterization  of  the  fluxes  of  momentum,  heat, 
and  moisture  in  the  PBL  was  developed  by  Blackadar  and  Tennekes  (1968)  and 
Zilitinkevich  (1969).  The  parametric  relations  result  from  the  matching  of 
mean  profiles  of  wind,  temperature,  and  moisture  predicted  by  surface  and 
outer- layer  similarity  theories  for  a  PBL  in  which  the  flow  is  assumed  to  be 
horizontally  homogeneous  and  quasi -stationary. 

20.  A  particularly  convenient  form  of  the  parameterization,  first  pro¬ 
posed  by  Deardorff  (1972),  expresses  the  PBL  fluxes  in  terms  of  layer -averaged 
mean  PBL  properties.  Deardorff' s  parameterization,  combined  with  Cardone 's 
roughness  parameterization,  has  been  found  to  provide  reasonable  results  for 
hurricanes.  The  parameterization  adapted  in  this  model  is  taken  from  Arya's 
(1977)  update  of  Deardorff' s  scheme. 

21.  The  general  form  of  the  parametric  relations  may  be  written: 


feu 

u7 


-(In  z.  +  A„) 


(13-15) 
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Ifl 


(13-16) 


kv 

u. 


~K 


A*'.'  M  -  -(Ini.  *Cm)  (13-17) 

hh L~-qi  -  -  (In  i.  +  Dm)  (13-18) 


where 


k 


u, 


v 


u. 

A 

Z. 

A..  Bm,  Cm,  Dm 

f/l  fl 

0. 


q. 

0. 


q 

q- 


-  von  Karmen's  constant 

-  vertically  integrated  (as  in  Equations  13-6  and  13-7) 
Horizontal  wind  components  (in  the  direction  of  the 
surface  shear  and  perpendicular  to  it) 

-  friction  velocity 

-  roughness  parameter  normalized  by  the  PBL  height  (zyh) 

-  universal  functions  of  dimensionless  similarity 
parameters 

-  sign  of  f  (+  or  -) 

-  mean  layer  virtual  potential  temperature  at 
desired  elevation 

-  potential  temperature  and  specific  humidity  at  z„ 

-  potential  temperature  scale  expressed  in  terms 
of  the  heat  flux 

-  specific  humidity  at  desired  elevation 

-  a  specific  humidity  scale  involving  the  moisture  flux 


22.  The  Monin-Obukov  length  L  may  be  expressed  in  terms  of  0,  and 
u.  since 


L  -  ~U:VC'  -  (13-19) 

kgH  g9. 

where  H  is  the  heat  flux. 

23.  Two  competing  theories  exist  for  the  form  of  the  universal 
functions.  In  one  theory,  known  as  the  Rossby-number  similarity  theory,  the 
boundary  layer  height  is  uniquely  determined  by  u./f  and  L.  For  the  near 
neutral  hurricane  PBL,  that  theory  predicts  the  PBL  to  increase  as  the  ratio 
u ,/f  increases  toward  the  center  of  storms,  in  apparent  contradiction  to 
observation. 
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24.  In  the  generalized  theory,  the  depth  of  the  PBL  h  is  specified  as 
an  independent  variable.  Arya  (1977)  presents  expressions  for  the  similarity 
functions  in  terms  of  this  generalized  theory  as  follows: 


ln(-  *15 

^  g  fh  e0 Jk/L 

u. 

1„(-  |)  .  3.7 

( unstable ) 

(13-20) 

-0.96*  +  2.5 

0.80*  +  1.1 

x  £  +2 

-2.00  *  +4.7 

(stable) 

(13-21) 

For  near-neutral  conditions,  -2  <  h/L  <2,  A_  ,  Bm  ,  and  Cm  are  assumed  to 
be  given  by  linear  interpolation  between  the  above  computed  values  at  h/L  - 
±2. 

25.  In  terms  of  the  similarity  relations  (Equations  13-15  through 
13-18),  the  drag  coefficient  with  respect  to  the  integrated  PBL  wind  is- 


_ ic* _ 

(In  z,  +  Ajf  +  Bl 


(13-22) 


while  the  angle  ft,  between  the  surface  wind  and  the  integrated  PBL  wind,  is 


P  *  tan  '(v/u) 


(13-23) 


26.  To  incorporate  the  similarity  theory  in  the  hurricane  model,  two 
cases  were  considered: 

a.  Land.  In  a  PBL  over  land,  the  following  parameters  are  pre¬ 
scribed:  f,  z„,  h,  0,  -  0O.  The  parametric  relations  then 
define  the  following  functions 


CD  •  F,(  \V\) 


f  ,z0,h,9,-90 


(13-24) 


fi 
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f  ,z„h,9-9t 


(13-25) 
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b.  Water.  Over  water,  the  roughness  parameter  is  not  known  but 
can  be  expressed  in  terms  of  u.  through  a  Charnock-type  rela¬ 
tion 

z.a^l  (13-26) 

*  g 

or  the  form  proposed  by  Cardone  (1969) .  The  parametric 
relations  can  then  be  solved  for 

C.-FJ(|*|)  (13-27) 

0-F.(l?i)  (13-28) 

i. 

27.  The  Charnock  constant  a  and  the  constant  k  are  assignable  in 
the  model.  However,  recommended  values,  based  on  the  work  of  Garratt  (1977) 
and  Cardone,  Greenwood,  and  Greenwood  (1992),  are 

a  -  0.035 

k  -  0.35 

28.  Some  assumptions  are  introduced  to  simplify  solution  of  the 
boundary  layer.  It  is  assumed  that 

*£  •  1  (13-29) 

u. 

so  that  the  similarity  functions  in  Equation  13-20  depend  only  on  the  ratio 
h/L.  Further,  it  is  assumed  that  over  water  the  air- sea  temperature 
difference  and  boundary  layer  height  can  be  considered  to  be  invariant  over 
the  domain  of  the  storm.  This  appears  to  be  a  reasonable  approximation. 

Except  for  hurricanes  crossing  major  ocean- current  boundaries,  the  assumption 
of  u.'  ~zontal  homogeneity  of  9 ,  -  60  also  seems  reasonable,  especially  for  Gulf 
of  Mexico  and  lower  U.S.  east  coast  hurricanes.  Little  is  known  about  the 
characteristics  of  the  PBL  in  hurricanes  over  land.  However,  given  the  high 
level  of  turbulent  mixing  in  hurricanes,  it  is  reasonable  to  assume  that  an 
adiabatic  lapse  rate  is  established  and  that  at  least  in  the  near-coastal 
zone,  the  boundary  layer  depth  is  close  to  that  assumed  for  the  over-water 
case. 

29.  Given  the  above  conditions,  the  parametric  relations  F, ,  F2,  F,,  and 
F4  can  be  found  once  for  a  given  storm  by  iteration,  and  expressed  in  terms  of 
tables.  In  practice,  the  upwind  and  crosswind  drag  coefficients,  the  ratio 
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iu/|p| ,  and  the  angle  0  between  the  surface  wind  and  the  integrated  wind  are 
computed  and  stored  in  a  table  for  |tf|  values  ranging  from  0.8  to  80  m/sec  in 
increments  of  0.8  m/sec.  Values  for  intermediate  wind  speeds  are  found  by 
linear  interpolation. 

Terrain  Transformation  Model 


30.  In  this  section  a  simple  model  is  developed  to  derive  the  surface 
wind  and  stress  over  terrain  of  arbitrary  roughness  from  the  numerical  wind- 
field  solution  computed  exclusively  from  the  revised  over-water  drag  formula¬ 
tion.  The  approach  is  to  employ  equilibrium  PBL  theory  to  relate  the  over¬ 
water  integrated  PBL  wind  to  the  flow  at  the  top  of  the  PBL  and  then  to  employ 
a  consistent  equilibrium  model  to  compute  the  surface  wind  stress  from  the 
wind  at  the  top  of  the  PBL  for  terrain  (including  lakes)  of  arbitrarily  speci¬ 
fied  roughness.  The  model  proposed  assumes  that  the  PBL  over  land  or  inland 
lakes  in  a  hurricane  is  neutrally  stratified. 

31.  The  transformations  are  derived  quite  simply  from  consideration  of 
the  alternate  forms  of  the  similarity  PBL  theory  adapted  in  this  study. 

Surface  drag  in  the  numerical  vortex  model  was  parameterized  using 
Equations  13-15  and  13-16,  which  relate  the  stress  to  the  integrated  PBL  wind. 
Alternatively  (Arya  1977),  surface  drag  may  be  referenced  to  the  wind  at  the 
top  of  the  PBL  with  components  uh  and  vk: 

^  -  (In  z,  +  A)  (13-30) 


kv> 

u. 


m 


or  to  the  surface  geos trophic  wind  components 


^  -  dn  z.  +  A,) 


u. 


(13-31) 


(13-32) 

(13-33) 


As  noted  by  Arya,  A„  and  B0  may  be  expected  to  differ  from  A  and  B  due  to  the 
presence  of  baroclinicity  and  also  in  very  low  latitudes  where  winds  are 
strongly  geostrophic.  In  hurricanes,  baroclinicity  in  the  PBL  may  be  ignored, 
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but  the  flow  at  the  top  of  the  PBL  is  aore  nearly  in  gradient  balance.  The 
effects  of  curvature  on  A  and  B  have  not  been  studied,  but  the  success 
achieved  with  the  similarity  PBL  theory  in  the  overwater  case  suggests  that 
such  effects  are  not  large.  In  this  model,  differences  between  A  and  B  and  A* 
and  B„  are  ignored. 

32.  The  relationship  between  A„  and  B0  and  A,  and  Bm  is  given  by 
Arya  (1977)  as  derived  from  the  equations  of  mean  motion  for  a  barotropic 
atmosphere  in  which  the  momentum  flux  is  assumed  to  vanish  at  z  -  h: 


Am  “  A. 


(13-34) 


K 


(13-35) 


For  the  restricted  case  of  fh/u.  —  1,  which  was  adopted  only  for  the  purposes 
of  attaining  a  workable  parameterization,  it  can  be  shown  simply  from 
Equations  13-15  through  13-18,  13-30,  13-31,  13-34,  and  13-35  that  the  flow  at 
the  top  of  the  PBL,  z  -  h,  is  related  to  the  vertically  integrated  flow 
through 

u,  «  u  (13-36) 

vt  ■  v  -  u.  (13-37) 

In  the  coordinate  system  adopted  (see  Figure  13-2),  vh  is  negative  and  u.  is 
always  positive,  so  that  the  wind  speed  at  the  top  of  the  PBL  is  always  larger 
than  the  mean  layer  wind  speed  P  and  turned  clockwise  (in  the  Northern 
Hemisphere)  with  respect  to  the  mean  layer  wind  ^  . 

33.  Given  the  wind  at  level  h  ,  the  consistent  similarity  theory 
defined  by  Equations  13-30  and  13-31  may  be  solved  for  the  surface  stress  and 
surface  layer  wind  in  a  neutral  PBL  over  terrain  of  specified  roughness  z„ . 
For  a  neutrally  stratified  PBL,  A,  and  B,  are  reduced  simply  to  constants 
(1.39  and  1.95  +  k,  respectively).  If  z0  is  a  constant,  as,  say,  over  a 
homogeneous  land  surface,  u.  may  be  obtained  directly  from  Equations  13-30 
and  13-31.  However,  since  roughness  over  a  lake  probably  depends  on  u. ,  z„ 
may  in  general  be  prescribed  in  terms  of  u.  using  the  general  form  proposed 
by  Cardone  (1969): 
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2.  -  C,±  .  C,  a’.  .  C,  (13-38) 

where  Ct,  C2,  and  C2  are  constants  to  be  chosen  to  impose  a  desired  drag  law. 
(For  example,  for  a  Chamock  law,  C,  -  C,  -  0,  C2  -  a/g\  for  land,  C,  -  C2  -  0, 

C3  is  the  roughness  parameter  for  the  terrain  type) . 

34.  Implementation  of  the  above  model  in  the  specification  of  hurricane 
surface  wind  fields  over  terrain  of  arbitrary  roughness  or  lakes  is  coded  as 
subroutine  UPDOWN,  which  allows  for  the  calculation  of  transformations  for  up 
to  six  terrain  categories.  Roughness  constants  C, ,  C2,  and  <7,  are  specified 
for  each  terrain  category.  The  procedure  followed  is: 


fl.  Given  the  integrated  boundary  layer  wind  u,  v,  and  the 

conditions  of  a  given  hurricane  over  water  (h,  f,  k,  a), 

compute  u.  from  the  revised  overwater  similarity 
parameterization. 

fe.  From  Equations  13-36  and  13-37,  compute  the  wind  speed  and 
direction  at  level  h,  the  top  of  the  PBL. 

£.  For  each  terrain  roughness  specified  in  terms  of  Equation  13- 
38,  the  neutral  similarity  model  (13-30  and  13-31)  is  used  to 
compute  the  friction  velocity  appropriate  to  the  terrain 
roughness  iu,,  the  ratio  iuS\V\  ,  and  the  angle  between  the 
surface  stress  and  the  integrated  wind.  The  ratio  and  the 
turning  angle  are  computed  for  |^|  ranging  from  0.8-80.0  m/sec 
in  increments  of  0.8  m/sec,  for  each  roughness  category  and 
stored  for  use  in  the  specification  of  surface  winds  in  a  given 
simulation  over  a  grid  covering  different  terrain  types. 
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35.  The  overall  behavior  of  the  transformations  is  exemplified  in 
Figure  13-3,  which  shows  the  ratio  of  surface  wind  speeds  at  20  m  (overland 
divided  by  oversea)  and  the  difference  between  the  overland  and  oversea  inflow 
angle  for  two  terrain  roughnesses:  0.04  m  and  0.32  m.  For  comparison.  Figure 
13-3  shows  the  results  for  the  wind- speed  ratio  derived  from  numerical  mixed- 
terrain  and  overwater  solutions  for  a  symmetric  stationary  vortex  (radius  and 
pressure  drop  as  in  Camille).  To  arrive  at  the  indicated  quantities,  the 
surface  wind  speed  and  direction  along  a  radial  extending  north  of  the  eye 
over  land  in  the  mixed  terrain  case  were  referenced  to  the  (symmetrical)  solu¬ 
tion  along  the  radial  for  open  ocean.  It  should  be  recalled  that  in  the 
mixed- terrain  solution,  the  wind  field  over  land  looked  reasonable.  Appar¬ 
ently,  that  solution  can  be  retrieved  simply  from  the  overwater  solution  using 
the  equilibrium  model  described  above.  It  is  also  interesting  to  note  that 
the  form  of  the  dependence  shown  in  Figure  13-3  conforms  closely  to  the 
empirical  wind- speed  ratio  derived  from  measurements  in  hurricanes  in  and 
around  Lake  Okeechobee  by  the  U.S.  Weather  Bureau,  Hydrometeorological  Section 
(Cardone,  Greenwood,  and  Greenwood  1992). 


Figure  13-3.  Behavior  of  overland/oversea  transformation  (from  Cardone, 

Greenwood,  and  Greenwood  (1992)) 
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PART  III:  GENERAL  STRUCTURE  OF  HURWIN 


36.  HURWIN  was ‘prepared  to  simplify  usage  of  the  tropical  storm 
planetary  boundary  layer  routines.  The  program,  in  conjunction  with  CHS 
library  routines,  assists  with  construction  of  input  files,  building  SCRIPT 
files  (job  control  language  on  the  Cray  supercomputer)  to  run  the  programs  and 
manipulate  files,  and  producing  various  output  products.  HURWIN  is  a  branch 
of  the  CHS  Spectral  Wave  Modeling  Hodule.  The  name  HURWIN  represents  the 
whole  collection  of  codes  needed  to  generate  hurricane  wind  fields  (wind 
direction  and  either  wind  speed  or  wind  stress  defined  as  u.*)  with  the 
planetary  boundary  layer  model.  An  option  is  available  to  save  pressure 
fields  as  well.  The  terms  "tropical  storm,"  "hurricane,"  and  "typhoon"  are 
used  interchangeably.  A  hurricane  is  merely  a  tropical  storm  with  maximum 
sustained  wind  speeds  exceeding  75  mph.  A  typhoon  is  a  tropical  storm  located 
west  of  the  date  line  in  the  western  Pacific  Ocean.  The  present  version  of 
HURWIN  is  not  generalized  or  verified  for  use  with  tropical  storms  in  the 
Southern  Hemisphere. 

37.  HURWIN  is  based  on  the  concept  that  a  tropical  storm  generally 
changes  structure  relatively  slowly.  For  example,  the  structure  of  a 
hurricane  traveling  over  open  water  can  usually  be  well -represented  by 
parameters  specified  at  intervals  of  from  6  to  24  hr.  These  representative 
states  are  referred  to  as  "snapshots”  in  HURWIN  nomenclature.  Although  the 
structure  typically  changes  slowly,  the  storm  position  can  change  relatively 
quickly.  The  differences  in  time  scale  between  storm  structure  and  position 
changes  are  accommodated  in  HURWIN  by  a  very  flexible  procedure.  The  user  is 
allowed  to  specify  the  position  of  the  storm  center  at  every  hour.  The 
appropriate  snapshot  or  pair  of  snapshots  for  that  hour  is  also  specified. 
Snapshot(s)  are  assumed  to  be  coincident  with  the  actual  storm  center.  If  a 
pair  of  snapshots  is  given,  the  user  specifies  the  relative  weighting  to  be 
assigned  to  each  snapshot  and  the  wind  field  is  linearly  interpolated  between 
the  two  snapshots.  The  user  has  additional  flexibility  through  an  option  to 
rotate  snapshots  through  any  angle.  This  option  changes  the  direction  of 
storm  approach. 

38.  The  operational  functions  of  the  HURWIN  branch  of  the  Spectral  Wave 
Hodeling  module  are  illustrated  in  Figure  13-4.  The  three  main  functions  are: 
input  file  building,  running  programs,  and  post-processing.  In  general,  the 
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Figure  13-4.  Operational  functions  of  HURU1N 


file -building  function  Bust  be  used  prior  to  the  functions  to  run  progress  and 
post-process  the  aodel  results.  The  routines  called  under  each  function  are 
represented  in  Figure  13-4.  Those  with  naaes  ending  in  ”.F"  are  written  in 
FORTRAN.  The  others  are  written  in  SCRIPT.  The  aain  purpose  of  each  routine 
is  as  follows: 

-  CMS_HURWIN_MN:  SCRIPT  file  to  select  the  desired 

function  (build  input  files,  run,  or  post-process)  and  run  the 
appropriate  SCRIPT  file. 

File-Building  Function 

CMS_HURWIN_FB :  SCRIPT  file  to  run  input  file-building  routine. 

-  HURINTERACT . F :  FORTRAN  prograa  to  create  input  data 

files  or  specify  input  files  that  are  already  available, 
including:  1)  namelists  (described  in  a  following  section); 

2)  output  grid  information;  3)  output  station  information; 
and  4)  hourly  time  history  information  on  the  storm. 

Run  Programs  Function 

-  CMS_HURWIN_RP:  SCRIPT  file  to  run  HURJ0BER1.F  and 
hurricane  wind  computation  programs  (HURWIN_RUN) . 

-  HURJ0BER1.F:  FORTRAN  program  to  create  SCRIPT  file 

to  run  winds,  including:  1)  create  HURWIN_RUN;  2)  prompt 
user  for  output  file  name  header;  and  3)  write  output  file 
names  to  a  plotting  file. 

-  HURWIN_RUN:  SCRIPT  file  to  run  hurricane  wind  computation 
programs . 

SNAP.F:  FORTRAN  program  to  compute  snapshot  wind  fields  at 
specified  times.  (Note:  There  are  two  different  programs  here. 
HURWIN  automatically  invokes  the  proper  program,  depending  on 
the  parameters  specified  by  the  user.  The  program  file 
"snap_adc.f"  is  used  for  the  standard  pressure  profile;  the 
program  file  "snap_hol.f"  is  used  for  the  generalized  pressure 
profile. ) 

HIST.F:  FORTRAN  program  to  interpolate  wind  fields 

in  space  and  time  to  provide  hourly  results  at  selected  stations 

and/or  output  grid  locations. 

BLEND. F:  FORTRAN  program  to  blend  wind  fields 
created  during  coarse  (background)  and  fine  resolution  runs  of 
SNAP.F  and  HIST.F  to  provide  winds  on  the  user-specified  output 
grid.  The  blended  wind  field  consists  of  winds  from  the  fine 
grid  over  its  coverage  area  and  winds  from  the  coarse  grid  (with 
25-km  cell  size)  for  other  output  grid  areas.  A  nearly 
identical  program  "blendp.f"  performs  the  same  function  for 
pressure  fields. 
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-  CMS_HURWIN_PP :  SCRIPT  file  to  run  post-processing  programs. 

-  VINDPLOT.F:  FORTRAN  program  to  plot  wind  fields. 

-  TIMEPLOT.F:  FORTRAN  program  to  plot  time  history  of 
winds  at  selected  stations  or  output  grid  locations. 

39.  Program  HURVIN  uses  a  number  of  files  as  inputs  or  outputs  for 
various  routines.  Some  of  the  files  may  be  pre-existing  files  provided  as 
input  by  the  user,  but  most  are  created  in  HURVIN.  Many  files  are  created  and 
used  internally  as  part  of  the  HURVIN  stream  of  programs,  and  are  not  of 
direct  interest-  to  the  user.  The  user  has  some  freedom  in  specifying  names  of 
external  input  and  output  files.  All  internal  files  are  pre-named.  Figures 
13-5  through  13-7  indicate  the  files  created  and  used  by  each  program  in 
HURVIN.  Each  file  has  associated  option  and  unit  numbers.  Unit  numbers  are 
arbitrary.  They  are  important  in  understanding  program  operation,  but  they 
are  irrelevant  to  a  routine  user.  Required  input  files  are  identified  with 
the  letter  "R."  Files  that  are  automatically  created  by  routines  are 
identified  with  the  letter  "A."  Files  that  the  user  has  an  option  of 
providing  or  not  providing  (though  they  may  be  required  later  in  the  program 
stream)  are  identified  with  the  letter  "0."  The' contents  and  format  of  data 
within  each  file  are  described  in  Part  V.  Coordination  between  programs  in 
the  three  functions  is  evident  in  that  files  created  in  one  function  are  often 
used  as  input  in  another  function. 

40.  The  organizational  procedure  for  a  typical  simulation  is  as 
follows : 

-  Run  the  input  file -building  routine  to  provide  input 
data  on  the  simulation  and  desired  output. 

-  Using  the  input  data,  run  the  computational  programs 

to  generate  wind  fields  at  desired  stations  and/or  on  a  desired 
output  grid. 

Plot  time  history  of  wind  parameters  at  selected 
stations  and/or  plot  full  wind  field  at  selected  times. 

-  Use  output  wind  fields  as  input  to  a  wave,  storm 

surge,  or  circulation  model  (e.g.,  the  SHALW  wave  branch  of 
the  Spectral  Vave  Modeling  Module) . 
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Figure  13-6.  Run  programs  function 
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PART  IV:  HURWIN  MODULE  FUNCTIONS  AND  ROUTINES 


41.  The  first  function  of  the  HURWIN  branch  of  the  Spectral  Wave 
Modeling  module  of  CMS  is  to  build  most  or  all  of  the  necessary  input  files 
for  the  user's  application.  The  second  function  is  to  run  the  computational 
programs.  The  third  function  is  to  display  graphically  or  statistically  the 
results  of  the  user's  most  recent  HURWIN  run.  Descriptions  of  the  operation 
of  each  function  and  routine  and  the  specific  information  required  by  each 
function  and  routine  are  presented  in  the  following  sections. 

File-Building  Function 

42.  The  first  function  of  the  HURWIN  module  is  to  assist  users  in 
building  the  input  files  required.  Users  are  prompted  for  input  through  a 
series  of  menu  selections  and  quest ion -answer  sessions.  It  is  recommended 
that  the  file-building  function  routines  be  used  at  all  times  to  organize  data 
and  files  for  a  HURWIN  run.  This  precaution  will  minimize  potential  file 
format  and  organizational  errors. 

43.  The  HURWIN  file-building  function  includes  options  to  use  any  pre¬ 
existing  input  files.  File  names  of  pre-existing  files  must  be  different  from 
the  file  names  used  in  HURWIN.  For  example,  a  pre-existing  file  must  not  be 
named  ' snapio.hur . '  The  user  can  specify  which  files  are  already  complete  and 
which  files  need  to  be  built.  Then  HURWIN  will  assist  the  user  in  building 
the  necessary  files.  It  will  also  show  the  user  the  contents  of  pre-existing 
files,  give  the  user  an  opportunity  to  change  any  values,  and  write  the  pre¬ 
existing  information,  with  any  user-specified  changes,  to  the  appropriately 
named  HURWIN  input  files. 

44.  The  input  file -building  routine  is  used  to  construct  the  required 
files  of  ..*>melists  1,  2,  and  3  (snapin.hur) ;  namelists  4  and  5,  station  and 
history  parameters  (histin.hur) ;  parameter  file  (filel.hur);  main  file  of  file 
names  (filenmi.hur) ;  output  grid  parameters  (gridparm.hur) ;  and  land/sea  table 
(lndsea.hur) .  The  latter  two  files  are  optional,  needed  only  when  an  output 
grid  of  winds  is  desired. 

45.  Descriptions  of  the  data  requested  by  the  input  file-building 
routine  (HURINTERACT. F)  are  provided  below.  The  descriptions  have  the 
following  format:  the  prompts  that  the  user  will  see  during  the  interactive 
file-building  session  are  provided  by  "Question:"  as  they  would  appear  on  the 


13-24 


user's  terminal  during  the  interactive  file-building  session  along  with  an 
identifying  number  for  reference  (e.g.  GF-1,  which  stands  for  general  £ile 
Building  question  number  1).  A  list  of  "Options:"  is  provided,  which  may  be 
used  to  answer  the  question.  If  the  question  allows  only  a  limited  range  of 
responses,  then  the  permitted  range  of  values  is  given.  Some  questions  allow 
the  user  to  choose  preset  responses,  defined  in  "Default:".  Also  the  FORTRAN 
name  of  the  variable  being  specified  is  given  in  parentheses  on  the  "Default:" 
line.  "Preselect:”  identifies  answers  to  previous  questions  that  are  required 
in  order  for  the  user  to  encounter  the  current  question.  For  example, 
question  GF-7  will  be  asked  only  if  the  user  previously  answered  "Yes"  to 
question  GF-6.  Each  question  is  explained  under  the  "Define:”  section,  and  a 
list  of  "Cautions : "  is  provided  when  appropriate  to  alert  users  to  special 
requirements  of  the  module.  Vhen  necessary,  a  "PATH:"  is  provided  that 
indicates  which  question  is  next  based  on  answers  to  previous  questions. 

While  all  of  the  possible  questions  are  listed  below,  the  user  will  not 
encounter  all  of  them  in  any  simulation. 


GF-1  Question 
Options 


Default 

Preselect 

Define 

Caution 

GF-2  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-3  Question 


Enter  DESIRED  SYSTEM  OF  UNITS  > 

1  -  METRIC  (m  &  km,  s,  m/s,  Deg-C) 

2  -  ENGLISH1  (ft  &  miles,  s,  mph,  Deg-F) 

3  -  ENGLISH2  (ft  &  nm,  s,  knots,  Deg-F) 

4  -  OCEANWEA  (units  of  Cardone  et  al.  1992) 

4  (IUNITS) 

None . 

Select  the  system  of  units  to  be  used  for  user-specified 
input . 

None. 

Enter  FOUR  CHARACTER  NAME  OF  STORM  > 

Any  group  of  up  to  four  characters  or  numbers. 

None.  (STORM) 

None . 

Specify  an  identifying  name  for  the  storm. 

None. 

Enter  NUMBER  OF  SNAPSHOTS  > 
value  >  0 

None .  (NZ) 

None. 


Options 

Default 

Preselect 


1  j-  25 


GF-4 


GF-5 


GF-6 


PATH: 

GF-7 


GF-8 


Define 

Caution 

Question 

Options 


Default  : 
Preselect  : 
Define  : 
Caution  : 
Question 
Options  : 
Default 
Preselect  : 
Define  : 

Caution  : 
Question  : 
Options  : 

Default  : 
Preselect  : 
Define  : 
Caution  : 
o  IF  "N" 
Question  : 
Options 
Default 
Preselect  : 
Define 

Caution  : 

Question 

Options 


Specify  the  number  of  snapshots  to  be  read. 

None. 

Enter  TIME  ZONE  FOR  STARTING  TIME  > 

1  -  (UTC)  Universal  Time  Coordinate  (formerly  GMT) 

2  -  (EST)  Eastern  Standard  Time 

3  -  (EDT)  Eastern  Daylight  Time 

4  -  (CST)  Central  Standard  Time 

5  -  (CDT)  Central  Daylight  Time 

1  (IZONE) 

None. 

Select  time  zone  for  input  and  output. 

None . 

Enter  STARTING  TIME  (YYMMDDHH)  > 

None.  (I START) 

None . 

Enter  the  year,  month,  day,  and  hour  for  the  first 
simulation  time. 

None. 

Change  parameters  in  roughness  law?  > 

Y  -  yes,  change  at  least  one  parameter; 

N  -  no,  use  default  values. 

None.  (REPLY) 

None. 

Provides  an  option  to  change  roughness  law  parameters. 
None. 

ON  GF-6,  GO  TO  GF-14. 

Enter  AIR-LAND  TEMP  DIFFERENCE  (Deg-C,Deg-F)  > 

-12.0  ^  value  <  12.0 

G  0  (DTHl) 

GF-6,  option  "Y". 

Specify  temperature  difference  between  air  and  land. 

(+  -  land  temperature  is  lower,  -  -  land  temperature  is 
higher) 

None. 

Enter  AIR-SEA  TEMP  DIFFERENCE  (Deg-C,Deg-F)  > 

-12.0  5  value  <  12.0. 
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(DTH2) 


Default 

Preselect 

Define 


Caution 

GF-9  Question 
Options 

Default 

Preselect 

Define 

Caution 

GF-10  Question 
Options 


Default 

Preselect 

Define 

Caution 

GF-11  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-12  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-13  Question 
Options 


-2.0 

GF-6,  option  "Y" . 

Specify  temperature  difference  between  air  and  water. 

(+  -  sea  temperature  is  lower,  -  -  sea  temperature  is 
higher) . 

None . 

Enter  BOUNDARY  LAYER  HEIGHT  OVER  VATER  (rn.ft)  > 

0.0  <  value  S  2000.0  m; 

0.0  <  value  i  6600.0  ft. 

500.0  m;  1640.4  ft.  (HH) 

GF-6,  option  "Y" . 

Specify  the  height  of  the  planetary  boundary  layer  over 
water . 

The  default  value  should  be  used  unless  the  user  has 
specific  knowledge  about  the  storm  being  modeled. 

Enter  POTENTIAL  TEMPERATURE  (deg  C,  deg  F,  deg  K)  > 

0.0  <  value  <  47.0  deg  C 

32.0  <  value  <  116.6  deg  F 

280  <  value  <  320  deg  K 

26.85  Deg-C;  80.33  deg  F;  300  deg  K.  (PTH) 

GF-6,  option  "YM. 

Specify  the  potential  temperature. 

The  default  value  should  be  used  unless  the  user  .has 
specific  knowledge  about  the  storm. 

Enter  VON  KARMAN'S  CONSTANT  > 

0 . 3  <  value  S  0.5 

0.35  (K35) 

GF-6,  option  "Y" . 

Specify  value  of  von  Karman's  constant. 

The  default  value  is  recommended. 

Enter  ROUGHNESS  LENGTH  OVER  LAND  > 

0.0  <  value  <  1.0 

0.08  (Z0LAND) 

GF-6,  option  "Y" . 

Specify  the  characteristic  roughness  length  over  land. 
None . 

Enter  CHARNOCK'S  CONSTANT  > 

0.01  <  value  <  0.05. 
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(GARR) 


GF-14 


PATH: 

GF-15 


GF-16 


GF-17 


Default 

Preselect 

Define 

Caution 

Question 

Options 

Default 


Preselect  : 
Define  : 


Caution 
o  IF  "N" 
Question  : 

Options  : 

Default  : 

Preselect  : 

Define  : 


Caution 


0.035 

GF-6 ,  option  "Y". 

Specify  Chamoclc's  constant. 

The  default  value  is  recommended. 

Do  you  want  to  change  any  of  the  Z0  coeff.T  > 

Y  -  yes,  change  at  least  one  set  of  coeff . ; 

N  -  no,  use  default  values. 

(ZCOEFF) 


Terrain 

£i 

O2 

& 

1)  Lake 

0.0 

0.0037 

0.015 

2)  Marsh 

0.0 

0.0 

0.04 

3)  Plains 

0.0 

0.0 

0.16 

4)  Woods 

0.0 

0.0 

0.32 

5)  Cities 

0.0 

0.0 

1.28 

None. 

Provides  an  option  to  change  roughness  law  coefficients 
for  non-open  water  terrains  (Equation  13-25).  Note  that 
the  coefficients  are  dimensioned  in  metric  units. 

None. 

ON  GF-14,  GO  TO  GF-19. 

Correct  which  number  ?  (1-5)  > 

1  S  value  £  5 
None. 

GF-14,  option  "Y". 

Specify  terrain  for  which  new  roughness 
coefficients  are  desired.  Terrain  type 
identification  numbers  are  given  with  GF-14. 

None. 


Question 

Options 

Default 

Preselect 

Define 

Caution 

Question 

Options 


Enter  value  Cl: 


None. 

GF-14,  option  "Y". 

Specify  value  for  C,  in  roughness  law 
(Equation  13-25). 

There  are  no  internal  checks  to  ensure  that 
reasonable  values  are  entered. 

Enter  value  C2: 


13-28 


Default 

Preselect 

Define 

Caution 

GF-18  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-19  Question 
Options 

Default 

Preselect 

Define 

Caution 

GF-20  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-21  Question 
Options 


Default 

Preselect 

Define 


None. 

GF-14,  option  *Y* . 

Specify  value  for  Cj  in  roughness  law 
(Equation  13-25). 

There  are  no  internal  checks  to  ensure  that 
reasonable  values  are  entered. 

Enter  value  C3: 

None . 

GF-14,  option  "Y"  . 

Specify  value  for  C}  in  roughness  law 
(Equation  13-25). 

There  are  no  internal  checks  to  ensure  that 
reasonable  values  are  entered. 

Enter  CHOICE  FOR  PRESSURE  SPECIFICATION  > 

1  -  standard  (single  exponential) 

2  -  generalized  (double  exponential) 

1 

None. 

Specify  form  of  pressure  profile. 

None. 

Enter  NUMBER  OF  COMPUTATIONAL  NESTS  (5-7)  > 

5  £  value  <  7 
7 

None. 

Specify  number  of  grid  nests  to  be  used  in 
computing  the  vortex  wind  fields. 

None . 

Enter  INNERMOST  GRID  SPACING  (km. miles ,nm)  > 

0.0  <  value  <20.0  km 

0.0  <  value  <12.4  miles 

0.0  <  value  <10.8  nm 

2.0  km;  1.2  miles;  1.0  nm. 

None. 

Specify  spacing  between  points  in  the 
innermost  (and  finest)  nested  grid  to  be  used 
for  computation.  DXX  is  converted  to  the 
parameter  DX  described  by  Cardone  et  al.  (1994) 


(IPSPEC) 


(NESTS) 


(DXX) 


Caution 


None . 


GF-22  Question 
Options 
Default 

Preselect 

Define 

Caution 

GF-23  Question 
Options 

Default 

Preselect 

Define 

Caution 


Enter  NUMBER  07  WIND  COMPUTATION  CYCLES  > 

0  <  value  <  4,000 

800  (GF-19 ,  option  "l"); 

3200  (GF-19,  option  "2") 

None. 

Specify  number  of  iterations  to  be  used  in 
computing  the  vortex  wind  fields. 

Default  values  are  recommended. 

Enter  PRESSURE  FIELD  QUADRANT  INDICATOR  > 

0  -  circularly  symmetric  pressure  field. 

1  -  first  quadrant  is  right  front. 

2  -  first  quadrant  is  forward. 

0 

GF-19,  option  "1". 

Specify  whether  the  pressure  field  is  to  be 
circularly  symmetric  or  not.  If  not,  specify 
the  location  of  the  first  quadrant. 

None. 


Questions  GF-24  through  GF-35  are  repeated 
for  each  snapshot. 


GF-24  Quest!  ^ 
Options 

Default 

Preselect 

Define 

Caution 

GF-25  Question 
Options 
Default 
Preselect 
Define 


Enter  SURFACE  GE0STR0PHIC  WIND , SGV 
(m/s , mph , knots  )> 

0.0  <  value  <50.0  m/sec. 

0 . 0  <  value  <  111.85  mph. 

0.0  <  value  <97.2  knots. 

None . 

None. 

Specify  the  speed  of  the  surface  geostrophic 
wind  or  steering  flow  (i.e.,  the  wind  which 
would  be  present  if  the  hurricane  were  not 
there) . 

None. 

Enter  DIRECTION  OF  SGW  (Deg)  > 

0.0  ^  value  <  360.0 
None. 

None. 

Specify  the  direction  from  which  the  surface 


(NM) 


(IQUAD) 


(SGW) 


(AN1) 
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Caution 


GF-26  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-27  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-28  Question 
Options 

Default 

Preselect 

Define 

Caution 

GF-29  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-30  Question 
Options 


geostrophic  wind  is  blowing,  in  deg  azimuth 
(meteorological  convention) . 

Some  old  parameter  files  may  be  based  on  other 
conventions,  such  as  direction  toward  which  the  wind  is 
blowing  and  deg  measured  counterclockwise  from  the 
x-axis  (Cardone,  Greenwood,  and  Greenwood  1992). 

Enter  NORTH  LATITUDE  OF  EYE  OF  STORM  (Deg)  > 

0.0  <  value  <  90.0 

None.  ( EYELAT ) 

None . 

Specify  latitude  coordinate  of  storm  eye 
(Northern  Hemisphere) . 

The  value  entered  should  best  characterize  the 
snapshot  (even  though  it  may  actually  be 
changing) . 

Enter  TRACK  DIRECTION,  clockwise  from  north  > 

0.0  ^  value  <  360.0 

None.  (DIREC) 

None . 

Specify  direction  toward  which  storm  is  moving, 
in  deg  azimuth. 

None. 

Enter  F0RVARD  SPEED  OF  STORM  (m/s ,mph, knots)  > 

0.0  <  value  <  50.0  m/s 
0.0  <  value  <  111 . 85  mph 
0.0  <  value  <97.2  knots 

None.  (SPEED) 

None. 

Specify  forward  translation  speed  of  storm. 

None. 

Enter  PRESSURE  AT  EYE  OF  STORM  (mb)  > 

850.0  5  value  <  1,000.0  mb 

None.  (EYPRES) 

None. 

Specify  central  pressure  of  storm. 

None . 

Enter  PRESSURE  PROFILE  SCALE 
RADIUS(km, miles ,nm)  > 

0.0  <  value  <  200.0  km 


13-31 


0.0  <  value  i  12S.0  miles 
0.0  <  value  :S  125.0  nm 

Default  None.  (RADIUS) 


Preselect  :  GF-19.  option  '1'  and  GF-23.  If  GF-23,  option  '0' 

is  chosen,  only  1  value  for  RADIUS  is  needed 
(other  3  default  to  0) .  If  '1'  or  '2'  option  is 
chosen,  four  values  for  RADIUS  are  needed. 

GF-19.  option  *2*.  Two  values  for  RADIUS  are  needed. 


Define  Specify  scaling  radius  (often  equivalent  to 

"radius  to  maximum  wind”)  for  storm. 


Caution  None . 

GF-31  Question  Enter  FAR  FIELD  PRESSURE  (mb)  > 

Options  950.0  <  value  <  1,050.0  mb 

Default  :  None.  (PFAR) 

Preselect  None . 

Define  :  Specify  atmospheric  pressure  outside  the 

hurricane  area. 


Caution  None . 


GF-32  Question 

Options 

Default 

Preselect 

Define 

Caution 

GF-33  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-34  Question 
Options 
Default 
Preselect 
Define 

Caution 


Enter  DISTANCE  FROM  AXIS  TO  0.5  SGff 
(km,  miles,  nm)> 

0  <,  value 

0.0  (ST12) 

GF-19,  option  »1* . 

The  default  value  is  recommended. 

Enter  EXPONENT  FOR  INNER  PRESSURE  RADIUS  > 

0.0  <  value  :£  3.0 

T.O  (Bl) 

GF-19,  option  '2'. 

Specify  Holland's  (1980)  B  coefficient  for  the 
inner  pressure  component. 

None. 

Enter  EXPONENT  FOR  OUTER  PRESSURE  RADIUS  >0.0 
0.0  £  value  ^  3.0 

0.0  (B2) 

GF-19,  option  '2'. 

Specify  Holland's  (1980)  B  coefficient  for  the 
outer  pressure  component. 

Automatically  set  to  0.0  if  DPI  not  specified 
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# 


GF-3S  Question 
Options 
Default 
Preselect 
Define 

Caution 


(GF-35)  or  if  DP1-PFAR- EYPRES .  If  B2-0.0,  the 
second  value  of  RADIUS  (GF-30)  is  set  to  0.0. 

Enter  PRESS.  DIFF.  BETS.  E7E  &  INNER  RADIUS  (mb)  > 

Q.O  S  value  3  (PFAR  -  EYPRES) 

PFAR  -  EYPRES  (DPI) 

GF-19,  option  '2'. 

Specify  pressure  anomaly  for  the  inner  pressure 
component . 

Automatically  set  to  PFAR  -  EYPRES  if  B2  -  0.0 
(GF-34) . 


*****  Questions  GF-24  through  GF-35  are  repeated 
*****  for  each  snapshot. 


GF-36 

Question 

Enter  OUTPUT  FIELD  OPTION  > 

Options 

1  -  surface  wind  speed  and  direction. 

2  -  surface  wind  stress. 

Default 

1  (IFIELD) 

Preselect 

None. 

Define 

Specify  the  type  of  output  desired.  The  wind 
stress  option  gives  u„2  as  output  instead  of 
wind  speed. 

Caution 

None. 

GF-37 

Question 

Enter  DESIRED  WIND  ELEVATION  (a, ft)  > 

Options 

0.0  <  value  S  500.0 

Default 

10.0  m;  33.0  ft  (GRIDHT) 

Preselect 

GF-36,  option  "1". 

Define 

Specify  the  elevation  above  the  water  at  which  wind 
output  is  desired. 

Caution 

None. 

GF-38 

Question 

Enter  DESIRED  OUTPUT  UNITS  >  2 

Options 

1  —  a/sec 

2  -  knots 

Default 

2  (KSPD) 

Preselect 

GF-36,  option  "1". 

Define 

Specify  the  desired  units  for  output  wind  speed. 

Caution 

None. 

GF-39 

Question 

Enter  DESIRED  DIRECTION  OUTPUT  UNITS  >  2 

Options 

1  -  meteorological  convention  (deg  azimuth,  coming  from) 
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GF-40 


GF-41 


GF-42 


PATH: 

GF-43 


Default 

Preselect 

Define 

Caution 

Question 

Options 

Default 

Preselect 

Define 

Caution 

Question 

Options 


Default 

Preselect 

Define 

Caution 

Question 

Options 

Default 

Preselect 

Define 

Caution 

o  IF  "N" 

Question 

Options 

Default 

Preselect 

Define 


2  -  WIS  convention  (Figure  13-8) 

2  (KDIR) 

None. 

Specify  the  desired  units  for  output  wind  direction. 

None . 

Enter  DESIRED  OUTPUT  OPTION  >  0 
0  —  suppress  pressure  field  output 
1  —  create  pressure  field  output 

0  (PRSOUT) 

None. 

Specify  whether  an  output  file  of  pressure  fields 
is  desired. 

None. 

Enter  DESIRED  OPTION  FOR  PRINTING  > 

0  -  suppress  printing 

1  -  print  pressures  and  initial  winds  (first  snapshot 
only) 

0  (IB) 

None. 

Specify  whether  or  not  detailed  printouts  are 
desired  for  the  first  snapshot. 

The  default  value  is  recommended . 

Do  you  want  winds  interpolated  to  output  grid?  > 

Y  -  yes,  interpolate  winds  to  output  grid; 

N  -  no,  do  not  use  output  grid. 

Y  (INTERP) 
None. 

Specify  whether  winds  are  to  be  interpolated 
onto  an  output  grid. 

None. 

ON  GF-42,  GO  TO  GF-52. 

Enter  INTERVAL  (hrs)  TO  PRINT  WINDS  ON  GRID  > 

0  £  value 

None.  (NPRT) 

GF-42,  option  Y. 

Specify  time  interval  in  hours  at  which 
printed  winds  are  desired.  NPRT  equal  to 
0  suppresses  printing. 
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Caution 


GF-44  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-45  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-46  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-47  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-48  Question 
Options 
Default 
Preselect 
Define 

Caution 


If  NPRT  equal  to  0,  printing  of  selected 
windfields  can  still  be  done  using  NRPTGRID 
(GF-69) . 

Enter  LATITUDE  FOR  NORTH  EDGE  OF  GRID  (Deg)  > 

0  s  value  <90.0 

None.  (LATN) 

GF-42,  option  Y. 

Specify  latitude  boundary  of  output  grid. 

None. 

Enter  LATITUDE  FOR  SOUTH  EDGE  OF  GRID  (Deg)  > 

0  <  value  £  90.0 

None.  (LATS) 

GF-42,  option  Y. 

Specify  latitude  boundary  of  output  grid. 

LATS  cannot  be  greater  than  LATN;  grid  must  be 
in  Northern  Hemisphere. 

Enter  LONGITUDE  FOR  VEST  EDGE  OF  GRID  (Deg)  > 

0  <  value  ^  360.0 

None.  (L0NW) 

GF-42,  option  Y. 

Specify  longitude  boundary  of  output  grid. 

If  the  grid  falls  west  of  the  International 
Date  Line,  longitudes  must  be  increased  above 
180  deg.  Negative  longitude  values  are  not 
permitted.  LONW  cannot  be  less  than  LONE. 

Enter  LONGITUDE  FOR  EAST  EDGE  OF  GRID  (Deg)  > 

0  <  value  ^  360.0 

None.  (LONE) 

GF-42,  option  Y. 

Specify  longitude  boundary  of  output  grid. 

LONE  cannot  be  greater  than  LONW. 

Enter  GRID  CELL  SPACING  (Deg)  > 
value  >0.0 

None.  (DXOUT) 

GF-42,  option  Y. 

Specify  spacing  between  points  on  the  output 
grid. 

Spacing  must  be  the  same  (in  terms  of  degrees)  in 
both  the  latitude  and  longitude  directions. 
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GF-49 


GF-50 


GF-51 


GF-52 


PATH: 


Question 

Options 

Default 

Preselect 

Define 

Caution 

Question 

Options 


Default 

Preselect 

Define 

Caution 

Question 

Options 

Default 

Preselect  • 

Define 


Caution 
Question  : 

Options  : 

Default  : 

Preselect 
Define  : 

Caution  : 

o  IF  VALUE 


Enter  the  LAND/VATER  MATRIX  OPTION  > 

1  -  enter  the  matrix  interactively. 

2  -  set  constant  land/water  par  ante  ter. 

None.  (IOPTION) 

GF-42 ,  option  Y. 

Specify  whether  or  not  the  land/water  matrix 
values  are  constant. 

None. 

Enter  TERRAIN  CODE  > 

1  -  open  water 

2  -  lake 

3  -  marsh 

4  -  plains 

5  -  woods 

6  -  city 

1  ( LANDS EA) 

GF-42,  option  Y;  GF-49,  option  2. 

Specify  a  single  terrain  code  to  be  used  for 
the  entire  land/water  matrix. 

None. 

Enter  TERRAIN  CODE  FOR  LAND/VATER  MATRIX  > 

Same  as  for  GF-50. 

None.  (LSTAB(NX.NY) ) 

GF-42,  option  Y;  GF-49,  option  1. 

Specify  the  terrain  code  for  each  point  in  the 
land/water  matrix.  The  user  is  prompted  for 
each  value. 

None . 

Enter  NUMBER  OF  MEASUREMENT  STATIONS  > 

0  <  value  :£  100 

None.  (NOGAGES) 

None. 

Specify  the  number  of  measurement  stations. 

None. 

-  0  ON  GF-52,  GO  TO  GF-60. 


Questions  GF-53  through  GF-59  are  repeated 
for  each  station. 
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4 


GF-S3  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-54  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-55  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-56  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-57  Question 
Options 


Enter  DEGREES  OF  NORTH  LATITUDE  > 

0  sS  value  <  90 
None . 

GF-52,  value  >  0 

Specify  station  latitude  coordinate. 
None . 

Enter  MINUTES  OF  LATITUDE  > 

0  ^  value  <  60 
None. 

GF-52,  value  >  0 

Specify  station  latitude  coordinate. 
None . 

Enter  DEGREES  OF  VEST  LONGITUDE  > 

0  ^  value  <  360 
None. 

GF-52,  value  >  0 

Specify  station  longitude  coordinate. 
None. 

Enter  MINUTES  OF  LONGITUDE  > 

0  <  value  <  60 
None. 

GF-52,  value  >  0 

Specify  station  longitude  coordinate. 
None. 

Enter  TERRAIN  CODE  FOR  STATION  > 

1  -  open  water 

2  -  lake  (continued) 


(NORLATD) 


(NORLATM) 


(WESLOND) 


(WESLONM) 
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3  -  marsh 


GF-58 


GF-59 


***** 

***** 

GF-60 


Default 

4  -  plains 

5  -  woods 

6  -  city 

1 

(ITERRAIN) 

Preselect 

GF-52 ,  value  >  0 

Define 

Specify  terrain  code  for  station. 

Caution  : 

None . 

Question 

Enter  STATION  HEIGHT  IN  METERS  > 

Options  : 

0.0  :S  value  s  50.0 

Default 

10.0 

(GAGEHGT) 

Preselect  : 

GF-52,  value  >  0 

Define 

Specify  elevation  of  station. 

Caution  : 

None . 

Question  : 

Enter  STATION  NUMBER  > 

Options  : 

0  s  value  5  999 

Default 

None. 

(IDGAGE) 

Preselect  : 

GF-52,  value  >  0 

Define  : 

Specify  station  identification  number. 

Caution 

None. 

Questions  GF-53  through  GF-59  are  repeated 


for  each  station. 


Question 

Enter  NUMBER  OF  HOURS  FOR  HISTORY  > 

Options 

0  <  value  <  300 

Default  : 

None. 

(NOHRS) 

Preselect  : 

None. 

Define 

Specify  number  of  hours  for  time -history  of 
storm. 
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4 


Caution  None. 


Question*  GF-61  through  GF-69 
for  each  history  hour. 


tpeated 


GF-61  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-62  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-63  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-64  Question 
Options 
Default 
Preselect 
Define 
Caution 

GF-65  Question 
Options 
Default 
Preselect 
Define 

Caution 


Enter  DEG  OF  NORTH  LATITUDE  OF  EYE  OF  STOR£> 

0  s  value  <  90 

None. 

None. 

Specify  latitude  coordinate  of  storm  eye 
(Northern  Hemisphere). 

None . 

Enter  MIN  OF  NORTH  LATITUDE  OF  EYE  OF  ST0RM> 

0  <  value  <  60 
None . 

None . 

Specify  latitude  coordinate  of  storm  eye. 

None . 

Enter  DEG  OF  WEST  LONGITUDE  OF  EYE  OF  STORM> 

0  <  value  <  360 

None. 

None . 

Specify  longitude  coordinate  of  storm  eye 
(Western  Hemisphere) . 

None. 

Enter  MIN  OF  VEST  LONGITUDE  OF  EYE  OF  ST0RM> 

0  ^  value  <  60 

None. 

None . 

Specify  longitude  coordinate  of  storm  eye. 
None . 

Enter  SEQ.  NO.  OF  1ST  SNAPSHOT  WIND  FIELD  > 

1  <,  value  <  NZ 
None . 

None . 

Specify  sequence  number  of  first  snapshot  to 
be  used  in  interpolating  winds  for  this  hour. 

None . 


(HISNLATD) 


(HISN1ATM) 


(HISWLOND) 


(HISWLONM) 


(N01SNAP) 


13-39 


GF-66  Question 
Options 
Default 
Preselect 
Define 

Caution 

GF-67  Question 
Options 
Default 
Preselect 

Define 


Caution 

GF-68  Question 
Options 
Default 
Preselect 
Define 


Caution 


Enter  SEQ.  NO.  OF  2ND  SNAPSHOT  WIND  FIELD  > 

NOISNAP  £  value  5  NZ;  or  value  -  0 

0  (N02SNAP) 

None . 

Specify  sequence  number  of  second  snapshot  to 
be  used  in  interpolating  winds  for  this  hour. 

Normally  N02SNAP  is  equal  to  NOISNAP  +  1. 

None. 

Enter  INTERPOLATION  DISTANCE  (1ST, 2ND  SNAPSHOTS)> 

0.0  <  value  <  1.0 

None.  (DISTANCE) 

GF-3,  value  >  1 
GF-66,  value  >  NOISNAP 

Specify  the  proportionate  weighting  to  be 
given  to  the  second  snapshot  wind  field  during 
interpolation.  A  weighting  of  1  -  DISTANCE 
will  be  given  to  the  first  snapshot  wind 
field. 

None. 

Enter  CLOCKWISE  ROTATION  OF  SNAPSHOT  (Deg)  > 

0  <  value  <  360 

0  (IROTATE) 

None. 

Specify  the  clockwise  rotation  to  be  given  to 
the  snapshot.  When  this  parameter  is  not 
equal  to  zero,  the  storm  direction  will  be 
rotated.  The  default  value  (no  rotation)  is 
recommended. 

None . 


GF-69  Question 
Options 

Default 

Preselect 


Enter  DESIRED  OPTION  FOR  PRINTING  > 

0  -  no  printing  of  output  grid  wind  field 
1  -  print  output  grid  wind  field 

0  (NPRTGRID) 

None . 


Define 

Caution 

PATH:  o 


Specify  desired  printing  option  for  each 
history  hour.  When  NPRTGRID  is  set  to  1,  it 
overrides  NPRT  (GF-43)  and  causes  printing. 

When  NPRTGRID  is  set  to  0,  printing  is 
controlled  by  NPRT. 

END  GENERAL  FILE  BUILDING  SESSION 
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46.  When  the  user  Is  ready  to  initiate  the  run  programs  function,  all 

of  the  necessary  input  data  must  be  composed  and  formatted  into  the  proper 
files.  After  selecting  the  run  programs  fr  the  user  is  asked  to 

specify  a  reference  name  under  which  the  oul^c  uata  will  be  stored.  The 
reference  name  must  conform  to  the  file  name  conventions  of  the  UNIX  operating 
system.  Up  to  10  characters  can  be  specified  in  the  name.  The  name  must  NOT 
have  a  file  name  extension  (e.g.  ".dat"  or  ".out")  because  this  part  of  the 
file  name  is  specified  by  the  program.  The  user-specified  reference  name  is 
attached  to  three  output  files,  each  with  its  own  extension. 

47.  One  output  file  created  by  HURWIN  is  the  Final  Output  Wind  Field 
File  (given  the  ".win"  extension)  which  contains  a  date  header  and  the  wind 
speed  or  stress  matrix  and  direction  matrix  for  each  point  on  the  user- 
specified  output  grid  and  each  hourly  output  interval  of  the  simulation. 
Another  output  file  is  the  Wind  Characteristics  File  (given  the  ".sum" 
extension)  which  contains  time  histories  of  wind  speed  and  direction  for  each 
user-specified  output  station  (questions  GF-52  through  GF-59).  Wind  speeds 
are  in  m/sec  or  knots  (as  selected  in  GF-38).  Wind  directions  are  given 
either  in  the  coordinate  system  used  for  input  to  the  SHALWV  wave  model 
(Figure  13-8)  or  following  a  meteorological  convention,  depending  on  the 
selection  in  GF-39.  The  third  (optional)  output  file  is  the  Final  Output 
Pressure  Field  File  (given  the  ".prs"  extension),  which  contains  a  date  header 
and  the  pressure  matrix  for  each  point  on  the  user- specified  output  grid  and 
each  hourly  output  interval  of  the  simulation.  Pressures  are  in  millibars. 
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NOTE:  When  grid  information  is  entered  in  matrix  form  (as 
from  a  file),  then  the  first  line  of  data 
corresponds  to  ROW  8  on  the  grid  above. 


Figure  13-8.  Orientation  convention  for  the  computational  grid 
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PART  V:  HURVIN  FILE  DESCRIPTIONS 


48.  The  organization  of  files  used  in  the  HURVIN  branch  of  the  Spectral 
Vave  Modeling  module  is  shown  in  Figures  13-5  through  13-7.  A  list  of  the 
files  (and  related  file  numbers)  is  provided  in  Table  13-1.  Reference  numbers 
in  the  table  correspond  to  file  numbers  given  in  the  figures.  In  general,  the 
creation,  organization,  and  usage  of  almost  all  files  in  the  CMS  Spectral  Vave 
Modeling  module  are  transparent  to  the  user.  The  module  creates  and  organizes 
the  files,  prompting  the  user  for  information  only  when  necessary.  For 
reference,  however,  a  brief  description  of  the  data  and  data  format  required 
by  each  file  is  provided  in  this  section.  Default  values,  available  in  HURVIN 
for  some  data  fields,  are  given  in  parentheses. 


Snapshot  Incut  File  (snapin.hur) 

49.  The  Snapshot  Input  File  is  one  of  the  files  built  by  the  Input 
File -Building  routine  (Part  IV) .  The  file  contains  data  needed  to  compute 
hurricane  snapshots.  RECORD  3  of  this  file  takes  two  different  forms, 
depending  on  the  selection  in  GF-19,  because  of  the  additional  parameters 
required  when  the  general  pressure  profile  is  used. 

a.  RECORD  1,  Namelist  Format  (NAME1) 

Field  1  -  Desired  printing  option  (0) . 

Field  2  -  Number  of  snapshots. 

Field  3  -  System  of  units  (OCEANVEA) . 

b.  PECORD  2,  Namelist  Format  (NAME2) 

Field  1  -  Air-land  temperature  difference  (0.0). 

Field  2  -  Air-sea  temperature  difference  (-2.0). 


Field  2 
Field  3 


Field  4 
Field  5 
Field  6 


-  Boundary  layer  height  over  water  (500.0  m; 
1,640.4  ft). 

-  Roughness  length  over  land  (0.08). 

-  Chamock's  constant  (0.035). 

-  Potential  temperature  (26.85  deg  C; 

80.33  deg  F;  300.0  deg  K) . 


Field  7  -  von  Karman's  constant  (0.35). 
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Table  13-1 


Snapshot  Input  File 

snapin.hur 

12 

Hourly  History  Input  File 

histin.hur 

16 

Output  Grid  Parameters  Input  File 

gridparm.hur 

18 

Land/Sea  Table  Input  File 

lndsea.hur 

14 

Version  of  "snap”  Files 

snapl.hur 

snap2.hur 

21 

22 

Parameter  File 

filel.hur 

9 

Main  File  of  File  Names 

filenmi.hur 

8 

Post -Processing  File  Names 

filenmpst.hur 

15 

Job  Control  File  To  Run  Wind 
Computational  Programs 

hurwin_run 

10 

Snapshot  Wind  Information  File 

fort. 13 

13 

General  Run  Information  from  "snap" 

fort. 25 

25 

General  Run  Information  from  "hist" 

fort. 26 

26 

Initial  Wind  Field  or  Pressure 
on  Fine  Grid 

windOl  or  presOl 

7 

m 

Initial  Wind  Field  or  Pressure 
on  Coarse  Grid 

wind02  or  pres02 

V 

Wind  Field  Output  File 

user- specified 
+  ".win" 

ii 

Wind  Characteristics  Output  File 

user* specified 
+  " . sum" 

9 

Pressure  Field  Output  File 

user- specified 
+  ".prs" 

11 

c.  RECORD  3  (standard  pressure  profile,  GF-19,  option  "1"),  Namelist 
format  (NAME3)  (repeated  for  each  snapshot) 

Field  1  -  Storm  name. 

Field  2  -  Surface  geostrophic  wind  speed. 

Field  3  -  Surface  geostrophic  wind  direction. 

Field  4  -  North  latitude  of  storm  haracteristic  or 

average  value) . 

Field  5  -  Track  direction. 

Field  6  -  Forward  speed  of  storm. 
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Field  7  -  Pressure  at  storm  eye. 

Field  8  -  Pressure  profile  scale  radius. 

Field  9  •  Far  field  pressure. 

Field  10  -  Number  of  wind  computation  cycles  (800) . 

Field  11  -  Innermost  grid  spacing  (2.0  km;  1.2  miles; 

1.0  nm) .  Note  this  is  equal  to  the  innermost 
grid  spacing  used  for  computation  (GF-21)  only 
if  NESTS  -  7  (GF-20)  (see  Table  1,  Cardone  et 
al.  1994). 

Field  12  -  Index  to  the  number  of  active  nests  for 

computation  (parameter  INSIDE  in  Table  1, 
Cardone  et  al.  1994). 

Field  13  -  Distance  from  axis  to  0.5  surface  geostrophic 

wind  (0.0). 


Field  14  -  Units  for  track  direction  (0) . 

Field  15  -  Pressure  field  quadrant  indicator  (0) . 

Field  16  -  Switch  for  creating  pressure  field  output  (0) . 


RECORD  3  (generalized  pressure  profile,  GF-19,  option  "2"), 
Namelist  format  (NAME3A)  (repeated  for  each  snapshot) 


Field  1 
Field  2 
Field  3 
Field  4 

Field  5 
Field  6 
Field  7 
Field  8 
Field  9 
Field  10 
Field  11 
Field  12 


Field  13 

Field  14 


Storm  name. 

Surface  geostrophic  wind  speed. 

Surface  geostrophic  wind  direction. 

North  latitude  of  storm  eye  (characteristic  or 
aveugt;  value) . 

Track  direction. 

Forward  speed  of  storm. 

Pressure  at  storm  eye. 

Pressure  profile  scale  radius,  inner  component 
Pressure  profile  scale  radius,  outer  component 
Far  field  pressure. 

Number  of  wind  computation  cycles  (3200). 

Innermost  grid  spacing  (2.0  km;  1.2  miles; 

1.0  nm) .  Note  this  is  equal  to  the  innermost 
grid  spacing  used  for  computation  (GF-21)  only 
if  NESTS  -  7  (GF-20)  (see  Table  1,  Cardone  et 
al.  1994). 

Index  to  the  number  of  active  nests  for 
computation  (parameter  INSIDE  in  Table  1, 
Cardone  et  al.  1994). 

B  coefficient  for  inner  pressure  component 

(1.0). 
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Field  15 
Field  16 
Field  17 
Field  18 


-  B  coefficient  for  outer  pressure  component 

(0.0). 

-  Pressure  anomaly  for  inner  pressure  component 
(PFAR  -  EYPRES) . 

-  Pressure  anomaly  for  outer  pressure  component 

(0.0). 

-  Switch  for  creating  pressure  field  output  (0) . 


Hourly  History  Input  File  (histin.hur) 

50.  The  Hourly  History  Input  File  is  one  of  the  files  built  by  the 
Input  File-Building  routine  (Part  IV).  The  file  contains  data  needed  to 
compute  the  hourly  wind  histories. 

a.  RECORD  1,  Namelist  Format  (NAME4) 


Field  1 
Field  2 
Field  3 
Field  4 
Field  5 
Field  6 
Field  7 
Field  8 
Field  9 


-  Storm  name. 

-  Starting  time. 

-  Time  zone  for  starting  time  (1  -  UTC) . 

-  Desired  output  grid  (0  -  no;  1  -  yes). 

-  Interval  (hours)  to  print  winds  on  grid. 

-  Output  field  option  (1  -  wind  fields). 

-  Desired  wind  elevation  (10.0  m;  33.0  ft). 

-  Desired  output  wind  speed  units  ("knots"). 

-  Desired  output  wind  direction  convention 
("WIS"). 


b.  RECORD  2,  Namelist  Format  (NAMES) 

Field  1  -  Roughness  law  coefficients  (non-open 

water  terrains,  GF-14  for  defaults). 

Field  2  -  Number  of  measurement  stations. 

Field  3  -  Number  of  hours  for  history. 


c. 


RECORD  3,  Format  (514 , F6 . 1 , IX, 13)  Station  Location  and  Data 


(repeated  for 

each  measurement  station) 

Field  1 

-  Degrees  of  north  latitude. 

Field  2 

•  Minutes  of  latitude . 

Field  3 

-  Degrees  of  west  longitude. 

Field  4 

-  Minutes  of  longitude. 

Field  5 

-  Terrain  code  (see  GF-57  for  codes). 

Field  6 

-  Station  height  (10.0  m) . 

Field  7 

-  Station  number  (for  identification  only) 
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d. 


RECORD  4,  Format  (614 , F8 .4 , 214)  Hourly  History  Data  (repeated  for 
each  hour) 


Field  1 
Field  2 
Field  3 
Field  4 
Field  5 

Field  6 

Field  7 

Field  8 
Field  9 


-  Degrees  of  north  latitude,  storm  eye  location. 

-  Minutes  of  latitude,  storm  eye  location. 

-  Degrees  of  west  longitude,  storm  eye  location. 

-  Minutes  of  longitude,  storm  eye  location. 

-  Sequence  number  of  first  snapshot  to  be  used  for 
this  hour. 

-  Sequence  number  of  second  snapshot  wind  field  to 
be  used  for  this  hour  (blank  if  no  interpolation 
this  hour) . 

•  Interpolation  distance  between  first  and  second 
snapshots  (blank  if  Field  6  is  blank) . 

-  Clockwise  rotation  of  snapshot  (deg) . 

-  Printing  flag  (-  0  to  suppress  printing  of  the 
output  grid  wind  field) . 


Output  Grid  Parameters  Input  File  (gridparm.hur) 

51.  The  Output  Grid  Parameters  Input  File  is  one  of  the  files  built  by 
the  Input  File-Building  routine  (Part  IV).  It  is  optional  and  is  needed  only 
when  an  output  grid  is  requested.  The  file  contains  data  needed  to  define  the 
output  grid. 


a.  RECORD  1,  Format  (4F6.1.F6.2.2I4) 


Field  1 
Field  2 
Field  3 
Field  4 
Field  5 
Field  6 
Field  7 


-  North  latitude  for  north  edge  of  grid  (deg) . 

-  North  latitude  for  south  edge  of  grid  (deg) . 

-  West  longitude  for  west  edge  of  grid  (deg) . 

-  West  longitude  for  east  edge  of  grid  (deg) . 

-  Grid  cell  spacing  (deg) . 

-  Number  of  rows  (y-values  or  latitudes) . 

-  Number  of  columns  (x-values  or  longitudes). 


Land/Sea  Table  Input  File  (lndsea.hur) 

52.  The  Land/Sea  Table  Input  File  is  one  of  the  files  built  by  the 
Input  File-Building  routine  (Part  IV).  The  file  contains  data  needed  to 
define  the  type  of  land  or  water  surface  (terrain)  at  each  point  on  the  output 
grid.  The  file  is  needed  only  if  an  output  grid  is  requested. 


a.  RECORD  1,  Format  (7811)  (repeated  for  each  row  of  the  output  grid, 
beginning  with  the  most  northerly) 
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Field  1 


-  Terrain  code  for  first  (most  westerly)  grid 
point  In  row  (see  GF-S7  for  possible  codes). 

Field  2  -  Terrain  code  for  second  grid  point  In  row. 

Field  N  -  Terrain  code  for  Nth  grid  point  in  row. 

Vers ion- of -"Snap"  Files  (snapl.hur  and  snap2.hur) 

53.  The  Version-of-”Snap"  Files  are  very  short  files  created  and  deleted 
automatically  each  time  HURWIN  is  run.  Their  only  purpose  is  to  indicate  to 
program  "snap"  whether  a  coarse  or  fine  output  grid  is  to  be  run. 

a.  RECORD  1,  Format  (AS) 

Field  1  -  Version  of  "snap"  to  be  run  (-  snapl  for  coarse 

grid;  -  snap2  for  fine  grid). 

Parameter  File  (filel.hur) 

54.  The  Parameter  File  is  created  by  the  Input  File-Building  routine. 
Values  in  the  parameter  file  are  used  to  dimension  variables  in  HURWIN  and  its 
supporting  programs.  The  programs  use  the  FORTRAN  command  "INCLUDE"  to 
include  the  Parameter  File  in  the  programs  during  compilation.  The  Parameter 
File  contains  a  FORTRAN  "PARAMETER"  statement.  The  variables  in  the 
"PARAMETER"  statement  are  defined  as  follows: 

IDMN  -  Number  of  columns  (I  values)  in  the  output  grid. 

JDMN  -  Number  of  rows  (J  values)  in  the  output  grid. 

IF  -  1 

IA  -  1 

KSST  -  0 

MKST  -  0 

NBPS  -  0 

NDIF  -  0  to  suppress  pressure  field  output;  -  1  to  activate. 

ITDM  -  Number  of  hours  for  history. 

NDMN  -  Number  of  measurement  (gage)  stations. 

MODEL  -  3,  specifies  HURWIN  model. 

INDM  -  1 

MAXI  -  Number  of  columns  (I  values)  in  the  output  grid. 

MAXJ  -  Number  of  rows  (J  values)  in  the  output  grid. 

MAXHRS  -  Number  of  hours  for  time  history. 

MAXSTA  -  Number  of  measurement  (gage)  stations. 
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Main  File  of  File  Names  (f ilenmi .hur) 

55.  The  Main  File  of  File  Nanes  contains  the  name  of  the  main  input 
files.  It  also  contains  the  file  name  of  the  "snap”  FORTRAN  program  to  be 
used.  The  format  is  as  follows: 

a.  RECORD  1,  Free  Format. 

Field  1  -  Snapshot  input  file  name. 

b.  RECORD  2,  Free  Format. 

Field  1  -  Hourly  history  input  file  name. 

c.  RECORD  3,  Free  Format. 

Field  1  -  Output  grid  parameters  file  name. 

d.  RECORD  4,  Free  Format. 

Field  1  •  Land/sea  (terrain)  table  file  name. 

e.  RECORD  5,  Free  Format. 

Field  1  -  File  name  of  the  "snap”  program  to  be  used 

(-  'snap_adc.f'  for  standard  pressure  profile; 

-  'snap_hol.f'  for  generalized  pressure  profile). 


Post-Processing  File  Name  File  (filenmpst.hur) 

56.  The  Post-Processing  File  Name  File  contains  the  names  of  the  output 
files  from  HURWIN.  The  format  is  as  follows: 

a.  RECORD  1,  Free  Format. 

Field  1  -  File  name  for  wind  characteristics  output  file 

(used  for  creating  time  history  plots  and 
statistics) . 

b.  RECORD  2,  Free  Format. 

Field  1  -  Blank. 

c.  RECORD  3,  Free  Format. 

Field  1  -  File  name  for  wind  field  output  file. 


Job  Control  File  To  Run  Wind  Computational  Programs  (hurwinrun) 

57.  The  file  of  job  control  language  (SCRIPT  file)  used  to  run  the  wind 
computational  programs  is  created  automatically  when  HURWIN  is  run.  It 
controls  file  handling  and  execution  of  the  FORTRAN  programs  needed  to  run 
HURWIN.  Because  of  its  specialized  nature,  it  is  not  discussed  in  detail 
here. 


Snapshot  Wind  Information  File  (fort. 13) 

58.  The  Snapshot  Wind  Information  File  contains  detailed  information  on 
each  snapshot  computed  in  program  "snap."  For  each  snapshot,  the  file 
includes  parameters  and  the  full  wind  field  on  the  nested  computational  grid. 
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This  file  is  created  and  deleted  automatically  when  HURWIN  is  run.  Because  it 
is  written  in  binary  form,  it  is  not  discussed  in  detail  in  this  report. 
Initial  Wind  or  Pressure  Field  on  Fine  Grid  File  (windOl  or  presOl) 

59 .  The  Initial  Vind  or  Pressure  Field  on  Fine  Grid  File  contains  wind 
or  pressure  fields  for  each  requested  hour.  They  are  computed  using  the  user- 
specified  grid  spacing  for  the  nested  computational  grid  (specified  in  GF-21 
for  the  innermost  computational  nest) ;  but  they  are  output  in  terms  of  the 
user-specified  output  grid.  This  file  is  created  and  deleted  automatically 
when  HURUIN  is  run  and  an  output  grid  has  been  requested.  Pressure  fields  are 
optional  (GF-40) . 

a.  RECORD  1,  Format  (110)  (Records  1-3  are  repeated  for  every  hour 
requested  in  the  hourly  history) 

Field  1  -  Date -time  for  the  wind  field  (YYMMDDHH) . 

b.  RECORD  2,  Format  (21F6.1)  for  wind  speed;  (21F6.3)  for  wind 
stress;  (21F6.1)  for  pressure  (repeated  for  successive  rows  in 
maximum  groupings  of  21  columns) 

Field  1  -  Wind  speed,  stress  (u.2;  see  GF-32),  or  pressure 

at  northwest  comer  of  output  grid  (knots)  . 

Field  2  -  Wind  speed,  stress,  or  pressure  at  first  point 

east  of  northwest  corner  of  output  grid. 


c.  RECORD  3,  Format  (21F6.1)  (repeated  as  with  RECORD  2  for  wind 
speed  or  stress;  omitted  for  pressure) 

Field  1  -  Wind  direction  at  northwest  corner  of  output 

grid  (deg;  0  deg  -  wind  blowing  toward  east; 
90  deg  -  wind  blowing  toward  north;  etc.;  see 
Figure  13-8). 


(same  repetition  as  with  RECORD  2) 

Initial  Wind  or  Pressure  Field  on  Coarse  Grid  File  (wind02  or  pres02) 

60.  The  Initial  Wind  or  Pressure  Field  on  Coarse  Grid  File  contains 
wind  or  pressure  fields  for  each  requested  hour.  They  are  computed  using  an 
initial  coarse  spacing  for  the  nested  computational  grid,  which  is  preset  in 
HURWIN  (25-km  spacing  for  the  inner  nest);  but  they  are  output  in  terms  of  the 
user-specified  output  grid.  This  file  is  created  and  deleted  automatically 
when  HURWIN  is  run  and  an  output  grid  has  been  requested.  The  record 
structure  and  formats  are  identical  to  those  for  the  Initial  Wind  or  Pressure 
Field  on  Fine  Grid.  They  are  not  repeated  here. 
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Wind  Field  Output  File  (user  specified  +  ".win") 

61.  The  Wind  Field  Output  File  is  the  file  of  final  hourly  hurricane 
wind  fields  on  the  user-specified  output  grid.  It  is  created  in  progran 
"blend"  by  aerging  the  initial  wind  fields  computed  on  the  coarse  and  fine 
grids  onto  the  output  grid.  Winds  froa  the  fine  grid  are  used  in  all  areas  of 
the  output  grid  covered  by  the  fine  grid.  Winds  froa  the  coarse  grid  are  used 
for  areas  of  the  output  grid  (if  any)  which  lie  outside  the  fine  grid 
coverage.  Record  structure  and  formats  are  identical  to  those  for  the  initial 
wind  fields  on  fine  and  coarse  grids.  They  are  not  repeated  here. 

Wind  Characteristics  Output  File  (user  specified  +  ".sum") 

62.  The  Wind  Characteristics  Output  File  contains  information  on  the 
hourly  time  history  of  winds  from  selected  special  output  (measurement) 
stations.  The  number  of  stations  and  descriptive  data  for  each  selected 
station  are  described  in  questions  GF-52  through  GF-59.  Records  are  organized 
as  follows: 


a. 

RECORD  1 

-  Station  1, 

Hour 

1 

b. 

RECORD  2 

-  Station  1, 

Hour 

2 

c. 

RECORD  3 

-  Station  1, 

Hour 

3 

n. 

RECORD  n 

-  Station  1, 

Hour 

NOHRS 

n+1. 

RECORD  n+1- 

Station  2 ,  Hour  1 

n+2. 

RECORD  n+2- 

Station  2,  Hour  2 

Each  record  contains  the  following  data: 

RECORD  n,  Format  (110 , 215 , 2F6 . 1) 

Field  1  -  Date  (YYMMDDHH) . 

Field  2  -  Hour. 

Field  3  -  Sequential  station  number. 

Field  4  -  Wind  speed  (knots) . 

Field  5  -  Wind  direction  (degrees,  see  Figure  13-8). 
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Pressure  Field  Output  File  (user  specified  +  ".prs") 

63.  The  Pressure  Field  Output  File  is  the  file  of  final  hourly 
hurricane  pressure  fields  on  the  user-specified  output  grid.  It  is  created  in 
program  "blendp"  by  merging  the  initial  pressure  fields  computed  on  the  coarse 
and  fine  grids  onto  the  output  grid.  Pressures  from  the  fine  grid  are  used  in 
all  areas  of  the  output  grid  covered  by  the  fine  grid.  Pressures  from  the 
coarse  grid  are  used  for  areas  of  the  output  grid  (if  any)  which  lie  outside 
the  fine  grid  coverage.  Record  structure  and  formats  are  identical  to  those 
for  the  initial  pressure  fields  on  fine  and  coarse  grids.  They  are  not 
repeated  here. 
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PART  VI:  EXAMPLE  APPLICATION 


64.  An  example  application  is  presented  in  this  section  to  illustrate 
the  use  of  the  HURWIN  model,  including  input  files  and  model  output.  The 
example  also  serves  as  a  benchmark  against  which  future  versions  of  the  model 
can  be  checked. 

65.  The  example  is  based  on  the  passage  of  a  typhoon  near  the  island  of 
Guam  in  the  south  Pacific  Ocean  (Northern  Hemisphere) .  Typhoon  Russ  impacted 
Guam  during  December  of  1990  The  basic  snapshot  information  for  Typhoon  Russ 
was  developed  by  WIS  from  the  1990  Joint  Typhoon  Warning  Center  publication, 
1990  Annual  Tropical  Cyclone  Report  (NOCC/JTWC  1990) . 

66.  The  fine  mesh  computational  grid  was  run  with  a  grid  spacing  of 
10  km  on  the  innermost  nest.  An  output  grid  was  defined  with  13  rows, 

13  columns,  and  a  grid  spacing  of  0.1  deg.  The  grid  boundaries  are  12.9-14.1 
deg  north  latitude  and  214.9-216.1  west  longitude  (143.9-145.1  deg  east 
longitude).  The  grid  encompassed  the  island  of  Guam.  To  represent  the  island 
and  surrounding  reef  areas,  three  points  on  the  grid  were  defined  as  "woods" 
(or  land)  and  five  points  were  defined  as  "lake."  While  the  "lake"  term  may 
seem  unusual  for  this  application,  the  corresponding  surface  roughness  values 
give  a  rough  approximation  of  the  actual  situation.  The  remaining  grid  points 
were  defined  as  "open  water."  Input  files  for  simulating  2  hr  (not 
consecutive)  during  the  passage  of  Typhoon  Russ  are  given  in  Figure  13-9. 

67.  During  the  first  hour  simulated  (1500  hr  on  20  December  1990) 
Typhoon  Russ  was  nearly  due  south  of  Guam.  Maximum  wind  speed  on  the  output 
grid  exceeded  95  mph.  Wind  speed  contours  are  shown  in  Figure  13-10.  Because 
the  storm  center  is  relatively  near  (but  outside)  the  grid  boundary,  wind 
contours  have  a  noticeable  curvature.  Wind  speeds  were  reduced  significantly 
at  Guam  because  of  the  increased  roughness  specified  for  those  grid  points. 
Wind  speed  and  direction  across  the  grid,  as  given  in  the  output  file 
"guam.win,"  are  tabulated  in  Figure  13-ii. 

68.  The  second  hour  simulated  is  0000  hr  on  21  December  1990.  (In  a 
normal  HURWIN  run,  wind  rieius  would  be  simulated  every  hour).  The  storm 
center  had  moved  to  a  position  due  west  of  Guam.  The  center  is  more  than 

1  deg  west  of  the  western  boundary  of  the  output  grid.  Maximum  wind  speed  on 
the  grid  is  58.6  mph.  Wind  speed  contours  are  reasonably  straight  across  the 
grid,  with  the  exception  of  distortions  around  Guam  (Figure  13-12).  Tabulated 
wind  speeds  and  directions  are  given  in  Figure  13-13. 
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SNAPIN.HUR  tout  FUe 

ANAME1  IB-0,  NZ-3,  UNITS  -  ’  METRIC’,  AEND 


ANAME2  DTH-O.O,  -2.0,  HH= 500.0,  ZOLAND-8.0E-O2, GARR- 1 .44E-02,  PTH -26.85,  K35-0.35,  AEND 


ANAME3STORM-’nu»’,  SGW-9.0,  ANl- 161.0,  EYELAT- 12.0,  EYLONG-215.0,  DIREC-289.0,  SPEED-6.17333, 
EYPRES-924.0,  RADIUS-37.04, 3*0.0,  PFAR- 1012.0,  NM— 800,  DX-10.0,  ST12-0.0,  TTRACK-0, 1QUAD-0,  AEND 

ANAME3  STORM-  ’roi»\  SGW-9.0,  ANl  - 161 .0,  EYELAT- 12.0,  EYLONG-215.0,  DIREC-289.0,  SPEED -7.202222, 
EYPRES —924.0,  RADIUS-37.04, 3*0.0,  PFAR-1012.0,  NM-800,  DX-10.0,  ST12-0.0,  TTRACK-0, 1QUAD-0,  A  END 

&NAME3  STORM  — ’niM’,  SGW-9.0,  ANl  -161.0,  EYELAT- 12.0,  EYLONG-215.0,  DIREC-289.0,  SPEED-8.231111, 
EYPRES— 924.0,  RADIUS-37.04, 3*0.0,  PFAR-1012.0,  NM-800,  DX-10.0,  ST12-0.0,  TTRACK-0,  IQUAD-0,  A  END 


H1STIN  .HUR  Incut  FUe 

ANAME4  STORM  =  ISTART-90122015,  ZONE-  ’UTM’,  ICNVRT-1,  NPRT-1,  IFIELD-1,  ORIDHT- 10.0,  A  END 


AN  AMES  ZCOEFF— 0.0, 3.7E-03,  1.5E-02, 2*0.0,  4.0E-02,  2*0.0, 0.16,  2*0.0,  0.32,  2*0.0,  1.28,  NOGAGES— 0, 
NOHRS-2,  AEND 

12  30  215  9  1  2  0.3333  0  0 
12  58  217  22  2  2  0.8333  0  0 


GRIDPARM.HUR  Input  File 
14.1  12.9  216.1  214.9  0.10  13  13 


LNDSEA-HUR  Iimut  File 

1111111111111 

1111111111111 

1111111111111 

1111111111111 

1111111111111 

1111111111111 

1111111125111 

1111111251111 

1111111251111 

1111111221111 

1111111111111 

1111111111111 

1111111111111 


Figure  13-9. 


Exanple  data  files,  Typhoon  Russ 
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256.6  254.8  253.0  250.8  248.7  246.3  243.8  241.1  238.4  235.6  232.6  229.6  226.6 

258.5  256.7  254.8  252.5  250.3  247.7  245.0  242.1  239.2  236.1  232.9  229.6  226.3 

260.6  258.7  256.7  254.3  252.0  249.2  246.3  243.2  240.0  236.6  233.1  229.6  225.9 

262.9  260.9  258.9  256.4  253.9  250.8  247.8  244.4  240.9  237.1  233.3  229.3  225.3 

265.3  263.2  261.2  258.5  255.9  252.6  249.3  245.5  241.7  237.6  233.4  229.0  224.6 

268.0  265.9  263.8  261.0  258.3  254.7  251.1  246.9  246.1  247.6  233.2  228.3  223.3 

270.7  268.6  266.5  263.5  260.6  256.8  252.9  251.6  253.1  238.4  233.1  227.5  222.0 
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280.8  278.9  276.5  273.7  270.1  265.7  260.3  253.9  246.4  238.0  229.0  219.5  210.6 
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288.8  287.4  285.3  282.7  279.2  274.5  268.2  259.5  248.0  233.8  217.6  201.7  188.8 


Figure  13-11.  Output  wind  speeds  and  directions, 
1500  hr,  20  Dec  90 
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Wind  Speed  Contours 
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Figure  13-12.  Wind  speed  contours,  0000  hr,  21  Dec  90 
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Figure  13-13.  Output  wind  speeds  and  directions, 
0000  hr,  21  Dec  90 
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CHAPTER  14 

WICM:  A  WAVE -INDUCED  CURRENT  MODEL 
THEORY  AND  PROGRAM  DOCUMENTATION 

PART  I:  INTRODUCTION 

1.  This  chapter  documents  the  Wave -induced  Current  Jjodel  (WICM) .  WICM 
is  a  two-dimensional,  depth-averaged  model  for  computing  wave- induced  currents 
and  water  surface  setup.  The  model  solves  finite  difference  approximations  of 
the  Navier-Stokes  (continuity  and  horizontal  momentum)  equations  for  water 
surface  displacement  (S)  and  unit  flow  rate  components  ( U  and  V) .  WICM  is  a 
two-dimensional  model;  therefore,  velocities  are  treated  as  depth- averaged 
quantities  (i.e.,  velocities  are  constant  through  the  water  column  in  both 
magnitude  and  direction).  WICM  can  simulate  flow  fields  induced  by  waves, 
winds,  river  inf lows/outf lows ,  and  tidal  forcing.  This  finite  difference 
model  is  developed  in  boundary -fitted  (curvilinear)  coordinates. 

2.  WICM  is  a  specialized  version  of  the  curvilinear  long  wave  hydrody¬ 
namic  model  CLHYD  documented  in  Chapter  6.  In  addition  to  the  forcing 
mechanisms  in  CLHYD,  WICM  includes  wave  stresses  (gradients  in  radiation 
stress)  caused  by  breaking  waves  to  force  nearshore  circulation  and  wave 
setup.  To  calculate  radiation  stress,  WICM  requires  input  wave  fields  from 
either  a  monochromatic  (Regional  Coastal  Processes  Wave  Propagation  model 
(RCPWAVE) ,  Chapter  5)  or  spectral  (Steady-State  Spectral  Wave  Energy  Propaga¬ 
tion  Model  (STWAVE),  Chapter  8)  wave  model. 

3.  The  user  must  have  a  thorough  understanding  of  the  model's  capabili¬ 
ties  and  limitations  before  applying  it  to  a  particular  study.  WICM  does  not 
provide  the  "total  solution"  to  any  given  hydrodynamic  problem.  The  user  must 
ensure  that  limitations  imposed  by  the  governing  equations  (e.g.,  linear  wave 
theory)  are  applicable  to  the  problem  being  investigated.  Furthermore,  the 
model  should  not  be  treated  as  a  "black  box" ;  the  engineer  or  scientist  must 
check  computed  results  and  assess  their  reasonableness. 

4.  This  chapter  is  divided  into  five  sections:  Part  II  presents  the 
governing  equations  and  computational  scheme  used  in  the  model,  Part  III 
defines  the  input  data  formats,  Part  IV  discusses  the  model's  input  data 
requirements,  and  Part  V  contains  several  illustrative  examples. 
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PART  II:  WICM  MODEL  FORMULATION 


5.  The  concept  of  radiation  stress,  introduced  by  Longuet-Higgins  and 
Stewart  (1962,  1963,  1964),  has  been  used  to  model  wave-induced  currents  and 
setup  in  one  dimension  (Bowen  1969;  Longuet-Higgins  1970a,  b;  and  Thornton 
1970)  and  two  dimensions  (Keely  and  Bowen  1977;  Ebersole  and  Dalrymple  1979; 
Vemulakonda  1984;  and  Wind  and  Vreugdenhil  1986).  Radiation  stress  is  a 
measure  of  the  momentum  flux  in  the  wave  field.  Gradients  in  radiation  stress 
due  to  decay  of  the  wave  height  (wave  breaking)  exert  a  stress  on  the  water 
column.  This  wave  stress  transfers  momentum  from  the  wave  field  to  the 
longshore  and  cross-shore  flow  field.  This  chapter  documents  the  two- 
dimensional  WICM  model  for  simulating  wave- induced  currents  and  water  surface 
setup.  In  two-dimensional  models,  the  governing  three-dimensional  equations 
are  integrated  over  the  water  depth  to  yield  vertically  averaged  velocities; 
that  is,  velocities  that  are  constant  through  the  water  column. 

Assumptions  and  Limitations 

6.  Proper  application  of  the  model  requires  a  clear  understanding  of 
the  physical  processes  occurring  in  a  study  area  and  comprehension  of  the 
capabilities  of  the  model  to  simulate  those  processes.  The  limitations  of  a 
model  define  its  range  of  applicability.  In  particular,  WICM  is  a  two- 
dimensional  depth- averaged  model;  therefore,  the  model  should  be  applied  only 
where  the  water  column  is  well  mixed  and  no  significant  vertical  variations 
occur.  Hydrostatic  pressure  conditions  are  assumed  in  the  model  formulation. 
Thus,  the  model  should  be  applied  only  where  there  is  no  significant  vertical 
acceleration  of  the  water. 

7.  Responsible  application  of  WICM  requires  knowledge  of  the  model's 
capabilities  as  well  as  its  limitations.  The  model  should  be  applied  such 
that  time  and  length  scales  associated  with  breaking  short-wave  and  long-wave 
(for  simulations  including  tide)  processes  can  be  resolved.  WICM  is  most 
applicable  for  open  coast  situations  where  the  primary  focus  is  surf  zone 
circulation.  Assumptions  in  the  model  include  linear  wave  theory  (calculation 
of  radiation  stress) ,  no  wave-current  interaction,  linear  bottom  friction,  and 
no  significant  velocity  variations  over  depth. 


14-2 


8.  Presently,  the  model  lacks  a  flood/dry  capability  and  cannot  treat 
submerged  or  overtopping  barriers.  A  thorough  comprehension  of  the  physical 
processes  simulated  by  the  model  is  necessary  to  ensure  that  the  model  is 
applied  to  appropriate  problems,  that  it  is  applied  correctly,  and  that 
accurate  results  are  produced. 

9.  A  discussion  of  the  hydrodynamic  equations  used  in  UICM  is  provided 
in  the  following  section.  It  is  recommended  that  the  reader  refer  to  Horikawa 
(1988)  for  a  detailed  discussion  of  coastal  hydrodynamics. 


Governing-  Equations 


10.  The  hydrodynamic  equations  used  in  WICM  are  derived  from  the 
classical  Navier-Stokes  equations  formulated  in  a  Cartesian  coordinate  system 
(Figure  14*1).  If  vertical  water  accelerations  are  assumed  to  be  small 
compared  with  gravitational  acceleration  (hydrostatic  pressure  conditions 
exist)  and  the  fluid  is  homogeneous,  the  depth -averaged  approximation  is 
appropriate  and  yields  the  following  two-dimensional  form  of  the  governing 
equations : 


^  *  4(^1 '  i(^)  *  - 
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VI 


unit  flow  rate  components  in  the  x-  and  y-directions , 
respectively  (discharge  per  unit  width) 

independent  space  and  time  variables 

total  water  depth  (h+S) 

static  water  depth  measured  from  the  same  vertical  datum 

water  surface  displacement  measured  relative  to  an  arbitrary 
vertical  datum 
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Figure  14-1.  Definition  sketch  for 
Cartesian  coordinate  system 


g 

f 


P 


gravitational  acceleration 
Coriolis  parameter 

external  surface  shear  stresses,  such  as  wind  stress,  in  the 
x-  and  y-directions ,  respectively 

fluid  density  (assumed  to  be  constant) 

wave  stress  (gradients  in  radiation  stress),  in  the  x-  and  y- 
directions,  respectively 


rfc,  rt,  “  bottom  friction  components  in  the  x-  and 
y-directions,  respectively 

A h  -  generalized  diffusion  coefficient 
Equations  14-1,  14-2,  and  14-3  represent  the  x-momentum,  y-momentum,  and 
continuity  equations,  respectively. 

11.  A  detailed  discussion  of  the  Navier-Stokes  equations  with  a 
rigorous  derivation  of  each  term  is  found  in  Harris  and  Bodine  (1977) .  Wave 
stress  and  friction  terms  are  derived  in  Longuet-Higgins  (1970a,  b) .  A  brief 
discussion  of  the  physical  significance  of  the  six  groups  of  terms  in 
Equations  14-1  through  14-3  is  given  in  Table  14-1  and  below.  The  Roman 
numerals  below  correspond  to  those  in  Equations  14-1  through  14-3: 
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Table  14-1 


Description  of  Terns  in  the  Governing  Equations 


Term 

Definition  and  Discussion 

dU  dv 

It'  ~5i 

Change  of  the  vertically  averaged  flow  per  unit  width 
with  respect  to  time .  Change  may  result  from  temporal 
acceleration  of  the  mean  flow. 

dS  „  dS 

gH~5x'  gHSP 

Pressure  gradient  terms;  describes  the  slope  of  the 
water  surface;  principal  driving  force  of  fluid  flow. 

Bottom  friction  terms;  stress  of  the  fluid  layer 
against  the  bottom  boundary;  serves  as  an  energy  dissi- 
pator . 

fu,  fv 

Coriolis  terms;  accounts  for  the  effect  of  the  Earth's 
rotation. 

d  I  UU\ 

^hrj 

Advective  (inertia)  terms;  describes  the  movement  of 
water  due  to  the  fluid  motion  itself. 

T«r»  xj y 

External  shear  stresses  (such  as  the  wind) ;  any  forcing 
function  that  serves  to  drive  the  fluid  motion. 

Xuxl  xtiy 

Wave  stresses;  forcing  by  gradients  in  radiation  stress 
caused  by  wave  breaking. 

A(JH£  + 

idx *  ay*  J 

Horizontal  diffusion  of  momentum;  describes  the  diffu¬ 
sion  of  momentum  due  to  fluid  motion;  A„  is  sometimes 
referred  to  as  the  eddy  viscosity. 

as 

Change  in  water  level  with  respect  to  time. 

dU  dU 

Hx'  Hy 

Rate  at  which  water  is  converging  or  diverging  horizon¬ 
tally  at  a  given  point  (x,y)  in  space. 
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I.  Local  flow  acceleration  (i.e.,  local  variation  of  momentum  with 
respect  to  time) . 

II.  Transport  of  momentum  by  advection  (i.e.,  spatial  acceleration). 

III.  Barotropic  pressure  forces  and  conservation  of  mass. 

IV.  Momentum  sources  and  sinks  due  to  Coriolis  force,  surface  wind 
stress,  and  wave  stress. 

V.  Momentum  sink  due  to  bed  friction. 

VI.  Horizontal  diffusion  of  momentum. 

12.  Various  formulations  of  the  terms  in  the  governing  equations  are 
permissible.  The  expressions  for  bottom  friction,  diffusion,  wind  stress, 
wave  stress,  and  the  Coriolis  coefficient  employed  in  WICM  are  given  below: 
Bottom  friction 

13.  The  bottom  shear  stress  impeding  fluid  motion  is  expressed  in  terms 
of  a  linear  friction  law  (Longuet-Higgins  1970a) : 

-  P C,  <|uorb|>  U 
*By=  PC,  <!uorJ)|>  V 


(14-4) 

(14-5) 


where 

C,  -  constant  friction  factor  (order  0.01) 

-  absolute  value  of  the  wave  orbital  velocity  at  the  bottom, 
averaged  over  one  wave  period 

From  linear  wave  theory,  the  absolute  value  of  the  wave  orbital  velocity  at 
the  bottom,  averaged  over  one  wave  period,  is  given  by: 


<Krti> 


2  H 

T  sinh(ic  \Ha\) 


(14-6) 


where 

H,  -  wave  height 

T  -  wave  period 

k  -  wave  number 

Equations  14-4  through  14-6  are  based  on  the  assumption  that  the  velocity 

components  U  and  V  are  small  compared  with  the  wave  orbital  velocity. 
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Diffusion  coefficient 

14.  The  diffusion  coefficient  Ah  is  treated  as  a  constant  in  UICM  and 
describes  the  rate  of  diffusion  of  momentum  due  to  fluid  motion.  Ah  is  a 
function  of  the  grid  spacing  with  typical  values  of  0.0  to  1.0  mJ/sec.  Due  to 
the  current  representation  of  the  advective  terms  in  WICH,  results  are  not 
extremely  sensitive  to  A„  . 

Wind  stress  coefficient 

15.  Wind  stress  r,  is  formulated  as: 

w\w  (14-7) 

where  is  air  density,  CD  is  the  wind  drag  coefficient  determined  from 

Garratt's  equation  (Garratt  1977): 

r  =  (0-75  +  0 . 067m)  (U.fl) 

D  1000 


where  w  is  the  resultant  wind  speed  (meters/second)  at  10  m  above  the  water 
surface,  and  W  is  the  wind  velocity. 

Wave  stress 

16.  Wave  stress  r„  is  formulated  in  terms  of  gradients  in  radiation 
stress : 


dS„  _  dS^ 
dx  dy 


(14-9) 


T 


►y 


(14-10) 


where  S„,  S^,  and  are  components  of  radiation  stress,  defined  as: 

=  E  [n  (cos2  (a)  +1.0)  -  0.5]  (14-11) 


-  E  n  cos (a)  sin(a)  (14-12) 

Syy  =  E  [(2 n  -  0.5)sin2(a)  +  (n  -  0.5)cos2(a)]  (14-13) 

where  E  is  wave  energy,  n  is  the  ratio  of  group  to  phase  velocity,  and  m 
is  the  local  wave  direction.  Wave  energy  is  calculated  from  linear  wave 
theory  using  monochromatic  wave  input  (from  RCPWAVE)  as: 
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E  »  0 . 125  p  g  (H,) 2 


(14-14) 


where  H is  significant  wave  height,  and  using  spectral  wave  input  (from 
STUAVE)  as: 

E  -  0.0625  p  g  ( H „)*  (14-15) 

where  H ^  is  the  energy-based,  zero-moment  wave  height.  Further  details  on 
the  calculation  of  wave  stress  based  on  radiation  stress  are  given  in  Longue t- 
Higgins  (1970a). 
fiailftlli  coefficient; 

17.  The  Coriolis  tens  accounts  for  the  fact  that  the  Earth  is  rotating, 
whereas  the  coordinate  frame  of  the  computations  is  fixed.  The  Coriolis 
parameter  f  is  expressed  as: 

f  -  2vsini  (14-16) 

where  v  is  the  angular  speed  of  the  Earth's  rotation  (7.292  x  10'5  rad/sec) 
and  A  is  the  latitude  of  the  study  area.  The  extent  of  nearshore  current 
problems  is  usually  small  enough  that  the  Coriolis  term  may  be  neglected. 

18.  This  completes  the  basic  model  formulation.  If  further  details  are 
desired,  the  reader  should  refer  to  Roache  (1976)  or  Horikawa  (1988). 

Grid  Systems 

19.  The  governing  equations  (Equations  14-1,  14-2,  and  14-3)  that 
describe  the  physical  processes  associated  with  nearshore  circulation  contain 
partial  derivatives  with  respect  to  time  and  space.  Model  UICM  uses  mathe¬ 
matical  (finite  difference)  approximations  to  represent  these  continuous 
equations.  The  continuum  is,  therefore,  represented  by  discrete  points  in 
time  and  space.  Discretization  of  the  horizontal  plane  is  accomplished  via  a 
computational  grid  composed  of  a  lattice  network  of  cells.  Each  cell  has 
certain  flow  field  parameters  associated  with  it.  In  the  case  of  UICM,  the 
water  surface  elevation  is  defined  at  the  center  of  each  cell  and  the  unit 
flow  rate  components  are  defined  on  the  cell  faces.  All  information  required 
as  input  to  the  model,  including  water  depths  and  external  forces  (such  as 
wind),  is  also  defined  at  each  cell  center. 
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20.  WICM  is  capable  of  using  a  uniform  Cartesian  grid,  a  stretched 
Cartesian  grid,  or  a  general  curvilinear  grid  (Figures  14-2  through  14-4). 

Each  of  these  grid  types  provides  a  successively  more  accurate  representation 
of  the  modeled  area;  however,  the  complexity  of  the  governing  equations 
increases  with  each  successive  grid  type.  The  grid  type  for  WICM  must  match 
the  grid  type  used  by  the  wave  model  providing  the  input  wave  field  for  wave 
stress  calculations  (uniform  Cartesian  grid  (RCPWAVE  and  STWAVE) .  The 
stretched  Cartesian  and  general  curvilinear  grid  options  in  WICM  are  reserved 
for  possible  future  applications. 

21.  Uniform  Cartesian  grids  simply  have  cells  of  equal  size  in  the  x- 
direction  and  equal  size  in  the  y-direction.  The  advantages  of  using  this 
type  of  grid  are:  (a)  the  simplicity  of  generating  a  grid;  (b)  the  simplicity 
of  formulating  the  governing  equations;  (c)  the  computational  requirements  are 
minimized  because  there  are  fewer  terms  in  the  governing  equations;  and 

(d)  the  uniform  grids  have  been  widely  used  and  thoroughly  tested. 


Figure  14-2.  Uniform 
Cartesian  grid 


Figure  14-3. 
Cartesian 


Stretched 

grid 


Figure  14-4.  General  curvilinear  grid 
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Disadvantages  of  using  a  uniform  Cartesian  grid  include:  (a)  inability  to 
increase  grid  resolution  in  areas  of  interest;  (b)  usually  uneconomical  to  use 
a  fine -resolution  grid  throughout  ntire  area;  and  (c)  inability  to  effi¬ 
ciently  resolve  an  irregular  shor  .,«s  thus  requiring  "stair-stepping, “  which 
can  lead  to  computational  inaccuracies  (Weare  1979) . 

22.  Stretched  Cartesian  grids  are  flexible  enough  to  increase  grid 
resolution  and,  therefore,  provide  a  more  economical  represertation  of  an  area 
than  a  uniform  Cartesian  grid.  In  other  words,  the  number  of  grid  cells 
required  to  represent  an  area  can  be  reduced  by  using  a  stretched  grid,  which 
in  turn  reduces  computational  costs.  However,  the  user  must  also  consider 
that  the  number  of  terms  in  the  governing  equati  - .s  increases  when  going  from 
a  uniform  to  a  stretched  grid. 

23.  Neither  the  uniform  nor  the  stretched  Ca>  ejian  grids  can  accurate¬ 
ly  depict  irregular  shorelines  and  must  therefore  resort  to  "stair-stepping." 
However,  nonorthogonal  curvilinear  (boundary- fitted)  grids  can  be  generated  to 
conform  to  bathymetric  features  and  provide  a  more  accurate  means  of  repre¬ 
senting  a  specific  study  area  (Figure  14-4).  These  grids  can  be  generated 
using  a  numerical  grid  generator  such  as  program  EAGLE  (Thompson,  Uarsi,  and 
Mastin  1985),  which  has  the  flexibility  to  concentrate  grid  lines  in  shal¬ 
low/deep  areas  or  in  areas  where  the  bathymetric  gradients  are  great.  With 
the  increase  in  accuracy  and  adaptability  comes  increased  complexity  of  the 
governing  equations  and  the  associated  computational  costs.  In  addition, 
general  curvilinear  grids  are  complex  networks  and  are  therefore  more  diffi¬ 
cult  to  generate.  However,  increased  accuracy  and  adaptability  are  often 
necessary  to  adequately  represent  features  in  coastal  areas. 

General  curvilinear  grid 

24.  For  areas  having  irregular  geometries,  WICM  is  capable  of  using  a 
nonorthogonal  curvilinear  grid.  Thompson  (1983)  developed  a  method  for 
generating  two-dimensional  boundary- fitted  grids  by  solving  elliptic  equa¬ 
tions  .  These  equations  relate  the  nonorthogonal  curvilinear  coordinate  system 
in  the  physical  plane  (x,y)  with  a  uniformly  spaced  coordinate  system  in  the 
transformed  plane  (£,»;)  (Figure  14-5).  The  elliptic  equations  are: 
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Figure  14-5.  Physical  and  transformed  planes 


(14-17) 

♦!»„-(?  (14-18) 

with  the  following  boundary  conditions: 

q  «  ^constant  }  OD  nozth  and  south  boundaries 

jj  “  ^°xfyfnt  }  on  ®astr  ^  west  boundaries 

where  the  functions  P  and  Q  may  be  chosen  to  obtain  the  desired  grid 
resolution  and  alignment. 

25.  Vhen  using  a  nonorthogonal  curvilinear  grid,  the  governing  equa¬ 
tions  must  be  transformed  into  curvilinear  coordinates.  A  straightforward 
method  is  to  transform  only  the  independent  variables  (x,y)  and  then  solve  for 
the  Cartesian  components  of  the  velocity  and  water  surface  displacement 
(Johnson  1982).  The  advantage  of  this  method  is  its  simplicity  in  generating 
the  transformed  equations  via  the  chain  rule;  although  the  resulting  equations 
are  quite  complex.  Disadvantages  are:  (a)  the  boundary  conditions  are 
complicated  because  the  Cartesian  velocity  components  are  generally  not 
aligned  with  the  grid  lines;  (b)  a  staggered  grid  (i.e.,  a  grid  in  which 
values  of  velocity  and  water  surface  elevation  are  defined  at  different 
points)  cannot  be  readily  used;  and  (c)  numerical  instabilities  may  develop 
(Sheng  1986). 
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26.  To  avoid  these  problems,  Sheng  transformed  the  dependent  (S ,U ,V)  as 
well  as  the  independent  variables  (x,y).  Vector  quantities  (e.g.,  U ,V)  are 
transformed  by  multiplying  the  vector  by  a  scale  factor.  Scalar  quantities 
(e.g.,  S )  in  the  physical  plane  are  the  same  in  the  transformed  plane; 
however,  the  spatial  derivatives  require  transformation  via  scale  factors. 
Equations  in  the  transformed  plane  (£,»j)  can  be  obtained  in  terms  of  the 
contravariant ,  covariant,  or  physical  velocity  components  via  coordinate 
transformations  (Thompson,  Warsi,  and  Mas tin  1985).  Sheng  recommended  use  of 
contravariant  components.  Contravariant  components  (g1 ,  §},  aJ)  are  normal  to 
coordinate  planes,  and  covariant  components  (4, ,  a^,  aj)  are  tangential  to 
coordinate  lines  (Figure  14-6).  The  three  components  are  identical  in  a 
uniform  Cartesian  coordinate  system  with  square  grid  cells,  but  differ  for 
other  systems.  Model  WICM  employs  contravariant  components  in  the  transfor¬ 
mation  of  the  governing  equations. 

27.  The  flow  rate  components  in  physical  space  (i.e.,  U(i)  and  V(j)) 
are  related  to  the  contravariant  components  (i.e.,  IP  ,  V1  ,  IP  ,  V1)  by  the 
following  equations: 


U(i) 


JZ±lu*  ♦  -?4 

1^1  1^1 


(14-19) 


VU) 


■St v> 


where : 


*1  ♦  y« 

+  y*y< 

**  +  yl 

or 


9 a  = 


9XX  9x2 
9i 1  922. 


and  |g|  is  the  determinant  of  the  metric  tensor 


(14-20) 


(14-21) 


(14-22) 


|g|  *  9xx922  ~  9x»92i 


(14-23) 
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Figure  14-6.  Covariant  and  contravariant 
components 


The  surface  slope  terms  are  transformed  as  follows: 


as 


»“*  *  ** 


where  g®  are  inverse  metric  tensor  components: 


(14-24) 


(14-25) 


or 


rU 


1 

xj  *  y,2  -  ixnxi  *  y,y<) 

m  1 

Sraa  “Sai 

i*r 

-  (x^  *  yty,)  xl  ♦  y\ 

i^r 

“?ia  ?ii . 

(14-26) 


sr 


U  = 


(14-27) 


For  details  of  these  transformations,  the  reader  is  referred  to  Thompson, 
Warsi,  and  Mastin  (1985). 
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Governing  Equations  in  General  Curvilinear  Coordinates 


28.  The  transformed  governing  equations  developed  by  Sheng  (1986)  are 
as  follows: 

£ -Momentum 


du 

dt 


+  Inertia *  ♦  gH 


Kf 


(14-28) 


V  +  2C,  H, 

p  H  T  sinh(ieH) 


\9u\ 


STli 

"I^iT 


V  ) 


+  Diffusion * 


0 


w- Momentum 
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-=*  +  ( "fir  ^  +  -^t-7  >  +  Diffusion *  *  0 

p  77  r  sinh(ieff)  IsrJ  |g-„l 


Continuity 


i  a 

Isr.l 


(bvl^)  =  o 


(14-30) 


where 

U,V  -  contravariant  unit  flow  rate  components  in  the  transformed 
plane  (superscripts  have  been  dropped  for  convenience) 

g'  -  inverse  metric  tensor  components  (where  i  and  j  are  1  or  2) 

gi}  -  metric  tensor  components  (where  i  and  j  are  1  or  2) 

|gu|  “  determinant  of  the  metric  tensor,  |g|  ,  at  a  U- face** 

7  -  average  y-direction  unit  flow  rate  at  a  17-face** 


rrf,  -  external  surface  shear  stresses,  such  as  wind  stress,  in  the 

£-  and  q-directions ,  respectively 


r»t  •  Tw, 


U 


wave  stress  (gradients  in  radiation  stress),  in  the 
£  -  and  »j -directions ,  respectively 

determinant  of  the  metric  tensor,  |g|  ,  at  a  7- face** 

average  x-direction  unit  flow  rate  at  a  7-face** 
determinant  of  the  metric  tensor,  \g\  ,  at  an  S-point** 


Inertia  and  diffusion  terms  are  given  in  Appendix  14- A. 
S-,  17- ,  7-positions  are  defined  in  Figure  14-7. 
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Figure  14-7.  Definition  of  variable  positions 
Finite  Difference  Methods 

29.  Although  analytic  solutions  of  the  governing  equations  exist  for 
simple  situations,  they  do  not  exist  for  situations  generally  encountered  in 
the  field.  Therefore,  it  is  necessary  to  use  numerical  approximations  of  the 
governing  equations  to  produce  a  general-purpose  model.  This  can  be  accom¬ 
plished  by  representing  the  derivatives  presented  in  Equations  14-28  through 
14-30  with  finite  difference  approximations. 

30.  In  the  finite  difference  approach,  the  continuous  problem  domain  is 
discretized  so  that  dependent  variables  are  defined  only  at  discrete  points. 
Derivatives  are  approximated  by  differences  resulting  in  an  algebraic  problem. 
Several  factors  determine  whether  the  solution  obtained  is  a  good  approxima¬ 
tion  to  the  exact  solution  of  the  original  partial  differential  equation  (PDE) 
including:  truncation  error,  numerical  stability,  roundoff  brror,  and 
discretization  error.  These  factors  will  be  discussed  later  in  this  section. 

31.  As  was  previously  mentioned,  the  continuous  problem  domain  is 
represented  with  a  finite  difference  grid  in  the  (x,y)  or  (£>»l)  plane 
(Figure  14-8).  It  should  be  noted  that  the  (x,y)  and  (£,ij)  planes  can  be 
used  interchangeably  here.  The  ( x,y )  plane  is  used  merely  for  explanation 
purposes.  If  a  variable  a(x,y)  is  defined  at  each  point  on  the  grid  as 
a(itsx,jAy)  ,  then  a  difference  equation  can  be  written  in  terms  of  the  general 
point  (i,j)  and  its  neighboring  points: 

atj  -  a(x,y)  aUuJ  »  a<x+A x,y)  atj =  a(x,y*Ay) 

ai-i.i  *  a(x-Ax,y)  -  a(x,y-Ay) 
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Figure  14-8.  Finite  difference  grid 


32.  Recall  that  the  definition  of  the  derivative  of  a  function  a(x,y) 
at  a  point  (x  -  xol  y  -  y0)  is  given  by: 

da|  _  a(x0*Ax.ye)  -  a(x0,y0) 

<&L  Ai"  - £x -  <u'n> 

This  approximation  improves  as  Ax  decreases.  Developing  a  Taylor  series 
expansion  for  a(x+Ax,y)  about  (x^.yj  yields: 

a(x0+Ax,y0)  *  a(x0.y0 )  +  -|||^  Ax*  -J^~|  ( *  *  H.O.T.  (14-32) 

where  H.O.T.  signifies  Higher  Order  £erms. 

A  forward  difference  is  formed  by  rearranging  Equation  14-32: 


dal  _  a(x0*Ax,y„)  -  a(x0,y„)  g>a  (Ax)  .  „  „  „ 


(14-33) 


or  by  simply  using  the  (i, j)  notation,  the  following  is  obtained: 

dal  ,  -  *±1  +  o(Ax)  (14-34) 

i.j  Ax 


where  O(Ax)  is  used  to  indicate  terms  of  order  Ax  .  The  backwards  differ¬ 
ence  can  be  obtained  similarly  as: 

M  .  **./-  a*-hl  ♦  0( Ax)  (14-35) 

ox\ i.j  Ax 
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The  central  difference  representation  Is  as  follows: 


da 

~^i.i 


~  *  0( Ax) 2 


2Ax 


(14-36) 


These  are  only  a  few  examples  of  the  ways  In  which  derivatives  can  be  approxi¬ 
mated.  Primarily  central  differences  are  used  to  approximate  spatial  deriva¬ 
tives  in  the  formulation  of  WICM. 

33.  The  accuracy  of  these  approximations  is  dependent  on  a  number  of 
factors.  Roundoff  error  is  inherent  to  repetitive  computations  that  are 
rounded  to  a  finite  number  of  digits,  such  as  the  continual  solution  of  the 
finite  difference  equations  in  WICM.  Truncation  error  is  defined  as  the 
difference  between  the  PDE  and  the  difference  approximation.  A  finite 
difference  representation  is  said  to  be  consistent  if  the  truncation  error 
vanishes  as  the  grid  is  refined.  A  more  difficult  problem  is  ensuring  thrt 
the  solution  scheme  is  stable.  A  stable  scheme  does  not  permit  errors,  such 
as  truncation  or  roundoff,  to  grow  as  the  calculations  proceed  from  one  time- 
step  to  another.  Discretization  e*ror  is  the  error  in  the  solution  to  the  PDE 
caused  by  replacing  the  continuous  problem  by  a  discrete  one,  and  it  is 
defined  as  the  difference  between  the  exact  solution  of  the  PDE  and  the  exact 
solution  to  the  finite  difference  equation.  In  other  words,  it  is  the  sum  of 
the  truncation  error  and  any  errors  introduced  by  the  treatment  of  the 
boundary  conditions. 

34.  The  finite  difference  approximations  incorporated  into  WICM  are 
based  on  a  Eulerian  system  where  the  velocities  and  water  surface  fluctuations 
are  computed  at  discrete  locations  within  the  flow  field.  A  network  of  grid 
cells  is  used  to  define  the  parameter  locations .  A  representative  grid  cell 
in  computational  space  (£,ij)  is  shown  in  Figure  14-7.  In  this  staggered 
grid,  the  water  surface  fluctuation  is  defined  at  the  cell  center  (i,j)  ,  £- 
direction  unit  flow  rates  (U)  are  defined  at  the  "west"  (i,j)  and  "east" 
(i+l,j)  cell  faces,  and  the  »j-dii  tion  unit  flow  rates  (V)  are  computed  at 
the  "south"  (i,j)  and  "north"  (i,j+l)  cell  faces.  Finite  difference 
approximations  of  the  governing  equations  follow.  Note  that  the  continuity 
equation  is  split  into  two  parts  (Equations  14-39  and  14-40)  to  be  used  in  the 
solution  scheme  described  in  paragraph  35.  The  sum  of  these  equations  is  the 
one  original  continuity  equation. 
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£ -Momentum 


(14-37) 


where 


-  ®2*  /V  -  2a  -  2a ♦  dlFRlOUlj  *  (1-6)  (FRlOUi'i  *  ®i,i  *  0 

IflU  Ifful  P  P 


n  -  previous  time  level 
*  -  intermediate  time  level 

f-1  -  solve  for  this  time  level 

6  -  weighting  factor  between  successive  time  levels 

I  -  inertia 

D  -  diffusion 


FRIC  *  2  cf  i  9ii  +  ffu  V  \ 

H  T  ainh(kH)  |firj  bu|  US  J 

evaluated  at  time  level  n 


Finite  difference  forms  of  the  inertia  and  diffusion  terms  are  given  in 
Appendix  14- A. 


i -Momentum 


Yhls^hl  +  !*.  +  gHg^  1 _  S1-U2.J-U2  j  + 


(14-38) 


*  (1-6)  gHg22|  S*'J  :^-l)  ♦  -|^fU  ♦ 


r  ~ 


2a  -  2a 
p  p 


+  9 (FRIO  V?,7  ♦  (1-6)  (FRIO  Vitj  +  Di^  =  0 
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*  •&}(i!Hr!k)  *  <->  *  <14-39) 

Iffvll  V‘l->  -  V‘A  _  „ 

"ral - Sn - 1  0 


-  s;.j  . ,  m  ( <}■.  -  ,,  „  m  (  k 

~Ti -  'iSrl - *n I  wl— 

\9v\l  Vl.i*x  ~  Vi.i\  _  n 

15jl - Si  I 


ti  -  r‘.i 

At)  ; 


(14-40) 


Rearranging  the  terms  in  Equations  14-37  through  14-40  yields: 

£ -Momentum 

o[-gHg1'^}slltJ  +  (1  +  6  ( FRIO  A  t)  Ujj  ♦  6(sr*Srll-ff  )5i.i  -  B4  (14-41) 


B,  B,  Bj 

where:  B4  represents  all  known  quantities 


O^-gHg33 -^-J S£j-i  +  (1  +  6  ( FRIO  A  t)  Vf*}  +  6[gHg33-^^sffi  «  Bt  (14-42) 


"(^w  Ti)*1  * [1’ Sl1  *  e(-fer  if)01-'  •  *• 


(14-43) 


where:  A4  represents  all  known  quantities 


•HSt*?)*0  * 111  *  e(w-^)^‘‘  ■  * 


(14-44) 
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35.  The  computational  procedure  used  in  W1CM  is  based  on  an  Alternating 
Direction  Implicit  (ADI)  scheme  (Roache  1976).  Using  this  method,  the  £-  and 
17-momentum  equations  are  solved  separately,  and  each  calculation  in  time  is 
made  in  two  stages  (Figure  14-9).  In  the  first  stage,  the  ^-continuity  and  £- 
momentum  equations  are  solved  along  each  row  of  the  grid  to  progress  from  time 
level  n  to  an  intermediate  time  level  *  .  The  £ -direction  unit  flow  rate 
components  and  water  surface  fluctuations  are  solved  implicitly,  and  the  17- 
direction  unit  flow  rate  components  are  supplied  from  time  level  n.  The  £- 
direction  unit  flow  rates  from  this  step  represent  those  at  time  level  n+I  , 
whereas  the  water  surface  fluctuations  are  only  an  approximation  to  those  at 
time  level  n+I  .  The  17 -direction  unit  flow  rate  components  remain  at  time 
level  n  .  In  the  second  stage,  the  »j- continuity  and  17 -momentum  equations  are 
solved  along  each  column  for  the  t)- direction  unit  flow  rates  and  the  water 
surface  fluctuations  at  time  level  n+I  .  ^-direction  unit  flow  rate  compo¬ 
nents  are  supplied  from  the  first  stage  calculations. 

36.  As  shown  in  the  finite  difference  approximations  to  the  governing 
equations,  a  weighting  factor  6  is  used  to  place  the  water  surface  slope  and 
bottom  friction  terms  between  time  levels  n  and  n+I  .  When  the  weighting 
factor  equals  0.0,  these  terms  are  evaluated  at  the  previous  time  level  n 
(explicit  treatment),  whereas  when  the  weighting  factor  equals  1.0,  they  are 
evaluated  at  the  new  time  level  n+1  (implicit  treatment).  Usually  a  value 
between  0.0  and  1.0  is  used.  Tests  have  shown  that  a  value  of  0.55  produces  a 
stable  solution  with  less  damping  of  the  solution  than  occurs  with  a  value  of 
1.0. 
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Figure  14-9.  Computational  procedure  for  ADI  scheme 
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Tridiaeonal  matrix 

37.  As  was  mentioned  earlier,  the  ^-direction  continuity  and  momentum 
and  the  ij-direction  continuity  and  momentum  equations  are  solved  separately. 
For  a  given  segment  (Figure  14-10),  the  ^-continuity  and  ^-momentum  equations 
are  alternately  solved  at  S-  and  0-points ,  respectively.  Similarly,  the  »j- 
continuity  and  q-momentum  equations  are  alternately  solved  at  S-  and  V-points. 
The  equations  corresponding  to  Figure  14-10  are: 


C 

M 

C 

M 

C 


A„Ub  +  Ajj  S(j  +  Ajj  Ui+,j 

B,.,+,Ss  +  Bj.i+iUj+ij  + 

A|,i+iUt+,j  +  Ali+1Si+1j  +  A^i+iUj+jj 

B|,l+2S|  +  |j  +  62,1+2^1+20  + 

^1,1+2^1+20  +  Ay+jS^jj  +  Aji+^Uj+jj 


—  A# 

“  84,1+1 
“  A44+1 
“  84,1+2 
“  A44+2 


where  C:  indicates  the  x- continuity  equation  and  M:  indicates  the  x-momentum 
equation,  or  in  matrix  form: 
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This  system  of  equations  is  in  the  form  of  a  banded  or  tridiagonal  matrix  and 
is  solved  using  the  Thomas  (1949)  algorithm  (Gaussian  elimination  and  back 
substitution  procedures). 
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Figure  14-10.  £ - sweep  solution  segment 


14-21 


38.  In  the  Gaussian  elimination  step,  all  lower  diagonal  elements  (A,'s 
and  Si's)  are  eliminated  and  are  incorporated  into  other  matrix  elements.  The 
system  of  equations  is  thus  reduced  to  an  upper  diagonal  matrix  with  only  two 
unknowns  per  equation.  The  back  substitution  sweep  then  involves  the  solution 
of  the  system  of  equations  for  all  variables,  U{  and  Sj  . 


Boundary  Conditions 


39.  The  matrix  equation  previously  shown  contains  N-2  simultaneous 
equations  containing  N  unknowns .  Boundary  conditions  provide  the  remaining 
information  needed  to  solve  the  system  of  equations.  The  boundary  conditions 
are  classified  into  three  general  categories:  open  boundaries,  land-water 
boundaries ,  and  subgrid  barriers . 

Open  boundary 

40.  Boundaries  classified  as  open  include:  a  seaward  edge  of  a 
computational  grid  defined  as  water,  or  a  channel  (river)  that  flows  across 
any  boundary.  Water  surface  (tide)  levels  or  unit  flow  (river  flow)  rates  can 
be  specified  along  an  open  boundary  as  functions  of  time.  An  open  boundary 
can  also  be  specified  with  a  uniform  flux  (zero  gradient)  boundary  condition. 
Research  has  been  conducted  to  add  a  radiation  boundary  condition  option, 
which  relates  water  levels  and  flow  rates,  into  the  model.  This  option  will 
be  incorporated  in  a  later  release  of  WICM. 

Tidal  boundary 

41.  A  tidal  boundary  condition  (TBC)  allows  for  fluctuations  in  the 
free  surface  and  is,  therefore,  applied  at  an  S-point  (Figure  14-11).  The 
computations  proceed  from  the  first  internal  77- face  to  the  last  S-point.  For 
the  case  of  a  tidal  boundary  on  the  west  (or  south) ,  a  special  formulation  of 
the  momentum  equation  at  the  next  internal  77- face  is  required: 


B2.i*\Ul*\.j  *  Bi.i*lSi*i.J  =  B*.  i*l  ~  B\,i*\Si,j 


(14-45) 


The  sweep  proceeds  from  this  17-point  to  the  easternmost  S-point. 

42.  For  a  tidal  boundary  on  the  east  (or  north),  a  special  formulation 
of  the  momentum  equation  at  the  last  internal  77-face  is  required: 


+  B2,  lmucUimax.J  ~  B*.imtx  ~  Bi. 


(14-46) 


The  sweep  proceeds  from  the  westernmost  S-point  to  this  77-point. 
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TBC  Applied  TBC  Applied 

Figure  14-11.  Tidal  boundary  condition  (TBC) 


River  boundary 

43.  A  river  boundary  condition  (RBC)  allows  flow  into  or  out  of  a  comp¬ 
utational  grid.  Being  a  flow  condition,  the  river  boundary  condition  is 
applied  at  a  U-  or  U-face  (Figure  14-12).  The  solution  proceeds  from  the  next 
internal  S-point  to  the  last  S-point.  For  a  river  boundary  on  a  west  (or 
south)  boundary  a  special  formulation  of  the  A,  in  the  continuity  equation  at 
the  first  S-point  is  required: 

~  Auui  (14-47) 

The  sweep  proceeds  from  this  S-point  to  the  last  S-point. 

44.  For  a  river  boundary  on  an  east  (or  north)  boundary,  a  special 
formulation  of  the  A,  in  the  continuity  equation  at  the  last  S-point  is 
required: 

A\,im*xUim*x.)  +  A1 .  Imt^lmvr.  j  ~  A*.imax  ~  Al.imucUlm*x+\.)  (14-48) 

The  sweep  proceeds  from  the  first  S-point  to  this  last  S-point. 
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R8C  Applied  RBC  Applied 

Figure  14-12.  River  boundary  condition  (RBC) 
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Uniform  flux  boundary  condition 

45.  A  uniform  flux  (or  zero  gradient)  boundary  condition  (UFBC)  implies 
that  flow  is  constant  across  such  a  boundary.  Similar  to  a  river  boundary, 
the  UFBC  is  applied  at  a  U-  or  F-face,  and  the  solution  begins  at  the  next 
internal  S-point  (Figure  14-13).  For  a  UFBC  on  a  west  (or  south)  boundary,  a 
special  formulation  of  the  A«  in  the  continuity  equation  at  the  first  S- 
point  is  required: 

♦  <*ii  -  -  **i  (14-49) 

The  sweep  proceeds  from  this  S-point  to  the  last  S- point. 

46.  For  a  UFBC  on  an  east  (or  north)  boundary,  a  formulation  of  the  A, 
in  the  continuity  equation  at  the  last  S-point  is  required: 

^*1.  imax  +  A3i  iawr)  Uiaax.j  *  *2,  lmv&imut.j  “  *4,i*Nur 

The  sweep  proceeds  from  the  first  S-point  to  this  last  S-point. 


SOLUTION  SMTP 


Land- 

water  boundaries 

47.  With  land-water  boundaries  (i.e.,  the  coastline),  water  velocities 
normal  to  the  land-water  interfaces  are  set  equal  to  zero  in  the  WICM  model. 
This  condition  means  that  the  boundary  is  closed  and  no  flow  is  permitted 
across  such  a  boundary.  Therefore,  this  model  lacks  the  capability  to 
simulate  flooding  and  drying  of  low-lying  areas.  This  capability  will  be 
added  with  a  future  release  of  WICM. 

48.  A  land-water  (closed)  boundary  is  applied  at  a  U-  or  F-face,  and 
the  solution  begins  at  the  next  internal  S-point.  For  a  land-water  boundary 
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on  a  west  (or  sout.h)  boundary,  a  special  formulation  of  the  At  in  the 
continuity  equation  at  the  first  5-point  is  required: 

+  ^31^1*1. j  m  i'*i  (14-51) 

The  sweep  proceeds  from  this  5-point  to  the  last  5-point. 

49.  For  a  land-water  boundary  on  an  east  (or  north)  boundary,  a  special 
formulation  of  the  At  in  the  continuity  equation  at  the  last  5-point  is 
required: 

Al.UmxVlmuc.J  *  "  A«.  imtx  (14-52) 

The  sweep  proceeds  from  the  first  S-point  to  this  last  S-point. 

Subgrid  barriers 

50.  Subgrid  barriers  (e.g.,  breakwaters)  are  treated  as  thin,  imperme¬ 
able  walls  that  extend  above  the  water  surface,  preventing  flow  from  one  cell 
to  the  next.  These  exposed  barriers  are  defined  along  specified  interior  cell 
faces  and  can  be  used  to  represent  any  nonovertopping  structure  in  a  study 
area.  They  can  also  be  used  to  represent  small  islands  when  using  a  coarse 
grid.  A  subsequent  release  of  WICM  will  include  submerged  and  overtopping 
barriers. 
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PART  III:  DEFINITION  OF  INPUT  DATA  FORMAT 


51.  The  input  data  set  format  was  designed  to  resemble  the  format 
required  by  the  series  of  models  released  by  the  USAE  Hydrologic  Engineering 
Center.  It  is  the  intent  that  this  structure,  being  familiar  to  Corps 
personnel,  will  reduce  the  time  needed  to  learn  this  system.  The  general 
format  of  the  input  data  set  records,  where  a  record  refers  to  one  line  of 
data,  is  presented  below: 


a.  Each  record  is  divided  into  10  fields  containing  8  columns 
each. 

&.  Field  1  (i.e.,  columns  1  through  8)  contains  a  mnemonic  identi¬ 
fication  label  that  describes  the  purpose  or  function  of  each 
record. 

c.  Fields  2  through  10  contain  data  that  may  be  real,  integer,  or 
character  in  type.  Integers  must  be  right- justified.  Real 
numbers  must  also  be  right- justified  if  the  decimal  point  is 
omitted.  Character  data  do  not  need  to  be  right-  or  left- 
justified. 

d.  Array  data,  such  as  depths,  are  read  with  DO  or  Implied  DO 
loops.  No  label  is  required  for  array  data  records.  However, 
a  general  specification  record,  such  as  BATHSPEC,  which  defines 
bathymetric  attributes,  must  precede  that  array. 


52 .  Spelling  of  record  identification  labels  and  alphanumeric  variables 
is  important.  Misspelled  entries  will  result  in  either  recognized  error 
conditions  that  force  the  model  to  abort  execution  or  bypassing  of  desired 
user-defined  operations,  such  as  saving  time-history  data. 

53.  Certain  records  and  variables  have  been  assigned  default  values  in 
the  model  for  minimizing  input  data  and  computer  resources.  Thus,  not  all 
input  data  records  will  be  needed  for  each  application,  and  only  those  records 
pertinent  to  the  simulation  or  required  by  the  model  should  be  included. 
Default  values  are  representative  of  those  chosen  in  previous  studies  per¬ 
formed  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station  Coastal  Engi¬ 
neering  Research  Center.  Although  these  quantities  may  not  be  applicable  to 
all  studies,  they  can  serve  as  a  guide  when  selecting  alternative  values. 
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54.  Default  values  are  processed  when  the  record  field  corresponding  to 
that  variable  is  blank.  Hence,  the  user  must  be  careful  when  leaving  fields 
blank  in  a  record;  blank  fields  will  not  necessarily  result  in  an  input 
variable  being  assigned  a  value  of  zero.  Input  variables  and  their  respective 
default  values  are  noted  in  Appendix  14-B.  The  following  discussion  pertains 
to  the  general  format  of  the  input  records  given  in  Appendix  14-B. 

55.  Each  record  is  presented  in  a  standardized  tabular  format  and  has 
as  its  heading  the  mnemonic  identification  label  or  name  with  a  brief  descrip¬ 
tion  of  its  function.  Following  its  name,  the  record  has  an  abbreviated  note 
documenting  whether  it  is  required  for  a  simulation.  These  abbreviations  have 
the  following  definitions: 

(Req)  Record  or  variable  is  required  for  each  simulation. 

(Opt)  Record  or  variable  is  optional.  Omitting  this  item 

results  in  either  the  default  value  being  used  or  the 
defined  operation  not  being  performed. 

(C-opt)  Record  or  variable  is  required  if  related  or  parent 
options  have  been  selected. 


For  example,  record  TIMESPEC,  presented  in  Appendix  14-B,  contains  the  note 
(Req) ,  meaning  that  this  record  must  reside  in  the  input  data  set  for  each 
simulation.  Record  CHNGBATH  contains  the  note  (Opt),  meaning  this  record  is 
optional  and  is  used  only  when  changes  to  the  bathymetric  data  are  desired. 
Record  XSTRETCH  contains  the  note  (C-opt),  meaning  this  record  is  required 
only  if  a  stretched  Cartesian  grid  has  been  specified  on  the  GRIDSPEC  record 

56.  Input  variables,  presented  in  column  2  of  each  table,  are  refer¬ 
enced  to  their  respective  record  fields  shown  in  column  1.  Generally,  data 
for  each  variable  occupy  a  single  8-column  data  field.  However,  variables 
assigned  titling  or  formatting  information  can  occupy  several  fields. 

57 .  Variable  attributes  are  presented  in  columns  3  through  6  of  each 
table.  Valid  data  types  are  listed  in  column  3  and  can  be  real,  integer,  or 
alphanumeric.  Abbreviations  presented  in  this  column  are  described  below: 
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Char*16  Alphanumeric  character  string  containing  up 
to  16  characters 

Char* 8  Alphanumeric  character  string  containing  up 

to  8  characters 

Integer  Integer  data 

Real  Real  (floating  point)  data 

58.  Column  4  of  each  table  defines  whether  the  respective  variable  must 

t  gned  a  value.  Abbreviations  listed  in  this  column  have  identical 

meat,  gs  as  those  for  the  records.  Default  values  are  listed  in  column  5.  A 

blank  entry  in  this  column  denotes  that  the  respective  variable  does  not  have 

a  default  value. 

59.  Column  6  of  each  table  lists  the  variables'  permitted  data  type  or 

all  valid  character  strings .  Variables  having  integer  or  real  data  types  are 

specified  with  the  following  notation: 

A  Alphanumeric  values 

+R  Positive  real  values 

R  Positive,  zero,  or  negative  real  values 

+1  Positive  integer  values 

I  Positive,  zero,  or  negative  integer  values 

60.  Variable  definitions  are  listed  in  column  7  of  each  table. 

Variables  whose  quantities  are  unit -dependent  contain  a  reference  to  that 
variable  designating  its  system  of  units.  For  example,  variable  TMAX  is 
assigned  a  value  having  units  defined  by  variable  TUNITS.  Variables  defining 
input  data  units  and  the  record  on  which  they  reside  are  presented  below. 


Definition 

BUNITS 

BATHS  PEC 

Bathymetry/topography  data 

FUNITS 

FUNCTION 

Tidal  and  discharge  boundary  data 

GUNITS 

GRIDSPEC 

Numerical  grid  data 

SUNITS 

GENSPECS 

Model  computations  and  output 

TUNITS 

TIMES PEC 

Time -dependent  variables 

WUNITS 

WINDSPEC 

Wind  velocity  data 
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PART  IV:  DISCUSSION  OF  INPUT  DATA  REQUIREMENTS 


61.  The  types  of .  data  processed  by  UICM  are  extensive  and  encompass  a 
wide  range  of  possible  applications.  Since  each  application  is  unique,  the 
type  of  input  data  required  for  each  specific  study  will  vary.  In  this 
discussion  of  model  input,  data  have  been  divided  into  seven  categories  to 
present  model  capabilities  and  data  requirements.  These  categories  are: 

g.  Model  control  specifications. 

}*.  Grid  description. 

£.  Physical  characteristics. 

«J.  Boundary  conditions. 

£.  Wind- field  specifications. 

£.  Wave  field  specifications. 

g.  Output  specifications. 

62.  Table  14-2  presents  WICM  input  data  records  pertaining  to  each 
category.  A  record  refers  to  one  line  of  data,  and  each  record  begins  with  a 
mnemonic  character  string  to  identify  one  record  type  from  another.  Record 
format  is  defined  in  Part  III.  Detailed  specification  of  each  record  follows. 
While  reading  this  part  of  the  report,  the  user  will  find  it  beneficial  to 
refer  to  Appendix  14-B,  where  descriptions  for  each  input  data  record  are 
given. 


Madsl  Cbhtrgl  gpqsificaupus 

63.  Model  control  parameters  define  how  the  user  controls  the  model's 
execution.  Records  contained  in  this  category  include  GENSPECS,  TIMESPEC, 
STARTUP,  and  ADDTERMS. 

64.  Record  GENSPECS  is  used  to  specify  the  general  title  of  the 
simulation  (TITLE)  and  the  system  of  units  (SUNITS)  used  for  model  computa¬ 
tions  and  displaying  model  results.  Variable  names  are  typically  given  in 
parentheses.  In  addition  to  the  general  title,  other  input  data  records  have 
provisions  for  titles.  Although  this  information  is  optional,  it  can  be  very 
helpful  when  reviewing  a  series  of  simulations.  A  title  should  specifically 
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Table  14-2 


Boundary  conditions 

XBOUNDRY 

YBOUNDRY 

FUNCTION 

CNRECORD 

CONSTIT 

TERECORD 

TFRECORD 
TABELEV 
TAB FLOW 

Wind- field  specifications 

WINDSPEC 

TABWINDS 

Wave  field  specifications 

WAVESPEC 

Output  specifications 

PRWINDOW 

RECGAGE 

RECSNAPS 

XRECRANG 

YRECRANG 
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state  data  set  attributes,  such  as  data  source  or  collection  date,  to  differ¬ 
entiate  it  from  data  used  in  other  simulations. 

65.  Model  output  is  presented  in  either  English  or  metric  units. 
However,  the  user  can  specify  a  different  system  of  units  for  the  input  data. 
For  example,  the  user  can  supply  wind  velocity  data  having  units  of  either 
miles  per  hour,  feet  per  second,  meters  per  second,  or  knots.  WICM  will 
perform  the  necessary  conversion  to  place  the  input  data  into  the  system  of 
units  used  for  computations. 

66.  Record  TIMESPEC  controls  the  processing  of  all  time-related  vari¬ 
ables.  Variable  DT  defines  the  time -step  or  time  interval  between  consecutive 
flow  field  computations  and  has  units  of  seconds.  Variable  TUN1TS  controls 
the  units  of  all  other  time -dependent  variables.  Valid  units  for  TUNITS  are 
hours ,  minutes ,  and  seconds . 

67.  For  meaningful  model  results,  the  selected  time-step  must  be 
consistent  with  the  Courant  criterion  C,  .  The  maximum  permissible  time-step 
size  for  a  given  value  of  C,  can  be  computed  from: 

c  =  VggL-  (14-53) 

r  Ax/ A  t 

where  g  is  the  gravitational  acceleration,  H  is  the  depth  at  that  cell, 

Ax  is  the  dimension  of  the  smallest  grid  cell  within  the  computational 
domain,  and  At  is  the  length  of  the  time -step.  Experience  indicates  that* 
the  value  for  C ,  should  be  less  than  5 .  Larger  values  of  the  Courant  number 
may  be  permissible  at  locations  distant  from  the  area  of  interest  within  the 
computational  domain;  however,  numerical  accuracy  will  be  affected. 

68.  Variable  IT1  defines  the  time-step  number  at  the  start  of  the 
simulation,  and  variable  IT2  defines  the  length  of  the  simulation  (in  time- 
steps).  Typically,  a  simulation  begins  at  time-step  1  and  has  a  duration 
ranging  from  one  to  several  hours.  This  duration  provides  the  hydrodynamic 
model  sufficient  time  in  which  to  develop  accurate  circulation  fields  and 
avoids  the  generation  of  numerical  instabilities  induced  by  instantaneously 
applying  forces  to  a  static  water  basin.  Longer  simulations  are  required  for 
tidal  forcing.  Additional  details  concerning  duration  are  discussed  with 
tidal  boundary  conditions. 

69.  Variable  NFREQ  specif' es  the  frequency  at  which  numerical  gauge 
time-history  data  (i.e.,  water  velocities  and  free  surface  elevations)  are 
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saved  for  subsequent  processing  by  the  post-processing  package  in  the  CMS. 
Numerical  gauge  locations  (locations  for  which  time  histories  are  saved)  are 
selected  with  the  RECGAGE  record,  which  is  discussed  in  the  model  output 
specifications  section.  Variable  ITBRKINC  defines  the  time  interval  for 
saving  field  data  for  use  in  a  follow-up  simulation  (see  STARTUP  record). 

70.  The  STARTUP  record  defines  the  flow- field  conditions  at  the  start 
of  a  simulation.  Normally,  the  initial  flow  field  is  assumed  to  be  static 
(i.e.,  water  surface  elevations  and  unit  flow  rates  are  set  to  zero). 

Although  this  assumption  may  not  be  entirely  correct,  it  is  often  necessary 
because  prototype  data  rarely  have  the  resolution  to  accurately  define  the 
flow  field  throughout  the  model  domain  at  the  start  of  a  simulation.  To  allow 
for  dampening  of  transients  associated  with  starting  from  static  initial 
conditions,  the  simulation  should  begin  well  before  the  time  period  of 
interest. 

71.  Initial  conditions  can  also  be  supplied  from  a  previous  simulation 
(SELEV  -  HOTSTART) .  These  simulations  must  be  sequential  in  time  since  the 
second  simulation  is  a  continuation  of  the  first.  Uhen  initial  conditions  are 
supplied  from  a  previous  simulation,  the  second  simulation  is  referred  to  as  a 
"hotstart."  Slight  computational  inaccuracies  also  exist  for  hotstart 
conditions  due  to  data  roundoff  incurred  by  storing  formatted  data.  These 
precision  errors  cause  transients  that  are  significantly  smaller  than  those 
generated  from  starting  with  static  conditions  and  can  be  minimized  to  a 
greater  extent  by  overlapping  the  simulations.  Overlapping  is  performed  by 
saving  the  hotstart  data  file  approximately  10  to  20  time -steps  before  the  end 
of  the  first  simulation.  The  second  simulation  is  then  run  with  this  hotstart 
file.  When  analyzing  model  results,  the  user  should  ignore  the  repeated 
(overlapping)  time-steps  in  the  second  simulation. 

72.  When  a  hotstart  is  selected,  variable  IT1  on  record  TIMESPEC  must 
be  set  to  the  time-step  when  the  desired  field  data  were  saved.  These  times 
are  printed  in  the  input  data  summary  file  of  the  preceding  run.  The  hotstart 
simulation  will  begin  at  time  IT1  and  end  at  time  IT2  +  IT1.  Therefore,  the 
user  should  not  incorporate  the  simulation  time  of  the  previous  run  into 
variable  IT2.  All  remaining  time  variables  should  reflect  the  change  in 
starting  time.  The  STARTUP  record  also  contains  variable  SECHO.  A  brief  or 
full  report  of  input  data  is  written  to  the  output  file  by  specifying  SECHO 
equal  to  SHORT  or  DETAILED,  respectively. 
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73.  Record  ADDTERMS  must  reside  in  the  input  data  set  if  either  the 
inertial  or  diffusion  terms  are  to  be  computed  in  the  Navier-Stokes  calcula¬ 
tions.  Variable  ADVF1AG  is  set  to  either  NO  for  no  inertial  terms  or  YES  to 
include  the  inertial  terms.  Diffusion  is  set  with  variable  DIFFIAG  equal  to 
either  NO  for  no  diffusion  or  YES  for  diffusion.  The  value  of  the 
diffusion  coefficient  is  defined  by  variable  AH  and  is  treated  as  a  constant 
in  WICM.  The  value  of  the  Coriolis  coefficient  is  defined  by  variable  COR  on 
record  ADDTERMS. 


Grid  Description 

74.  The  study  area  is  defined  in  the  model  via  a  computational  grid. 

The  grid  is  composed  of  rectilinear  or  curvilinear  cells,  where  each  cell  is 
assigned  a  two-dimensional  index.  The  first  index  (i)  corresponds  to  the 
x-coordinate ,  and  the  second  index  (j)  corresponds  to  the  y-coordinate .  The 
grid  index  system  is  presented  in  Figure  14-14.  All  flow  field  data,  such  as 
depths,  are  assigned  and  referenced  to  their  respective  grid  cells  with  this 
system.  Guidelines  for  developing  grids  are  discussed  in  Appendix  A  of  the 
1991  edition  of  this  manual. 

75.  WICM  permits  uniform  grids;  rectilinear  stretched  grids,  where 
successive  grid  cells  can  smoothly  vary  in  width;  or  general  curvilinear 
grids.  Procedures  for  constructing  a  grid  are  presented  in  documentation  of 
package  CMSGRID  (Appendix  A,  1991  edition  of  this  manual) .  CMSGRID  guidelines 
for  grid  generation  pertaining  to  hydrodynamic  models  are  therefore  useful  for 
WICM  applications. 
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Figure  14-14.  Definition  of  coordinate  system 
and  grid  cell  convention 
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76.  Selection  of  grid  coordinate  systems  is  controlled  by  variable 
GRTYPE  on  record  GRIDSPEC.  A  uniform,  or  constant,  grid  cell  size  is  selected 
by  assigning  the  character  string  RECTANG  to  variable  GRTYPE,  whereas  string 
RSTRETCH  selects  a  stretched  rectilinear  grid.  If  the  stretched  grid  option 
is  selected,  record  GRIDSPEC  must  be  followed  by  a  set  of  XSTRETCH  and 
YSTRETCH  records  (Table  14-3).  These  records  are  generated  by  the  interactive 
grid  generation  program  MAPIT  contained  in  package  CMSGRID  and  can  be  placed 
directly  into  the  input  data  set  without  modification.  The  reader  can  refer 
to  Appendix  A  in  the  1991  edition  of  this  manual  for  more  information  on 
XSTRETCH  and  YSTRETCH. 

77.  A  general  curvilinear  grid  is  selected  by  assigning  the  character 
string  CURVILIN  to  variable  GRTYPE.  If  the  curvilinear  grid  option  is 
selected,  record  GRIDSPEC  must  be  followed  by  a  GRIDCORN  record.  This  record 
specifies  the  format  of  the  x-  and  y- coordinate  (XCT  and  YCT)  data  produced  by 
the  grid  generation  program  EAGLE;  the  data  directly  follow  the  GRIDCORN 
record.  (The  CMSGRID  package  will  contain  the  EAGLE  software.) 

78.  Variable  GUNITS  on  record  GRIDSPEC  controls  the  system  of  units  for 
the  computational  grid.  Valid  units  are  feet  and  meters.  WICM  will  convert 
the  data  to  the  system  of  units  for  computations  (SUNITS)  internally. 

Variables  ICELLS  and  JCELLS  specify  the  number  of  grid  cells  in  the  x-  and 
y-directions ,  respectively.  Variables  A1XREF  and  ALYREF  on  record  GRIDSPEC 
specify  a  grid's  overall  length  in  the  x-  and  y-directions,  respectively. 

The  variable  XMAP  is  assigned  the  value  of  the  map  scale  if  relative  map 
distances  were  used  in  generating  the  grid.  For  example,  a  map  with  a  scale 
of  1:80,000  is  used  to  generate  a  stretched  rectangular  grid.  If  the  physical 

Table  14-3 

Sample  XSTRETCH  and  YSTRETCH  Records 


XSTRETCH 

1 

4 

-.696333695 

.696333695 

.933460501 

XSTRETCH 

4 

9 

-6.69031772 

5.80606023 

.277821490 

XSTRETCH 

9 

19 

.143110486 

.710267924 

.770050578 

XSTRETCH 

19 

26 

.167892340 

.730987967 

.809764231 

* STRETCH 

1 

17 

- . 500000000 

. 500000000 

1.00000000 

YSTRETCH 

17 

26 

.177832310 

.755987967 

.823454231 
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map  distances  (i.e.,  distances  computed  using  the  appropriate  scaling)  are 
used  to  generate  the  mapping  coefficients,  then  XMAP  is  1.0  because  no 
conversion  is  necessary.  However,  if  the  map  distances  (i.e.,  grid  distances 
measured  in  map  inches)  are  used  to  generate  the  mapping  coefficients,  then 
XMAP  must  be  6666.67  to  convert  map  inches  to  physical  distances  (e.g.,  in 
feet) . 


Physical  Characteristics 

79.  A  study  area's  physical  characteristics  include  its  (a)  topography/ 
bathymetry,  (b)  bottom  friction  coefficients,  and  (c)  if  applicable,  any 
barriers  or  obstructions  influencing  tidal  circulation  or  storm  surge  levels. 
Records  pertaining  to  each  of  the  above  topics  are  described  individually  in 
the  following  sections. 

Topography /bathymetry 

80.  Each  grid  cell  must  be  assigned  a  water  depth  or  land  elevation. 
Topography /bathymetry  data  are  referenced  to  an  arbitrary  datum.  Typically, 
the  map  datum  from  which  the  depths  are  taken  is  used.  Water  cells  are 
designated  by  positive  values,  whereas  land  cells  are  zero  or  negative  values. 

81.  One  BATHSPEC  record  is  required  for  defining  the  general  charac¬ 
teristics  of  the  topography /bathymetry  array  and  must  precede  this  array. 
Variable  BUNITS  defines  the  units  of  topography /bathymetry  data.  Valid  units 
are  feet,  meters,  or  fathoms.  The  input  sequence  for  reading  this  array  is 
controlled  by  variable  BSEQ.  Eight  options  for  the  input  sequence  are 
available  for  reading  the  array  data  and  are  documented  in  Table  14-4.  As  an 
example,  for  the  first  input  sequence  (Figure  14-15),  the  depths  are  read 
along  the  x-direction,  then  y  is  incremented  to  a  value  of  2,  and  again  the 
sweep  in  the  x-direction  takes  place.  This  procedure  is  repeated  until  the 
entire  array  is  read.  The  input  format  for  reading  this  array  can  be  speci¬ 
fied  by  the  user  with  variable  BFORM. 

82.  The  maximum  water  depth  is  specified  with  variable  DLIMIT,  and  any 
array  values  deeper  than  DLIMIT  are  set  to  DLIMIT  (in  BUNITS).  Grid-wide 
adjustments  to  land  elevations  contained  in  the  topography /bathymetry  array 
can  be  made  with  variable  LDATUM.  The  value  assigned  to  this  variable  is 
subtracted  from  all  land  cells  in  the  grid.  Positive  LDATUM  values  will 
increase  land  elevations,  whereas  negative  values  will  decrease  land  eleva- 
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Table  14-4 


Input  Sequence  for  Array  Data 


Ha 

Seauence 

Description 

l 

XY 

1 

DO  1  J-1,J CELLS 

READ (LUN, BFORM)  (VAR(I , J) , 1-1 , ICELLS) 

2 

-XY 

2 

DO  2  J— 1 , JCELLS 

READ(LUN , BFORM)  (VAR(I , J) , I-ICELLS ,1,-1) 

3 

X-Y 

3 

DO  3  J-JCELLS ,1,-1 

READ (LUN, BFORM)  (VAR(I ,J) , 1-1 , ICELLS) 

4 

-X-Y 

4 

DO  4  J-JCELLS, 1,-1 

READ (LUN, BFORM)  (VAR(I ,J) , I-ICELLS , 1 , -1) 

5 

YX 

5 

DO  5  1-1, ICELLS 

READ (LUN, BFORM)  (VAR(I ,J) ,J-1 .JCELLS) 

6 

-YX 

6 

DO  6  1-1, ICELLS 

READ (LUN, BFORM)  (VAR(I ,J) .J-JCELLS , 1 , -1) 

7 

Y-X 

7 

DO  7  I-ICELLS, 1,-1 

READ (LUN, BFORM)  (VAR(I ,J) ,J-1 .JCELLS) 

8 

-Y-X 

8 

DO  8  I-ICELLS, 1,-1 

READ (LUN, BFORM)  (VAR(I ,J) .J-JCELLS , 1 , -1) 

■  *»  x 

Figure  14-15.  Input  sequence  -  option  1 
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Cions.  Similarly,  grid-wide  adjustments  to  water  depths  can  be  made  with 
variable  WDATUM.  The  value  assigned  to  this  variable  is  added  to  all  water 
cells.  Since  these  cells  have  positive  values,  positive  WDATUM  values  produce 
deeper  basins,  and  negative  WDATUM  values  produce  shallower  basins. 

83 .  Changes  to  the  topography /bathymetry  array  can  also  be  made  to 
individual  cells,  or  a  group  of  cells  with  record  CHNGBATH.  This  record 
allows  the  user  to  quickly  change  values  assigned  to  the  bathymetry  array 
(using  variable  BATH)  without  editing  the  array  itself.  It  should  be  noted 
that  (a)  values  of  the  variable  BATH  on  the  CHNGBATH  record  are  assumed  to 
have  units  consistent  with  those  selected  for  bathymetry/topography  (i.e., 
variable  BUNITS  on  record  BATHSPEC) ,  and  (b)  LDATUM  and  WDATUM  are  not  applied 
to  cells  specified  with  record  CHNGBATH;  therefore,  the  effect  of  nonzero 
LDATUM  and  WDATUM  must  be  included  in  the  value  of  variable  BATH. 

84.  Variables  X1INDX  and  X2INDX  on  the  CHNGBATH  record  specify  the 
minimum  and  maximum  cell  numbers  in  the  x-direction,  respectively,  where  the 
bathymetry/topography  value  will  change.  Similarly,  variables  Y1INDX  and 
Y2INDX  on  the  CHNGBATH  record  specify  the  minimum  and  maximum  cell  numbers  in 
the  y- direction,  respectively,  where  the  bathymetry/topography  value  will 
change.  More  than  one  CHNGBATH  record  is  permitted. 

Bottom  friction  coefficient 

85.  Several  choices  are  available  for  specifying  bottom  friction 
coefficients.  Bottom  friction  coefficients  may  be  specified  as  constant  C, 
values  (FRDEF  -  CONSTANT) ,  as  individual  values  for  each  cell  (FRDEF-TABULAR) , 
or  as  functions  of  depth/elevation  (FRDEF  -  VARYBATH) .  For  the  constant  value 
option,  separate  friction  coefficients  can  be  assigned  to  land  and  water 
cells.  Friction  coefficients  for  land  are  assigned  to  variable  FRLAND,  and 
coefficients  for  water  cells  are  specified  by  variable  FRWATR.  A  coefficient 
is  assigned  to  each  cell  based  on  its  topography /bathymetry  value  at  the  start 
of  a  simulation  and  will  remain  constant  throughout  the  simulation. 

86.  The  tabular  friction  option  (FRDEF-TABULAR)  requires  friction 
coefficients  for  each  individual  cell.  The  friction  array  must  follow 
immediately  after  the  FRICTION  record.  The  input  sequence  for  reading  the 
friction  array  is  controlled  by  variable  FSEQ  on  record  FRICTION.  Eight 
options  for  the  input  sequence  are  available  for  reading  the  array  data  and 
are  documented  in  Table  14-4.  The  input  format  for  reading  the  friction  array 
can  be  specified  by  the  user  with  variable  FFORM  on  record  FRICTION. 
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87.  The  depth-variable  friction  option  requires  at  least  two  FRICTABL 
records,  which  must  follow  immediately  after  the  FRICTION  record  in  the  input 
data  set.  Each  FRICTABL  record  must  contain  a  depth/elevation  (FDEPTH)  and  a 
corresponding  friction  value  (FRICT)  for  creating  a  tabular  friction  function. 
Cells  with  a  water  depth  greater  than  or  equal  to  FDEPTH  are  assigned  a 
friction  value  of  FRICT.  Cells  whose  depths  are  deeper  than  the  deepest  depth 
entered  in  the  table  are  assigned  the  friction  coefficient  corresponding  to 
the  deepest  table  depth. 

88.  In  deep  water,  bottom  friction  has  a  negligible  effect  on  long-wave 
hydrodynamics.  Variable  FDMAX  on  record  FRICTION  defines  the  maximum  depth  at 
which  bottom  friction  influences  hydrodynamics.  It  is  provided  for  reducing 
computer  memory  requirements  when  the  depth-variable  option  is  requested  and 
has  a  default  value  of  300. 0  ft. 

89.  Friction  coefficients  can  also  be  assigned  to  individual  cells,  or 
a  group  of  cells  with  record  CHNGFRIC.  Record  CHNGFRIC  overrides  the  above 
options  for  all  cells  entered  on  this  record.  This  feature  permits  the  user 
to  change  a  cell's  friction  coefficient  in  certain  areas,  such  as  rivers, 
where  parameters  differ  from  the  norm.  Variables  (X1INDX,  YlINDX)  and 
(X2INDX,  Y2INDX)  define  the  "patch"  of  cells  where  friction  values  are  to  be 
changed.  More  than  one  CHNGFRIC  can  be  used  in  a  simulation.  These  records 
must  follow  FRICTION  and,  if  applicable,  FRICTABL  records. 

Barriers 

90.  Subgrid  barriers  refer  to  flow  field  obstructions,  such  as  breakwa¬ 
ters,  whose  widths  are  much  narrower  than  the  widths  of  adjacent  grid  cells, 
but  have  lengths  that  are  equal  to  or  greater  than  those  of  a  grid  cell. 
Barriers  are  treated  as  exposed,  impermeable  structures  in  model  WICM. 

91.  Barriers  that  prevent  or  impede  flow  in  the  x-direction  (i.e., 
barriers  perpendicular  to  the  x-axis)  are  specified  with  XBARRIER  records 
(Figure  14-16).  A  barrier  for  a  particular  cell  is  defined  at  the  cell  face, 
where  the  velocity  is  also  defined.  Barrier  grid  locations  are  defined  with 
variables  BRP0S1,  BRP0S2,  and  BRP0S3.  Variable  BRP0S1  is  assigned  the 
barrier's  grid  cell  X- index,  and  variables  BRP0S2  and  BRP0S3  are  assigned  the 
barrier's  starting  and  ending  grid  cell  Y- indices,  respectively. 
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Figure  14-16.  Definitions  of  XBARRIER  and  YBARRIER 

92.  Similarly,  YBARRIER  records  define  those  barriers  that  prevent  or 
impede  flows  in  the  y-direction.  Variable  BRP0S1  is  assigned  the  barrier's 
grid  cell  Y- index.  Variables  BRP0S2  and  BRP0S3  are  assigned  the  barrier's 
starting  and  ending  grid  cell  ^-indices. 


Boundary  Conditions 

93.  WICM  is  generally  applied  in  regions  where  three  grid  boundaries 
lie  in  open  water  (e.g.,  open  coast  regions).  In  these  cases,  the  flow  field 
is  typically  driven  or  forced  by  prescribing  a  time  series  of  water  surface 
levels  along  this  boundary  throughout  a  simulation.  Regions  may  also  contain 
a  river  whose  flow  influences  local  circulation  patterns.  These  river  flow 
rates  can  be  specified  as  a  function  of  time. 

94.  To  encompass  as  many  modeling  scenarios  as  possible,  three  types  of 
open  boundary  conditions  are  incorporated  into  this  model,  including: 

fi.  Tidal  elevation. 

£>.  Discharge. 

£.  Uniform  flux  condition. 

Discussion  of  these  open  boundary  conditions  is  presented  in  three  parts. 
First,  an  overview  of  each  type  of  open  boundary  describes  its  application  and 
attributes.  Second,  input  data  requirements  for  defining  a  boundary's  grid 
location  are  presented.  Third,  input  data  requirements  for  entering  time 
series  data  are  described. 
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Overview  of  open  boundary  conditions 

95.  tfith  the  tidal  elevation  boundary  condition,  the  flow  field  is 
driven  by  specifying  the  water  surface  elevations  as  a  function  of  space  and 
time  along  a  grid's  outer  boundary.  Water  surface  elevations  are  measured 
relative  to  the  topography /ba thyme try  datum  and  can  be  supplied  from  tidal 
constituent  parameters  or  tabulated  data. 

96.  Two  options  are  available  for  specifying  the  spatial  distribution 
of  data  along  a  boundary.  The  first  option  applies  a  spatially  constant 
elevation  function  across  a  boundary  segment  where  each  cell  in  the  segment 
acts  in  unison,  having  identical  water  surface  elevations  at  each  time-step. 
This  option  is  normally  applied  along  a  grid’s  lateral  boundary  (i.e., 
perpendicular . to  the  shoreline)  in  areas  of  deep  water.  It  can  also  be  used 
along  a  seaward  boundary  if  no  significant  phase  or  elevation  differences 
occur  along  this  segment. 

97.  The  second  option  is  a  spatially  variable  elevation  forcing  boun¬ 
dary  where  separate  time-series  functions  are  specified  at  each  end  of  a 
boundary  segment.  Elevations  for  each  cell  within  the  boundary  segment  are 
obtained  by  linearly  interpolating  the  two  forcing  functions  (at  each  time- 
step)  .  This  option  is  typically  applied  along  the  grid's  seaward  boundary.  A 
boundary  can  be  treated  as  one  continuous  segment,  or  as  several  segments. 

98.  Discharge  boundary  conditions  can  be  specified  at  any  interior  grid 
cell  face.  An  error  message  will  be  printed  and  execution  terminated  if  this 
boundary  condition  is  specified  at  the  edge  of  the  computational  grid. 
Furthermore,  depths  on  each  side  of  the  boundary  cell  faces  must  be  greater 
than  zero,  signifying  water  cells.  Discharges  are  required  to  have  units  of 
flow  per  unit  width  (i.e.,  cubic  feet/second/foot  or  cubic  meters/second/ 
meter),  and  flow  trajectories  are  referenced  relative  to  the  grid  axes. 
Positive  discharge  rates  define  flows  directed  in  the  positive  x-  or 
y-directions,  whereas  negative  discharge  rates  define  flows  directed  in  the 
negative  x-  or  y-direction. 

99.  A  uniform  flux  boundary  condition  is  similar  to  a  river  discharge 
boundary  condition  in  that  a  flux  is  applied  at  the  boundary  cell  face. 
However,  the  flux  through  the  cell  face  is  computed  internally  by  the  model, 
as  opposed  to  a  user-defined  flux.  This  boundary  condition  assumes  that  the 
flux  through  the  boundary  cell  face  equals  the  flux  through  the  opposing  cell 
face. 
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100.  The  uniform  flux  boundary  condition  is  the  most  common  lateral 
boundary  condition  for  wave- induced  current  applications.  There  are  three 
situations  where  it  can  be  applied.  First,  it  can  be  applied  across  the 
entire  lateral  boundaries  for  current  fields,  in  relatively  shallow  water, 
driven  within  the  model  domain  by  wave  forcing. 

101.  Second,  uniform  flux  boundary  conditions  can  be  applied  across  the 
lateral  boundaries  near  the  shoreline  (Figure  14-17).  Tidal  signals  generated 
from  constituents,  which  are  normally  prescribed  as  a  forcing  function  along 
the  seaward  boundary,  should  not  be  used  in  this  reach  because  bottom  friction 
effects  attenuate  the  tidal  signal.  Hence,  forcing  a  tidal  signal  in  this 
reach  may  not  accurately  depict  the  physical  processes  occurring  within  the 
study  area.  Third,  the  uniform  flux  boundary  condition  may  be  applicable  to 
steady-state,  storm  surge  analyses  of  lakes.  However,  wind  must  be  the  only 
driving  force  causing  surge,  and  resulting  seiches  cannot  be  simulated  with 
this  boundary  condition. 

102.  Because  the  uniform  flux  boundary  condition  provides  an  estimate 
of  the  actual  flow  entering/exiting  the  grid,  its  use  will  introduce  error 
into  the  flow  field  solution.  Bottom  friction  effects  can  dampen  these 
inaccuracies,  ensuring  accuracy  of  results  in  the  area  of  interest,  if  the 
uniform  flux  boundaries  are  placed  a  suitable  distance  away  from  the  area  of 
interest.  Generally,  this  distance  ranges  from  1.5  to  2.0  times  that  distance 
measured  from  the  shoreline  to  the  seaward  boundary. 


Shoreline 


Figure  14-17.  Uniform  flux  boundary  condition 
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Spatial  description  of  open  boundaries 

103.  Each  distinct  open  boundary  must  be  declared  by  either  an  XBOUNDRY 
or  a  YBOUNDRY  record.  Record  XBOUNDRY  defines  a  boundary  that  induces  or 
forces  flow  in  the  x-direction  (i.e.,  boundary  is  perpendicular  to  the  x-axis) 
(Figure  14-18).  Conversely,  a  YBOUNDRY  record  defines  a  boundary  that  is 
perpendicular  to  the  y-axis. 

104.  Boundary  types  are  declared  on  records  XBOUNDRY  and  YBOUNDRY  with 

variable  BNDTYP  and  can  be  assigned  any  of  the  three  options  previously 
discussed:  (a)  spatially  constant  tide  boundaries  are  selected  with  character 

string  CONSTELV,  (b)  spatially  variable  tide  boundaries  are  chosen  with  string 
INTRPELV,  (c)  spatially  constant  discharge  boundaries  are  selected  with  string 
CONSTDIS,  and  (d)  uniform  flux  boundaries  are  selected  with  string  UNIFLUX. 

105.  This  model  employs  a  staggered  grid  cell  system  where  water 
surface  levels,  x-direction  velocities,  and  y-direction  velocities  are  defined 
at  different  locations  on  a  cell.  These  positions  were  defined  in 

Figure  14-7.  Elevations  are  referenced  relative  to  a  cell's  center,  and  river 
discharge  boundaries  are  specified  along  cell  faces. 

106.  Boundary  grid  positions  are  defined  with  variables  BNP0S1,  BNP0S2, 
and  BNPOS3 .  For  XBOUNDRY  records,  variable  BNP0S1  is  assigned  the  boundary's 
X- index,  whereas  variables  BNP0S2  and  BNP0S3  are  assigned  the  Y- indices  of  the 
starting  and  ending  cells,  respectively  (Figure  14-18).  For  YBOUNDRY  records, 
variable  BRPOS1  is  assigned  the  boundary's  P- index,  and  variables  BRPOS2  and 
BRP0S3  are  assigned  the  X-indices  of  the  starting  and  ending  cells,  respec¬ 
tively.  Boundaries  can  be  identified  by  assigning  each  one  a  boundary  name 
with  variable  BNDNAM  on  the  XBOUNDRY  or  YBOUNDRY  records. 
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Figure  14-18.  Definition  of  XBOUNDRY  and  YBOUNDRY 
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Tine-series  description  of  open  boundaries 

107.  Time-series  daCa  are  required  for  driving  the  flow  field  for  all 
of  the  boundary  options:  (a)  spatially  constant  elevation,  (b)  spatially 
variable  elevation,  and  (c)  spatially  constant  discharge  boundary  options. 
XBOUNDRY  and  YBOUNDRY  variables  BNDFN1  (for  all  three  options)  and  BNDFN2  (for 
spatially  variable  elevation  boundaries  only)  define  which  group  of  time- 
series  data  is  to  be  used  to  drive  the  boundary,  and  correspond  to  the 
function  index  number  assigned  to  variable  FUNNO  on  record  FUNCTION.  For 
example,  an  application  may  have  a  spatially  variable  elevation  boundary  and 
two  spatially  constant  elevation  boundaries  (Figure  14-19).  This  application 
would  require  two  FUNCTION  records  (one  at  each  end  of  the  spatially  variable 
boundary)  to  define  the  time  series  at  two  locations.  These  values  are  also 
used  along  the  spatially  constant  boundaries. 

108.  The  FUNCTION  record  is  required  for  time-series  data  and  defines 
the  character  of  the  forcing  function  driving  the  boundary.  Variable  FUNTYP 
on  record  FUNCTION  is  used  to  determine  if  the  boundary  is  forced  with 
harmonic  constituents  (HARMCNST) ,  tabular  elevation  data  (TABELEVS) ,  or 
tabular  discharge  data  (TABFLOWS).  Variable  FUNITS  defines  the  units  used  for 
the  given  elevations  or  flows.  The  forcing  function  can  be  adjusted  by 
shifting  it  in  time  (phase  shift)  or  in  the  vertical  (datum  shift) ,  or  by 
multiplying  the  entire  time-series  by  an  adjustment  factor.  Variables  FSHIFT, 
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Figure  14-19.  Example  showing  use  of  BNDFN1, 
BNDFN2,and  FUNNO  variables 
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FDATUM,  and  FMULT  on  record  FUNCTION,  respectively,  are  used  to  perform  these 
tasks.  The  forcing  function  can  be  gradually  spline-fit  from  initial 
conditions  (usually  zero)  to  the  given  function  value  over  the  time  period 
assigned  to  variable  FEATHR  defined  on  record  FUNCTION. 

Harmonic  constituents 

109.  A  boundary  driven  with  harmonic  constituents  (FUNTYP  -  HARMCNST) 
requires  at  least  two  records  following  the  FUNCTION  record:  one  CNRECORD 
record  and  at  least  one  CONSTIT  record.  First,  the  CNRECORD  record  specifies 
the  physical  coordinates  and  timing  associated  with  a  constituent  forcing 
function.  The  calendar  date  (variables  RYEAR,  RMONTH,  and  RDAY) ,  time  at  the 
start  of  the  simulation  (variable  RHOUR) ,  and  the  west  longitude  of  the  grid 
cell  where  the  forcing  function  is  applied  (RLONG)  are  specified  on  record 
CNRECORD.  Dates  and  times  can  be  referenced  to  Greenwich  Mean  Time  (GMT)  or 
local  time.  GMT  is  assumed  if  the  longitude  (RLONG)  equals  zero,  and  nonzero 
longitude  values  signal  a  local  time  specification.  CNRECORD  also  contains 
variable  RNAME  for  specifying  a  record  name  for  the  forcing  function. 

110.  One  CONSTIT  record  is  needed  for  each  constituent  used  in  generat¬ 
ing  the  time- series  data.  Variable  CNAME  on  record  CONSTIT  defines  which 
constituents  are  selected.  Table  14-5  provides  a  list  of  the  37  constituents 
available.  A  maximum  of  37  constituents  can  be  selected  (requiring  37  CONSTIT 
records)  and  are  referenced  with  the  naming  convention  developed  by  the 
National  Ocean  Survey.  Variable  CNAME  on  record  CONSTIT  is  used  to  select  any 
of  the  constituents  listed  in  Table  14-5.  Other  data  entered  on  this  record 
include  the  constituent's  amplitude  (CAMP)  and  epoch  (CEPOCH) . 

Tabular  elevation  time  series 

111.  Tabular  elevation  data  can  be  used  to  specify,  or  force,  a 
boundary.  WICM  can  access  components  of  program  TIDEGEN  in  package  CMSUTIL 
(Appendix  B  of  the  1991  edition  of  this  manual)  to  internally  compute  the 
tabular  elevation  time-history  data  from  tidal  constituents.  Each  tabular 
tide  boundary  requires  a  TERECORD  to  precede  the  actual  water  surface  eleva¬ 
tion  data.  The  elevation  data  can  be  given  at  even  (constant)  time  increments 
(RRINT  -  Real  number)  or  at  user-specified  (irregular)  times  (RRINT 

-  IRREGINT). 

112.  Constant  time-interval  data  are  entered  as  a  one -dimensional  array 
immediately  following  record  TERECORD.  The  time  interval  between  entries  is 
specified  with  variable  RRINT  (in  TUNITS)  on  record  TERECORD.  The  total 
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Table  14-5 

List  of  Available  Constituents 


* 


_ Description _ 

M2 

Semidiurnal 

S2 

Semidiurnal 

N2 

Semidiurnal 

K1 

Diurnal 

M4 

Shallow- water 

quarter  diurnal 

01 

Diurnal 

M6 

Shallow-water 

sixth  diurnal 

MK3 

Shallow-water 

terdiumal 

S4 

Shallow -water 

quarter  diurnal 

MN4 

Shallow-water 

quarter  diurnal 

NU2 

Semidiurnal 

S6 

Shallow-water 

sixth  diurnal 

MU2 

Semidiurnal 

2N2 

Semidiurnal 

001 

Diurnal 

LAMBDA2 

Semidiurnal 

SI 

Diurnal 

Ml 

Diurnal 

J1 

Diurnal 

MM 

Long-period 

SSA 

Long-period 

SA 

Long-period 

MSF 

Long -period 

MS 

Long-period 

RH01 

Diurnal 

Q1 

Diurnal 

T2 

Semidiurnal 

R2 

Semidiurnal 

2Q1 

Diurnal 

PI 

Diurnal 

2SM2 

Shallow-water 

semidiurnal 

M3 

Terdiumal 

L2 

Semidiurnal 

2MK3 

Shallow-water 

terdiumal 

K2 

Semidiurnal 

M8 

Shallow-water 

eighth  diurnal 

MS4 

Shallow-water 

quarter  diurnal 
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number  of  tabular  elevation  data  entries  in  the  one -dimensional  array  is 
specified  with  variable  RENT.  The  user  must  also  specify  which  of  the  array 
elements  corresponds  to  the  start  of  the  simulation  with  variable  RSTART. 
Typically,  this  element  is  the  first  value  in  the  array.  However,  this 
variable  must  be  updated  for  hotstart  simulations.  Variable  RFCRM  is  used  to 
specify  the  format  for  reading  the  one -dimensional  array,  and  variable  RNAME 
is  used  for  naming  the  tabular  tide  record. 

113.  For  irregularly  spaced  (in  time)  tabular  tide  data,  one  TERECORD 
and  a  series  of  TA6ELEV  records  are  required.  Variable  RR1NT  on  record 
TERECORD  is  set  to  IRREGINT  for  irregularly  spaced  (in  time)  tide  data.  The 
remaining  variables  on  this  record  are  ignored  with  this  option. 

114.  One  TABELEV  record  is  needed  for  each  tabular  tide  entry.  The 
elevation  (FMAG  in  FUNITS)  and  time  of  occurrence  (FHR,  FMIN,  FSEC  in  hours, 
minutes,  and  seconds)  are  required  on  each  TABELEV  record.  Elevations  are 
linearly  interpolated  at  those  time -steps  occurring  between  entries;  there¬ 
fore,  at  least  two  TABELEV  records  are  required.  The  first  TABELEV  record 
must  specify  the  elevation  at  or  before  the  start  of  the  simulation,  whereas 
the  last  TABELEV  record  must  specify  the  elevation  at  or  after  the  end  of  the 
simulation. 

Tabular  rjvey  discharge.,  time  geyies 

115.  River  discharge  boundary  conditions  can  be  specified  at  any 
interior  grid  cell  face.  An  error  message  will  be  printed  and  execution 
terminated  if  this  boundary  condition  is  specified  at  an  outer-boundary  face. 
Furthermore,  depths  on  opposite  sides  of  the  river  boundary  cell  face  must  be 
greater  than  zero,  signifying  a  water  cell.  Discharges  are  required  to  have 
units  of  flow  per  unit  width  (i.e.  cubic  feet/second/foot  or  cubic 
meters/second/meter) . 

116.  Each  river  discharge  boundary  function  requires  one  TFRECORD 
(analogous  to  the  tabular  tide  function  record  (TERECORD)  discussed  previous¬ 
ly)  to  precede  the  flow  data.  The  flow  data  can  be  given  at  even  (constant) 
time  increments  (RRINT  -  real  number),  or  at  user-specified  (irregular)  times 
(RRINT  -  IRREGINT). 

117.  Constant  time -interval  data  are  entered  as  a  one -dimensional  array 
immediately  following  record  TFRECORD.  The  time  interval  between  entries  is 
specified  with  variable  RRINT  (in  TUNITS)  on  record  TFRECORD.  The  total 
number  of  tabular  flow  data  entries  in  the  one -dimensional  array  is  specified 
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with  variable  RENT.  The  user  must  also  specify  which  of  the  array  elements 
corresponds  to  the  start  of  the  simulation  with  variable  RSTART.  Typically 
this  element  is  the  first  value  in  the  array.  However,  this  variable  must  be 
updated  for  hotstart  simulations.  Variable  RFORM  is  used  to  specify  the 
format  for  reading  the  one -dimensional  array,  and  variable  RNAME  is  used  for 
naming  the  tabular  flow  record. 

118.  For  irregularly  spaced  (in  time)  tabular  flow  data,  one  TFRECORD 
and  a  series  of  TABFLOW  records  are  required.  Variable  RR1NT  on  record 
TFRECORD  is  set  to  IRREGINT  for  irregularly  spaced  (in  time)  flow  data.  The 
remaining  variables  on  this  record  are  ignored  with  this  option. 

119  One  TABFLOW  record  is  needed  for  each  tabular  flow  entry.  The 
flow  rate  (FMAG  in  FUNITS)  and  time  of  occurrence  (FHR,  FMIN,  FSEC  in  hours, 
minutes,  and  seconds)  are  required  on  each  TABFLOW  record.  Flow  rates  are 
linearly  interpolated  at  those  time-steps  occurring  between  entries;  there¬ 
fore,  at  least  two  TABFLOW  records  are  required.  The  first  TABFLOW  record 
must  specify  the  flow  rate  at  or  before  the  start  of  the  simulation,  whereas 
the  last  TABFLOW  record  must  specify  the  flow  rate  at  or  after  the  end  of  the 
simulation. 

Uniform  flux  boundary  condition 

120.  The  uniform  flux  boundary  condition  is  similar  to  a  river  dis¬ 
charge  boundary  condition  in  that  a  flux  is  applied  at  the  boundary  cell  face. 
However,  the  flux  through  the  cell  face  is  computed  internally  by  the  computer 
as  opposed  to  specifying  the  flux  for  a  river  boundary.  This  boundary 
condition  assumes  that  the  flux  through  the  boundary  cell  face  is  equal  to  the 
flux  through  the  opposing  cell  face. 

121.  Since  the  UFBC  does  not  drive  a  flow  field,  its  application  is 
limited.  Therefore,  care  must  be  taken  in  applying  it.  Selecting  an  inappro¬ 
priate  reach  (length)  of  the  boundary  may  lead  to  erroneous  results.  The  user 
should  test  the  sensitivity  of  choice  by  trying  other  segment  lengths  that 
band  the  chosen  reach. 


Wind-Field  Specifications 

122.  WICM  can  simulate  hydrodynamics  influenced  by  uniform,  non- 
uniform,  steady,  and  time-varying  wind  fields.  However,  wind-field  calcula¬ 
tions  are  only  performed  if  record  WINDSPEC  and  related  records  reside  in  the 
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input  data  set.  Uniform  (spatially  and  temporally  constant)  wind  fields  are 
specified  by  assigning  character  string  UNIFSTRS  or  UNIFSPED  to  variable 
WNTRVL  on  record  WINDSPEC.  The  WINDSPEC  record  is  followed  by  one  TABWINDS 
record  specifying  the  x-  and  y-components  of  constant  wind  stress  or  speed 
(Table  14-6). 

123.  Nonuniform  (spatially  variable)  wind  fields  are  specified  by 
assigning  character  string  SPVRSTRS  or  SPVRSPED  to  variable  WNTRVL  on  record 
WINDSPEC.  The  WINDSPEC  record  is  followed  by  one  TABWINDS  record  specifying 
the  sequence  (TSEQ)  and  format  (TWFORM)  of  the  wind- stress  or  wind- speed 
arrays,  similar  to  the  bathymetry  sequence  (BSEQ)  and  format  (BFORM)  (see 
Figure  14-15).  The  TABWINDS  record  is  then  followed  by  the  wind-stress  or 
wind- speed  arrays. 

124.  Time -variable,  spatially  constant  wind  fields  are  specified  by 
assigning  character  string  TVRUNISP  or  TVRUNIST  to  variable  WNTRVL  on  record 
WINDSPEC.  The  WINDSPEC  record  must  be  followed  by  at  least  two  TABWINDS 
records  that  specify  the  x-  and  y- components  of  the  wind  stress  or  speed  at 
specific  times  during  the  simulation. 

125.  Wind  fields  that  vary  in  time  and  space  are  specified  by  assigning 
character  string  TSPVRSTR  or  TSPVRSPD  to  variable  WNTRVL  on  record  WINDSPEC. 
The  WINDSPEC  record  is  followed  by  at  least  two  TABWINDS  records,  and  each 
TABWINDS  record  is  followed  by  wind-stress  or  wind-speed  arrays. 

126.  For  the  time-variable  uniform  wind  options,  one  TABWINDS  record  is 
typically  entered  for  each  hour  of  the  simulation.  (More  frequent  data  rarely 
exist.)  The  model  will  then  linearly  interpolate  to  obtain  wind-speed  and 
direction  values  at  the  time  interval  specified  by  WINTRP  on  record  WINDSPEC. 
It  should  be  noted  that  WINTRP  must  be  assigned  a  value  that  is  an  integer 
multiple  of  the  time -step  DT. 

127.  Variable  WUNITS  declares  the  system  of  units  for  wind  magnitudes. 
Valid  units  are:  (a)  miles/hour,  selected  with  character  string  MPH; 

(b)  feet/second,  chosen  with  FPS;  (c)  meters/second,  specified  with  MPS;  and 
(d)  knots,  selected  with  KNOTS. 

128.  Strong  winds  instantaneously  applied  to  a  static  water  basin,  such 
as  at  the  start  of  a  simulation,  will  generate  extraneous  water  oscillations 
that  can  lead  to  the  model  becoming  unstable  or  can  remain  in  the  model 
solution  for  a  significant  portion  of  a  simulation.  These  problems  can  be 
alleviated  by  gradually  increasing  wind  magnitudes  from  zero  at  the  start  of 
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Table  14-6 

Options  for  Wind-Field  Specifications 


Cases  1  and  2  Constant  wind  field  with  respect  to  time  and  space 


WINDSPEC 

UNIFSTRS 

MPH 

WINTRP 

WIND  DATA  NAME 

or 

UNIFSPED 

TABWINDS 

TAUX 

TAUY 

Cases  3  and  4  Spatially  variable  wind  field 


WINDSPEC 

SPVRSTRS 

MPH 

WINTRP 

or 

SPVRSPED 

TABWINDS 

WIND  DATA  NAME | 
TSEQ 


TWFORM 


(Array  of  TX1) 
(Array  of  TY1) 


Cases  5  and  6  Time -variable  wind  field 


WINDSPEC 

TVRUNISP 

MPH 

WINTRP 

or 

TVRUNIST 

TABWINDS 

IDAY(l) 

IHOUR(l) 

TAUX(l) 

TAUY(l) 

TABWINDS 

IDAY(2) 

IHOTJR(2) 

TA\JX(2) 

TAUY(2) 

TABWINDS 

IDAY ( 3 ) 

IHOUR(3) 

TAUX(3) 

TAUY(3) 

WIND  DATA  NAME 


Cases  7  and  8  Time  and  spatially  variable  wind  field 


WINDSPEC 

TSPVRSTR 

MPH 

WINTRP 

or 

TSPVRSPD 

TABWINDS 

IDAY(l) 

IHOUR(l) 

(Array  oi 

f  TX1) 

(Array  of  TY1) 

TABWINDS 

IDAY(2) 

IH0UR(2) 

(Array  of  TX1) 

(Array  of  TY1) 

WIND  DATA  NAME 


TSEQ 


TWFORM 


the  simulation,  to  the  actual  magnitude  several  hours  into  the  simulation. 
This  task  can  be  accomplished  by  adjusting  wind  magnitudes  entered  on  the 
TABWINDS  records . 

129.  All  TABWINDS  records  must  immediately  follow  the  WINDSPEC  record. 
Furthermore,  times  entered  on  these  records  must  correspond  to  the  simula- 
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tion's  starting  time  defined  by  variable  IT1  on  record  TIMESPEC.  The  first 
TABVINDS  record  must  occur  no  later  than  the  start  of  the  simulation.  For  the 
time-variable  option,  the  last  TABVINDS  record  must  occur  no  earlier  than  the 
end  of  the  simulation. 

130.  For  example,  if  the  simulation  begins  at  0  hr  0  min  and  has  a 
duration  of  36  hr,  the  first  TABUINDS  record  must  specify  wind  conditions 
occurring  at  (or  before)  0  hr  0  min,  and  the  last  record  must  occur  at  36  hr 
0  min  (or  later) .  Should  the  simulation  be  continued  with  the  hotstart 
option,  all  times  entered  on  the  TABVINDS  record  are  referenced  to  the 
original  starting  time.  Assuming  the  hotstart  begins  at  hour  36  and  has  a 
duration  of  24  hr,  the  first  TABUINDS  record  must  have  a  time  no  later  than 
36  hr  0  min.  The  last  TABUINDS  record  must  have  a  time  of  60  hr  0  min  (or 
later) . 


Uave-Field  Specifications 

131.  UICM  can  simulate  hydrodynamics  driven  by  steady,  nonuniform, 
monochromatic  or  spectral  wave  fields.  Uave  stress  calculations  are  performed 
only  if  record  UAVESPEC  resides  in  the  input  data  set.  Uave  fields  are 
specified  by  assigning  character  string  RCPVAVE  or  STUAVE  to  variable  VAVMODEL 
on  record  UAVESPEC.  VAVMODEL  specifies  the  type  of  wave  input  (monochromatic 
for  RCPVAVE  or  spectral  for  STUAVE).  Values  of  wave  height,  period,  and 
direction  for  each  grid  cell  follow  the  UAVESPEC  record  in  the  format  produced 
by  the  RCPVAVE  or  STUAVE  wave  model.  That  is,  if  variable  VAVMODEL  is 
RCPVAVE,  then  UICM  reads  wave  data  in  RCPVAVE  output  format.  Similarly,  if 
variable  VAVMODEL  is  STUAVE,  then  UICM  reads  wave  data  in  STUAVE  output 
format . 

132.  Strong  wave  stress  applied  to  a  static  water  basin  at  the  start  of 
a  simulation  will  generate  extraneous  water  surface  oscillations  that  can  lead 
to  the  model  becoming  unstable  or  remain  in  the  model  solution  for  a  signifi¬ 
cant  portion  of  a  simulation.  These  problems  can  be  alleviated  by  gradually 
increasing  the  wave  stress  from  zero  at  the  start  of  the  simulation  to  the 
actual  magnitude  over  several  time-steps.  This  task  is  accomplished  by 
specifying  the  number  of  time -steps  over  which  to  build  the  wave  stress  to  its 
actual  value  with  variable  IBUILD. 
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Output  Specification 


133.  U1CM  provides  several  options  for  displaying  output,  including: 

(a)  numerical  gage  time -histories  at  selected  grid  cells  (e.g.,  elevations, 
wind,  and/or  water  velocities),  (b)  wind  or  water  velocity  vector  fields  at 
user-specified  times  during  a  simulation,  and  (c)  an  output  listing  containing 
an  input  data  summary  and  a  printout  of  field  arrays  (e.g.,  free  surface 
elevations,  water  ^cities,  wind  velocities,  bathymetry).  CMSMODEL  will 
prompt  the  user  u  rovide  a  name  for  each  specific  output  file  (i.e,  hydro¬ 
graph  file  name  for  RECGAGE  data,  snapshot  file  name  for  field  array  (REC- 
SNAPS)  data).  These  data  files  can  be  used  as  input  to  other  hydrodynamic 
models  as  well  as  wave  hindcast  models  contained  within  the  CMS.  Data  can  be 
displayed  in  either  tabular  or  graphical  form  with  the  post-processing  package 
CMSPOST  discussed  in  Appendix  C. 

134 .  The  output  listing  containing  a  summary  of  the  input  data  set  is 
generated  for  every  simulation.  Diagnostic  error  and  warning  messages  are 
also  contained  in  this  listing.  A  sample  output  listing  containing  a  summary 
of  the  input  data  set  is  presented  in  Figure  14-20. 

135.  Each  input  record  is  summarized  in  tabular  form  with  a  heading 
containing  its  record  identification  label  followed  by  a  brief  description  of 
that  record's  function.  A  table  is  composed  of  each  variable's  name,  a 
description  of  that  variable  (including  its  units,  when  applicable),  and  an 
error  diagnostic  note. 

136.  WICM  contains  error  diagnostic  features  that  inspect  an  input  data 
set  for  possible  errors.  These  features  include:  (a)  comparing  an  input 
value  against  a  range  of  quantities  that  are  representative  for  that  variable, 

(b)  checking  for  misspelled  character  data,  and  (c)  checking  for  missing  data. 
The  error  diagnostic  note  can  be  assigned  one  of  three  character  strings, 
which  are:  (a)  "FATAL"  for  errors  where  the  model  cannot  execute  given  the 
value  supplied,  (b)  "WARN"  for  data  that  are  outside  the  range  of  values 
typically  selected  for  that  variable,  and  (c)  a  null  string  for  instances 
where  an  error  condition  has  not  been  identified.  Although  this  model 
contains  error  diagnostic  capabilities,  the  user  should  thoroughly  inspect  the 
input  data  summary  to  ensure  that  the  data  are  correct. 

137.  Field  arrays  (e.g.,  free  surface  elevations,  water  velocities, 
wind  velocities,  bathymetry)  are  printed  along  with  the  input  data  summary  by 
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COASTAL  MODELING  SYSTEM  (CMS):  UICN  ,  VERSION  1.0 

-  VICM  SIMULATION  NO.  1:  PLANE  BEACH,  H*1.«M,T-4S,A-30DEG  - 

*****  GENSPECS  CARD:  SPECIFICATION  Of  TITLE  AND  GENERAL  SYSTEM  OF  UNITS 
VARIABLE  DESCRIPTION  OF  USAGE:  VALUE:  NOTES: 

UNITS  UNITS  SYSTEM  USED  IN  COMPUTATIONS  METRIC 


*****  GRIDSPEC  CARD:  SPECIFICATION  OF  THE  TYPE  OF  FINITE-DIFFERENCE  GRID  USED 


VARIABLE 

DESCRIPTION  OF  USAGE: 

VALUE:  NOTES: 

GRTYPE 

TYPE  OF  FINITE-DIFFERENCE  GRID 

RECTANG 

GUNITS 

SYSTEM  OF  UNITS  USED  FOR  THE  GRID 

METRIC 

XCELL 

NUMBER  OF  GRID  CELLS,  X  DIRECTION 

60 

YCELL 

NUMBER  OF  GRID  CELLS,  Y  DIRECTION 

8 

ALXREF 

TOTAL  LENGTH  IN  X  DIRECTION 

900. 

ALYREF 

TOTAL  LENGTH  IN  Y  DIRECTION 

240. 

*****  TIME SPEC  CARD:  SPECIFICATION  OF  COMPUTATIONAL  TINE  UNITS 
VARIABLE  DESCRIPTION  OF  USAGE:  VALUE:  NOTES: 


DELT  COMPUTATIONAL  TIME  STEP  (SECONDS)  S. 
TUNITS  TIME  UNITS  FOR  ALL  INPUT  VARIABLES  SECONDS 
IT1  TIMESTEP  AT  START  OF  MOOEL  SIMULATION  1 
IT2  TOTAL  NUMBER  OF  TINESTEPS  SIMULATED  ABO 
NFREQ  TIMESTEP  INTERVAL  TO  SAVE  GAGE  DATA  6 


ADOTERMS  CARD:  SPECIFICATION  OF  ADDITIONAL  TERMS  IN  GOVERNING  EQUATIONS 


VARIABLE 

DESCRIPTION  OF  USAGE: 

VALUE:  NOTES 

ADVFLAG 

GOVERNS  USE  OF  ADVECTION  TERMS 

NO 

DIFFLAG 

GOVERNS  USE  OF  DIFFUSION  TERMS 

NO 

DIFCOF 

DIFFUSION  COEFFICIENT 

0.0 

CORIOLIS 

CORIOLIS  COEFFICIENT 

0.0 

Figure  14-20.  Sample  output  listing 
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including  one  or  more  PRWINDOW  records.  Variable  VPRVAR  on  record  PRWINDOV  is 
used  Co  specify  which  field  arrays  are  Co  be  princed  (e.g.  VPRVAR  -  E  for 

i 

water  surface  elevaCions).  The  user  can  conCrol  princing  of  field  arrays 
Chrough  Che  course  of  a  simulation  by  specifying  Che  sCarcing  (WPRSTR)  and 
ending  (WPREND)  Cime-sCeps,  CogeCher  wich  Che  Cime  inCerval  (WPRINT) ,  ac  which 
Che  daCa  will  be  prinCed.  FurChemore,  Che  user  can  specify  princing  of 
subgrid  regions  (or  windows)  as  opposed  Co  Che  encire  grid.  This  is  done  by 
specifying  Che  x-  and  y-boundaries  of  a  subgrid  region  wich  variables  VXCEL1, 
WXCEL2,  WYCEL1 ,  and  WYCEL2 . 

138.  Gage  time-histories  of  field  arrays  are  saved  wich  record  RECGAGE. 
One  RECGAGE  record  is  required  for  each  desired  ouCpuC  location  where  Che 
field  array  data  are  Co  be  saved.  A  maximum  of  120  gages,  each  containing  up 
Co  1,000  points,  can  be  processed.  All  gage  time -histories  are  stored  in  one 
file  and  are  processed  by  programs  HYDADD,  HYDLST,  and  HYDPLT  in  package 
CMSPOST . 

139.  Variables  1ST  and  JST  on  record  RECGAGE  define  the  grid  location 
where  field  arrays  are  saved.  The  Cime  inCerval  at  which  these  data  are 
stored  is  controlled  by  variable  MFREQ  on  record  TIMESPEC. 

140.  Record  RECSNAPS  allows  the  user  to  store  wind  and  water  velocity 
data  for  generating  vector  plots.  The  stored  data  can  also  be  used  as  input 
to  other  models.  Data  can  be  saved  at  regular  time  intervals  (SNPTYP  - 
INTERVAL),  at  specific  time-steps  during  a  simulation  (SNPTYP  -  TIMES),  or 
both.  For  data  storage  at  regular  time  intervals,  SNPSTR  defines  the  time- 
step  at  which  data  storage  is  to  begin,  SNPEND  defines  the  time -step  at  which 
data  storage  is  to  end,  and  SNPINT  is  the  regular  time-step  interval  at  which 
data  storage  occurs.  For  data  storage  at  user-specified  times,  one  RECSNAPS 
record  is  required  for  each  time  (SNPTIM)  at  which  a  snapshot  is  to  be 
recorded.  More  than  one  RECSNAPS  record  can  be  used,  and  up  to  100  snapshots 
can  be  saved  in  a  file.  These  data  can  be  processed  with  programs  SNAPCON, 
SNAPLST,  and  SNAPVEC,  which  reside  in  package  CMSPOST. 

141.  Records  XRECRANG  and  YRECRANG  allow  the  user  to  save  data  describ¬ 
ing  discharges  across  arbitrary  transects.  These  data  allow  the  user  to  check 
flow  through  various  tributaries  or  to  determine  the  tidal  prism.  One  record 
is  required  for  each  desired  transect.  A  maximum  of  120  ranges,  each  contain¬ 
ing  up  to  1,000  points,  can  be  processed.  All  range  time-histories  are  stored 
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in  one  file  and  are  processed  by  programs  HYDADD,  HYDLST,  and  HYDPLT  in 
package  CMS POST. 

142.  Similar  to  barriers  and  boundaries,  variable  RP0S1  defines  the  x- 
cell  position  (for  XRECRANG)  or  the  y-cell  position  where  data  are  to  be 
recorded.  Variables  RP0S2  and  RP0S3  define  the  range,  or  extent,  of  the 
recording  range  location.  For  XRECRANG,  records  RP0S2  and  RP0S3  indicate  the 
y-extent,  and  for  YRECRANG,  records  RP0S2  and  RP0S3  indicate  the  x-extent. 

The  time  interval  at  which  these  data  are  stored  is  controlled  by  variable 
NFREQ  on  record  TIMESPEC.  RNAME  is  used  to  name  the  range  data. 
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PART  V:  ILLUSTRATIVE  EXAMPLES 


* 


143.  Three  illustrative  examples  are  included  in  this  section  to 
demonstrate  WICM's  capabilities.  The  first  example  is  a  simple  case  of  waves 
breaking  on  a  plane  beach.  RCPVAVE  is  used  to  supply  the  input  wave  field. 
The  other  two  examples  involve  application  of  UICM  to  Leadbetter  Beach,  CA, 
using  wave  fields  calculated  from  both  RCPVAVE  and  STWAVE.  Examples  demon¬ 
strating  tidal  forcing,  wind- induced  setup,  and  river  discharge  are  given  in 
Chapter  6  (CLHYD  model). 


Plane  Beach  Simulation 

144.  The  plane  beach  simulation  was  executed  on  a  60-  by  8-cell  grid. 
The  grid's  spatial  resolution  was  15  m  in  the  on-offshore  (x-)  direction  and 
30  m  in  the  longshore  (y-)  direction  for  an  overall  grid  size  of  900  by  240  m. 
Water  depth  varied  linearly  from  0.15  m  along  the  shoreline  cells  to  9.00  m  at 
the  offshore  boundary.  The  duration  of  the  simulation  was  40  min. 

145.  The  wave  field  was  calculated  with  RCPVAVE  for  the  plane  beach 
simulation.  RCPVAVE  was  run  with  input  wave  parameters  of  1.4  m  wave  height, 

4  sec  period,  and  30  deg  wave  angle.  Bathymetric  data  used  for  RCPVAVE  are 
the  same  as  required  for  WICM  (see  Table  14-7).  Figure  14-21  shows  the  wave 
height  calculated  from  RCPVAVE  and  the  plane  beach  profile.  Wave  height* is 
approximately  constant  up  to  incipient  breaking,  then  it  decays  rapidly. 
(Instructions  for  RCPVAVE  implementation  are  given  in  Chapter  5.) 

146.  The  WICM  input  data  set  is  given  in  Table  14-7  and  is  discussed 
below.  Wave  input  is  specified  by  the  WAVESPEC  record  (wave  information  from 
RCPVAVE,  built  from  zero  to  the  full  value  over  20  time-steps).  Record 
FRICTABL  specifies  the  constant  friction  coefficient.  The  location  and  type 
of  boundary  forcing  functions  are  specified  with  XB0UNDRY  and  YBOUNDRY 
records.  FUNCTION,  TERECORD,  and  TABELEV  records  give  specifics  about  the 
offshore  boundary  (elevations  in  tabular  form) .  Record  BATHSPEC  is  used  to 
describe  the  bathymetric  data,  and  RECGAGE  records  show  the  locations  of  the 
recording  gages .  The  PRVIND0V  record  was  used  to  print  elevations  and 
velocities  for  the  entire  grid  to  the  output  summary  file  at  10-min  time 
intervals  (120  time-steps). 
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Table  14-7 

Incut  Data  Set  for  the  Plane-Beach  Simulation 


GENSPECS 

WICM  SIMULATION 

NO.  1: 

TIMESPEC 

5. 

SECONDS 

1 

GRIDSPEC 

METRIC 

60 

PRWINDOV 

ADDTERMS 

NO 

NO 

WAVES PECRCPWAVE 

20 

RECGAGE 

1 

4 

RECGAGE 

2 

4 

RECGAGE 

3 

4 

RECGAGE 

4 

4 

PLANE  BEACH,  H-l . 4M , T-4S , A-30DEG  METRIC 


480 

8 


6 

900. 

120 


240. 


1. 


EV 


RECGAGE 

45 

4 

RECGAGE 

50 

4 

RECGAGE 

55 

4 

RECGAGE 

59 

4 

FRICTABL 

0.01 

10.0 

XBOUNDRYCONSTELV 

60 

1 

8 

YBOUNDRYUNI  FLUX 

1 

1 

59 

YBOUNDRYUNI FLUX 

8 

1 

59 

FUNCTION 

1TABELEVS 

TERECORD 

2 

1 

2400. 

TABELEV 

0.0 

0.0 

0.0 

0.0 

TABELEV 

0.0 

40.0 

0.0 

0.0 

BATHS PEC 

1 

0.15 

0.15 

0.15 

0.15 

2 

0.30 

0.30 

0.30 

0.30 

3 

0.45 

0.45 

0.45 

0.45 

4 

0.60 

0.60 

0.60 

0.60 

5 

0.75 

0.75 

0.75 

0.75 

57 

8.55 

8.55 

8.55 

8.55 

58 

8.70 

8.70 

8.70 

8.70 

59 

8.85 

8.85 

8.85 

8.85 

60 

9.00 

9.00 

9.00 

9.00 

1  BDRYX 

BDRY1 
BDRY2 


YX  (8X.9F8.2) 


0.15 

0.15 

0.15 

0.15 

0.15 

0.30 

0.30 

0.30 

0.30 

0.30 

0.45 

0.45 

0.45 

0.45 

0.45 

0.60 

0.60 

0.60 

0.60 

0.60 

0.75 

0.75 

0.75 

0.75 

0.75 

8.55 

8.55 

8.55 

8.55 

8.55 

8.70 

8.70 

8.70 

8.70 

8.70 

8.85 

8.85 

8.85 

8.85 

8.85 

9.00 

9.00 

9.00 

9.00 

9.00 
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Figure  14-21.  RCPWAVE  results  for  plane  beach 

simulation 


147.  The  simulation  was  run  for  480  time-steps  of  5  sec  each  (40  min), 
and  elevation  and  velocity  values  were  saved  every  6  time -steps  at  the  20 
recording  gage  locations.  At  the  conclusion  of  the  simulation,  calculated 
velocities  and  elevations  were  examined  to  check  model  results  (Figures  14-22 
and  14-23).  The  evolution  of  the  water  surface  elevation  and  the  cross-shore 
and  longshore  velocities  with  time  at  grid  cell  (10,4)  are  shown  in  Figure  14- 
22.  Wave  stress  is  built  up  over  20  time-steps  (100  sec),  and  it  takes 
approximately  240  time-steps  for  the  solution  to  reach  steady  state.  The 
small  oscillations  in  5  and  U  can  be  reduced  by  applying  the  wave  forcing 
more  gradually,  but  this  is  unnecessary  if  only  steady- state  results  are 
required.  Figure  14-23  shows  S  ,  U  ,  and  V  for  a  cross-section  of  the 
grid  for  the  final  time-step.  The  cross-shore  velocity  is  zero  because  the 
cross -shore  wave  stress  is  balanced  by  the  water  surface  setup  against  the 
shoreline  and  the  bathymetry  is  homogenous  alongshore.  The  peak  longshore 
velocity  is  2.5  m/sec  at  the  point  of  maximum  wave  height  decay.  The  sharp 
discontinuity  in  the  longshore  current  profile  is  due  to  the  absence  of 
diffusion  (Ah  -  0.0).  The  water  surface  elevation  shows  the  typical  depres¬ 
sion  (setdown)  at  the  break  point  and  superelevation  (setup)  near  the  shore¬ 
line  . 
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148.  In  sumary,  a  longshore  current  was  simulated  for  monochromatic 
wave  breaking  on  a  plane  beach,  neglecting  mixing.  This  simple  example  is  a 
good  test  case  for  first- time  users. 

Leadbetter  Beach.  CA.  Simulation  with  RCPWAVE 

149.  Leadbetter  Beach,  CA,  was  selected  as  an  example  because  field 

data  are  available  at  this  site  from  the  National  Sediment  Transport  Study 
(NSTS)  (Thornton  and  Guza  1986)  for  model  verification.  This  example  uses 
data  collected  on  4  February  1980.  NSTS  data  consist  of  root -mean -square  wave 
height  (  -  0.56  m) ,  peak  period  (14.2  sec),  and  mean  direction  (9.0  deg) 

measured  at  a  water  depth  of  3.8  m  for  model  input  and  wave  height  and 
longshore  current  measurements  across  the  almost -plane  profile  at  Leadbetter 
Beach.  To  apply  the  input  wave  conditions  to  RCPWAVE,  the  wave  height  and 
direction  were  linearly  transformed  to  deep  water  using  Snell's  law,  with  H 
set  equal  to  H ^  .  Also,  the  incipient  breaking  criterion  in  RCPWAVE  was 
changed  to  a  height- to -depth  ratio  of  0.45.  (RCPWAVE  is  formulated  for 
incipient  breaking  of  a  monochromatic  wave.  Since  a  random  wave  field 
consists  of  a  distribution  of  wave  heights,  a  reasonable  incipient  breaking 
criterion  is  an  to -depth  ratio  of  0.4  to  0.5  (Thornton  and  Guza  1986).) 

The  beach  profile  was  assumed  to  be  homogeneous  in  the  alongshore  direction. 
Figure  14-24  shows  the  beach  profile  and  wave  height  across  shore  calculated 
from  RCPWAVE. 

150.  A  WICM  grid  was  constructed  with  a  5-m  grid  spacing  offshore 
(18  cells)  and  10-m  spacing  alongshore  (10  cells).  A  time-step  of  3  sec  was 
used,  and  the  simulation  length  was  400  time-steps.  Wave  stress  was  built  up 
over  20  time-steps.  The  bottom  friction  coefficient  was  set  to  0.007  and  the 
diffusion  coefficient  was  set  to  0.1.  The  WICM  input  data  set  is  given  in 
Table  14-8.  Figure  14-25  shows  the  calculated  velocities  and  surface 
elevations  along  across-shore  row  5  at  the  end  of  the  simulation.  The  lines 
represent  model  results  and  the  symbols  are  longshore  current  measurements. 

The  model  longshore  velocity  results  compare  well  with  the  measurements.  The 
model  was  "calibrated"  to  the  measurements  by  adjusting  the  friction  factor 
and  diffusion  coefficient  (within  reasonable  limits)  to  best  fit  the  data. 

The  current  profile  is  much  smoother  in  this  example  than  in  the  plane  beach 
simulation  because  diffusion  is  included. 
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Table  14-8 

Input  Data  Set  for  the  Leadbetter  Beach  Simulation 


GENSPECS 

WICM  SIMULATION 

NO.  2: 

LEADBETTER,  4 

OCT  1980, 

RCPWAVE 

METRIC 

TIMESPEC 

3. 

SECONDS 

1 

400 

10 

GRIDSPEC 

METRIC 

18 

10 

90. 

100. 

1. 

PRWINDOW 

100 

EV 

ADDTERMS 

NO 

YES 

0.1 

WAVESPECRCPWAVE 

20 

RECGAGE 

1 

5 

RECGAGE 

2 

5 

RECGAGE 

3 

5 

RECGAGE 

4 

5 

RECGAGE 

15 

5 

RECGAGE 

16 

5 

RECGAGE 

17 

5 

RECGAGE 

18 

5 

FRICTABL  0. 

001 

5.0 

XBOUNDRYCONSTELV 

18 

1 

10 

1 

BDRYX 

YBOUNDRYUN I FLUX 

1 

1 

17 

BDRY1 

YBOUNDRYUN I FLUX 

10 

1 

17 

BDRY2 

FUNCTION 

1TABELEVS 

TERECORD 

2 

1 

1200. 

TABELEV 

0.0 

0.0 

0. 

0 

0.0 

TABELEV 

0.0 

20.0 

0. 

0 

0.0 

BATHS PEC 

YX 

(6X.10F6. 

2) 

1 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

3 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

4 

0.72 

0.72 

0.72 

0.72 

0.72 

0.72 

0.72 

0.72 

0.72 

0.72 

5 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

6 

1.14 

1.14 

1.14 

1.14 

1.14 

1.14 

1.14 

1.14 

1.14 

1.14 

7 

1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

8 

1.44 

1.44 

1.44 

1.44 

1.44 

1.44 

1.44 

1.44 

1.44 

1.44 

9 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

10 

1.77 

1.77 

1.77 

1.77 

1.77 

1.77 

1.77 

1.77 

1.77 

1.77 

11 

1.96 

1.96 

1.96 

1.96 

1.96 

1.96 

1.96 

1.96 

1.96 

1.96 

12 

2.15 

2.15 

2.15 

2.15 

2.15 

2.15 

2.15 

2.15 

2.15 

2.15 

13 

2.35 

2.35 

2.35 

2.35 

2.35 

2.35 

2.35 

2.35 

2.35 

2.35 

14 

2.61 

2.61 

2.61 

2.61 

2.61 

2.61 

2.61 

2.61 

2.61 

2.61 

15 

2.89 

2.69 

2.89 

2.89 

2.89 

2.89 

2.89 

2.89 

2.89 

2.89 

16 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

17 

3.53 

3.52 

3.52 

3.52 

3.52 

3.52 

3.52 

3.52 

3.52 

3.52 

18 

3.84 

3.84 

3.84 

3.84 

3.84 

3.84 

3.84 

3.84 

3.84 

3.84 
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Leadbetter  Beach.  CA.  Simulation  with  STWAVE 


151.  To  apply  the  input  wave  conditions  to  STWAVE,  the  input 

measured  at  the  depth  of  3.8  m  was  multiplied  by  a  factor  of  1.4  to  give  the 

input  energy-based  wave  height  H ^  required  by  STWAVE  (for  a  Rayleigh 
distribution  of  wave  heights,  H ^  ~  1.4  Hrms  well  outside  the  surf  zone). 
Input  wave  period,  wave  direction,  and  bathymetric  data  were  the  same  as  used 
in  the  RCPWAVE  simulation.  Instructions  for  STWAVE  implementation  are  given 
in  Chapter  8.  Figure  14-26  shows  the  beach  profile  and  wave  height  across 
shore  calculated  from  STWAVE.  Note  that  the  calculated  wave  height  is  H 
and  the  measured  wave  height  is  H ^  in  Figure  14-26,  so  it  is  expected  that 
calculated  values  will  exceed  measured  values . 

152.  The  WICM  grid  used  for  this  simulation  is  the  same  as  discussed  in 

paragraph  150.  The  WICM  input  data  set  is  similar  to  the  one  given  in 

Table  14-8,  except  "RCPWAVE"  is  replaced  with  "STWAVE"  on  the  WAVESPEC  record 
and  the  bottom  friction  coefficient  is  0.006.  Figure  14-27  shows  calculated 
velocities  and  surface  elevations  along  across-shore  row  5  at  the  end  of  the 
simulation.  The  lines  represent  model  results  and  the  symbols  are  the 
longshore  current  measurements.  Model  results  compare  well  with  measurements. 
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Depth  (m) 


Figure  14-27.  WICM  results  for  Leadbetter  Beach 
simulation  (driven  by  STWAVE) 
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APPENDIX  14 -A:  TRANSFORMATION  OF  ADVECTION  AND  DIFFUSION  TERMS 
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The  advective  (inertia)  terns  in  the  x-  and  y-nonentun  equations  are 
given  as: 


X-Monentun 

d(UU\  +  d  tUV\ 
H  I  *  dy(ff) 


Y-Momentum 

±im  ♦  ±i™\ 

These  expressions  are  represented  in  the  transformed  plane  as: 


X- Momentum 
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Finite  Difference  Representation 

The  upwind  treatment  of  the  advective  terms,  otherwise  known  as  the 
"pigpen  method"  (Roache  1976),  is  implemented  in  the  tfICM  model.  In  this 
formulation,  the  effect  of  a  perturbation  is  carried  only  in  the  direction  of 
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Che  velociCy.  D.  B.  Spalding  called  this  Che  pigpen  method  because  a 
concentrated  quantity  would  be  "smelled”  only  downwind  of  the  source.  This 
translates  into  the  model  as  meaning  that  a  point  in  the  grid  is  affected  only 
by  what  occurs  upwind  of  it.  Advantages  of  this  method  are  its  simplicity  in 
formulation  and  stability  of  the  solution  due  to  its  diffusive  nature.  The 
stability  can  also  act  as  a  disadvantage  if  the  diffusion  becomes  excessive. 

For  flow  in  the  positive  x-direction,  the  finite  differences  are  taken 
between  cells  i-2  and  i  in  accordance  with  the  pigpen  method: 


jk 

SJ* 


°l 

El*  +Z 

d£  B  dr\J 

U  dy  +  V  3y' 

h  31  *  Hd n, 

*1 

l  A* 

E 

H 

1 

+  u( 

__*L _ 1 

gfi 


Udx 
H  di 


For  flow  in  the  negative  x-direction,  the  finite  differences  are  taken  between 
cells  i+1  and  i  in  accordance  with  the  pigpen  method,  and  the  equations 
are  similar  to  those  shown  above. 


Y-Monentua 


For  flow  in  the  positive  y-direction,  finite  differences  are  taken 
between  cells  j-1  and  j  in  accordance  with  the  pigpen  method: 
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Diffusion  terms  in  the  x-  and  y-momentum  equations  are  given  as: 


X- Momentum 
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Y-Mwnfra 


.  ( a*y  &v) 

'la*2  +  ay’J 


These  expressions  are  represented  in  the  transformed  plane  as: 
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Hq  -  inverse  metric  tensor  components 
Of  -  Christoffel  symbols  of  the  second  kind 
(D^LO{  -  derivative  with  respect  to  ( 

(DfU)^  =  derivative  with  respect  to  t) 

(DyV)(  -  derivative  with  respect  to  5 
(Dy  10,  ■  derivative  with  respect  to  i) 

Y- Momentum 

VuHn  ♦  2VlnHn  ♦  V^Hn 

♦  (D^U)tHn  *  (D^U)nHn  ♦  (D^U)^  ♦  (D^U)^ 

♦  (D^V\Hn  +  {D2aV\Hu  ♦  ( n&V)fHtl  ♦  (Dl,V)nHn 

♦  Ui(D1nHn  +  DljHJ  ♦  l/,(DX  ♦  I>X> 

♦  F<0>X  ♦  dX  -  d}1Hu  -  2DljHu  -  I>X) 

-  wX  ♦  0X2) 

+  ia>Xu  ♦  dIp\i  -  *>Xu  -  dWi2>h12  *  (i>X  -  z>XXj2i^ 

♦  ia>Xi  -  *  v>M  *  j>Xu  -  *>X?.  - 


Treatment  of  the  advective  and  diffusive  terms  is  different  near 
boundaries.  The  velocity  cross  derivative  is  assumed  to  be  zero  near  the 
boundary,  and  a  slip  condition  exists.  This  "slip"  boundary  is  assumed  for 
both  the  advection  and  diffusion  terms. 
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APPENDIX  14-B:  WICM  DATA  SPECIFICATION  RECORDS 
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Model  Control  Specifications 


(Req)  GENSPECS 
(Req)  TIMES  PEC 
(Opt)  ADDTERMS 
(Opt)  STARTUP 


(Req)  GRIDSPEC 
(C-Opt)  GRIDCORN 
(C-Opt)  XSTRETCH 
(C-Opt)  YSTRETCH 


(Req)  BATHSPEC 
(Req) 

(Opt)  CHNGBATH 
(Opt)  XBARRIER 
(Opt)  YBARRIER 
(Opt)  FRICTION 
(Opt)  FRICTABL 
(Opt)  CHNGFRIC 


(Opt)  XBOUNDARY 

(Opt)  YBOUNDARY 

(C-Opt)  FUNCTION 
(C-Opt)  CNRECORD 
(C-Opt)  CONSTIT 
(C-Opt)  TERECORD 
(C-Opt)  TFRECORD 
(C-Opt)  TABELEV 


Specify  general  title  and  system  of  units 
Specify  time -related  controlling  variables 
Specify  optional  terms  in  governing  equations 
Specify  initial  conditions 

Grid  Description 

Specify  general  grid  characteristics 
Specify  x-  and  y  comer  points  for  grid 
Specify  x-coordinates  to  create  stretched  grid 
Specify  y-coordinates  to  create  stretched  grid 

Physical  Characteristics 

Specify  characteristics  of  bathymetry/topography 
Two-dimensional  array  of  bathymetric/topographic 
data 

Specify  changes  to  the  bathymetric/topographic  data 

Specify  barrier  perpendicular  to  x-axis 

Specify  barrier  perpendicular  to  y-axis 

Specify  character  of  bottom  friction 

Specify  entry  for  depth-variable  friction  table 

Modify  the  friction  values  at  selected  locations 

Boundary  Conditions 

Specify  driving  boundary  perpendicular  to  x-axis 
Specify  driving  boundary  perpendicular  to  y-axis 
Specify  driving  boundary  forcing  function 
Specify  attributes  of  constituent  forcing 
Specify  harmonic  constituent  forcing  function 
Specify  attributes  of  tabular  elevation  forcing 
Specify  attributes  of  tabular  velocity  forcing 
Specify  irregularly  spaced  tabular  elevations 
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(Req) 

PRWINDOW 

Specify 

(Opt) 

RECGAGE 

Specify 

(Opt) 

RECSNAPS 

Specify 

(Opt) 

XRECRANG 

Specify 

(opt) 

YRECRANG 

Specify 

location  and  timing  of  a  print  window 
location  of  recording  gage  in  grid 
snapshot  time(s)  for  recording 
discharge  range  perpendicular  to  x-axis 
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CMS  Data  Specification:  GENSPECS  Record:  (Req) 

Purpose:  Specify  general  title  and  system  of  units. 
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CMS  Daca  Specification:  TIMESPEC  Record:  (Req) 

Purpose:  Specify  time-related  controlling  variables. 
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CMS  Data  Specif taction:  ADDTERMS  Record:  (Opt) 

Purpose:  Specify  optional  terms  in  governing  equations. 
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CMS  Data  Specification:  STARTUP  Record:  (Opt) 

Purpose:  Specify  inital  conditions. 


a 


o 


as 


XI 

C 

4 

•a 


TJ 

u 

o 

o 

4 

a: 


I 


3 


3 


2  • 
«  o 

H  U 
« 

C  N 
O 

•H  O 
XI  XI 

4 

>  tJ 
V  0) 
r-l  N 

at  —i 


•O 

4 

e 

u 


« 

x> 


o 

4J 


•o 

§ 


4  4 

O  N-t 

4  XI 


U  C 

3  --I 
« 

« 

Xl  Xl 
a>  4 

4J 

4  n 
>  « 


i-4  4J 
4  -4 
—4  O 
4J  O 


C  4 
M  > 


O* 
.  4 
Ml  (4 


4J 

& 


1 7 

cd 

C 


i 


4 

M 

(X 


XI  • 

3  C 
CX  4 

e 

<*4  XI 


•o 

§ 


<XI 

O 


« 

c 

o 


-  c 

M  o 


XI 
U  4 
4  r- 1 
*— 4  2 
jo  S 
Cd 

•H  W 

u 

Cd  AJ 

►  s 

"O  XI 
i-H  4 
4 

•»4  O 
fn  4-1 


XI 

Cl  • 

c  c 

■H  4 
XI 
<4-1  XI 
O  -4 
XI 

tJ  ^ 

O  4 
CX  JO 
4 

X4  O 
4-1 
XI 

Xl  4 
O  4-1 

^■g 


XI  4  XI 

Xl  rCs  -H 

O  T3 

a  o  c 

4  XI  O 

X4  O 

4 

/-v  xi  a 

•o  4  3 

4  *0  XI 

Xl 
4 
XI 

4-14  W 

4  C 

•O  *H  <X4 

0  O 


4 


CO 


fH  0) 
£  & 


•o 

4 

H 

a 

Xl  4 

04 

06 

u 

XI  XI 

<3 

■J 

•H  4 

P 

* 

H 

H 

M 

a  a 

06 

06 

M 

06 

< 

Xi 

< 

H 

O 

H 

4 

H 

O 

X 

W 

PL< 

M 

X 

M 

a 

4c 

< 


I 

X 

M 


XI 

& 


s* 


oo 

00 

\o 

oo  m 

* 

4c 

4c 

4c 

X4  r-4 

Xl 

Xi 

Xi 

4  4 

Xl  4 

4 

4 

jC  4 

O  06 

S€ 

6 

€ 

§ 

£ 

CO 


CN 


CO 


O 

H 

i 

•4* 


14-78 


CMS  Data  Specification:  GRIDSPEC  Record:  (Req) 
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Real  (Req)  R*  Map  sea 


CMS  Data  Specification:  GRIDCORN  Record:  (C-opt) 

Purpose:  Specify  general  characteristics  of  the  curvilinear  grid  coordinate  data. 
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CMS  Data  Specification:  XSTRETCH  and  YSTRETCH  Records:  (C-opt) 

Purpose:  Specify  the  data  to  create  grid  coordinates  in  a  stretched  rectilinear 

cartesian  coordinate  system. 
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9-10  BNAME  Char  *16  (Opt)  A*  Nana  of  bathymetry/ topography  set. 
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CMS  Data  Specification:  XBAKRIER  and  YBARRIER  Records:  (Opt)  < 

Purpose:  Specify  the  location  and  characteristics  of  a  subgrid- scale  barrier. 
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CMS  Data  Specification:  FRICTABL  Record:  (C-Opt) 
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CMS  Data  Specification:  XBOUNDRY  and  YBOUNDRY  Record:  (OPT) 

Purpose:  Specify  location  and  character  of  a  driving  boundary. 
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location  within  the  grid:  Boundary 

BNP0S2  Integer  (Req)  +1*  extends  from  (and  includes)  cells 

BNP0S2  to  BNP0S3  along  cell  BNP0S1. 

BNP0S3  Integer  (Req)  +1* 


BNDFN1  Integer  (C-opt)  +1*  Integer  index  of  forcing  function 

(tabular  or  harmonic  constituent)  for 
CONSTELV,  INTRPELV  or  CONSTDIS  type 
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CMS  Data  Specification:  FUNCTION  Record:  (C-opt) 

Purpose:  Specify  index  number  and  character  of  driving  boundary  forcing  function. 
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FDATUM  Real  (Opt)  0.0  R*  The  function  values  are  added  to  this 

"datum"  quantity. 


CMS  Data  Specification:  CNRECORD  Record:  (C-opt) 

Purpose:  Specify  physical  coordinates  and  timing  of  a  constituent  forcing  function. 
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CMS  Data  Specification:  CONSTIT  Record:  (C-opt) 

Purpose:  Specify  and  quantify  a  harmonic  constituent  for  a  boundary  forcing  function. 
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This  record  must  follow  a  FUNCTION  record  if  TABELEVS  or  TAB FLOWS  specified  for  FUNTYP. 

If  the  tabular  data  are  regularly- spaced  (in  time),  a  1-d  array  must  follow  this  record. 

If  the  tabular  data  are  provided  at  irregular  time  intervals,  one  or  more  TABELEV  or  TABFLOW 
records  must  follow  this  record. 
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CMS  Data  Specification:  WINDSPEC  Record:  (Opt) 

Purpose:  Specify  the  character  of  wind  field  data. 
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TABWINDS  record  is  followed  by  wind 
speed  arrays . 


WUNITS  Char  *8  (Opt)  FPS  MPH  Units  for  wind  values. 

FPS 
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CMS  Data  Specification:  VAVESPEC  Record:  (Opt) 

Purpose:  Specify  the  character  of  wave -field  data. 
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CMS  Data  Specification:  PRWINDOW  Record:  (Opt) 

Purpose:  Specify  location  and  timing  of  a  print  window. 
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CHS  Data  Specification:  RECGAGE  Record:  (Opt) 

Purpose:  Specify  location  and  character  of  a  recording  gage  in  the  grid. 


CMS  Data  Specification:  RECSNAPS  Reocrd:  (Opt) 

Purpose:  Specify  snapshot  tiae(s)  for  recording. 
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Alternate  form  for  specific  times  - 

Permitted 
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APPENDIX  14 -C:  INPUT  DATA  SET  FOR  THE  INDIAN  RIVER  INLET  EXAMPLE 
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GFWSPECS  INDIAN  RIVER  INLET  ENGLISH 


GRIDSPECCURVILIN 

ENGLISH 

84 

84  66666.6 

66666.6  1.0 

GRIDCORN  XY 

3333.33301 

(F12.5.4F16 

5000.00000 

■  5> 

6666.66699 

8333.33301 

10000.0000 

11666.6699 

13333.3301 

15000.0000 

16666.6699 

18333.3301 

20000.0000 

21666.6699 

23333.3301 

25000.0000 

26666.6699 

28333.3301 

30000.0000 

31666.6699 

33333.3281 

35000.0000 

36666.6602 

38106.6719 

39183.3281 

40000.0000 

40666.6602 

41250.0000 

41770.0000 

42226.6719 

42636.6719 

43003.3281 

43333.3281 

43633.3281 

43906.6602 

44160.0000 

44396.6719 

44613.3281 

44813.3281 

45000.0000 

45186.6602 

45390.0000 

45610.0000 

45846.6602 

46100.0000 

46373.3281 

46666.6602 

TIMES PEC 

20.  SECONDS 

1  11340 

90 

ADDTERMS 

YES 

NO 

0.0 

0 . 00009 

YBOUNDRYCONSTELV 

1 

70  83 

1 

XBOUNDRYCONSTELV 

84 

1  84 

1 

YBOUNDRYCONSTELV 

84 

68  83 

1 

FUNCTION 

1TABELEVS 

FEET 

0.5 

TERECORD 

127 

1 

1800.  (15F5.2) 

1.80  1.90 

2.00  2.11 

2.42 

2.42  2.39  2.23 

1.90  1.54  1.12  0.72  0.13  -.56-1.15 

-1.61-2.07- 

-2.30-2.42- 

-2.39- 

2.29-1.97-1.74- 

-1.41-0.98-0.49-0.10  0.30  0.59  0.69 

0.76  0.66 

0.06  0.36 

0.13 

-.26  -.62-1.18- 

-1.61-2.07-2.36-2.49-2.52-2.26-2.00 

-1.61-1.18 

-.69  -.10 

0.56 

1.12  1.77  2.20 

2.62  2.78  2.88  2.79  2.49  2.20  1.77 

1.28  0.66 

0.16  -.46 

-1.12- 

1.67-2.07-2.26- 

-2.29-2.19-1.90-1.64-1.32  -.89  -.26 

0.16  0.59 

0.98  1.12 

1.25 

1.25  1.18  0.95 

0.66  0.46  0.00  -.43-1.05-1.54-1.94 

-2.19-2.26- 

-2.32-2.10 

-1.74- 

■1.35  -.85  -.39 

0.23  0.85  1.51  2.03  2.55  2.88  3.05 

2.92  2.82 

2.62  2.30 

1.84 

1.35  0.62  -.16 

-.89-1.44-1.94-2.16-2.19-2.23-2.00 

-1.77-1.38 

-.98  -.69 

-.10 

0.26  0.76  0.00 

0.00 

WINDSPECTVRUNISP 

FPS 

TABWINDS 

0 

0 

-2.201  11.077 

TABWINDS 

0 

1 

-2.096  12.884 

TABWINDS 

0 

2 

-0.763  13.472 

TABWINDS 

0 

3 

2.750  13.210 

TABWINDS 

0 

4 

5.101  8.570 

TABWINDS 

0 

5 

3.952  7.200 

TABWINDS 

0 

6 

5.934  5.022 

TABWINDS 

0 

7 

5.988  3.088 

TABWINDS 

0 

8 

7.783  2.120 

TABWINDS 

0 

9 

6.150  0.348 

TABWINDS 

0 

10 

4.494  -0.692 

TABWINDS 

0 

11 

6.008  -1.361 

TABWINDS 

0 

12 

4.059  -0.625 

TABWINDS 

0 

13 

3.290  -0.745 

TABWINDS 

0 

14 

3.857  -2.117 

14-110 


TABWINDS 

0 

15 

5.194 

-2.396 

TABUINDS 

0 

16 

2.875 

2.031 

TABWINDS 

0 

17 

-3.090 

0.930 

TABUINDS 

0 

18 

-6.266 

2.501 

TABWINDS 

0 

19 

-5.526 

5.876 

TABWINDS 

0 

20 

-6.123 

7.497 

TABUINDS 

0 

21  - 

10.286 

3.899 

TABWINDS 

0 

22 

-8.928 

6.916 

TABWINDS 

0 

23 

-5.832 

4.200 

TABWINDS 

1 

0 

-3.526 

4.317 

TABWINDS 

1 

1 

10.423 

1.696 

TABWINDS 

1 

2 

11.877 

0.257 

TABWINDS 

1 

3 

10.183 

3.755 

TABWINDS 

1 

4 

6.334 

2.720 

TABWINDS 

1 

5 

4.485 

1.860 

TABWINDS 

1 

6 

4.708 

1.643 

TABUINDS 

1 

7 

7.489 

3.372 

TABWINDS 

1 

8 

6.713 

-0.677 

TABWINDS 

1 

9 

4.413 

1.987 

TABWINDS 

1 

10 

0.535 

3.924 

TABWINDS 

1 

11 

7.816 

3.357 

TABWINDS 

1 

12 

5.941 

4.044 

TABWINDS 

1 

13 

6.841 

3.371 

TABWINDS 

1 

14 

7.104 

2.341 

TABWINDS 

1 

15 

8.083 

1.460 

TABWINDS 

1 

16 

7.517 

-3.309 

TABWINDS 

1 

17 

3.155 

-6.129 

TABWINDS 

1 

18 

-6.163 

-1.913 

TABWINDS 

1 

19 

-7.929 

-3.082 

TABWINDS 

1 

20 

-9.665 

-2.999 

TABWINDS 

1 

21 

-7.295 

1.652 

TABWINDS 

1 

22 

-7.072 

3.881 

TABWINDS 

1 

23 

-6.039 

1.819 

TABWINDS 

2 

0 

-6.660 

4.286 

TABWINDS 

2 

1 

-3.646 

8.490 

TABWINDS 

2 

2 

-2.020 

7.961 

TABWINDS 

2 

3 

3.413 

4.406 

TABWINDS 

2 

4 

8.600 

1.866 

TABWINDS 

2 

5 

10.069 

-1.016 

TABWINDS 

2 

6 

4.721 

-5.990 

TABWINDS 

2 

7 

5.169 

-4.780 

TABWINDS 

2 

8 

-0.551 

-11.427 

TABWINDS 

2 

9 

-3.668 

-11.454 

TABWINDS 

2 

10 

-8.492 

-13.718 

TABWINDS 

2 

11 

-5.476 

-9.200 

TABWINDS 

2 

12 

-2.071 

-7.340 

TABWINDS 

2 

13 

-3.808 

-8.250 

TABWINDS 

2 

14 

-8.427 

-9.776 

TABWINDS 

2 

15 

-7.454 

-9.624 

BATHSPEC 

FEET 

0.0 

0.0 

YX(10F8 . 1) 

8.5 

0.0 

0.0 

0.0 

9.5 

9.4 

9.6  10.1 

7.4 

10.0 

8.5 

8.5 

8.5 

8.5 

8.5 

8.4  8.4 

7.9 

7.8 

7.5 

7.5 

7.4 

.4 

7.4 

7.4 

7.9  7.8 

7.8 

7.8 

7.8 

7.8 

7.7 

7.7 

7.7 

7.7 

7.7  7.6 

7.6 

7.6 

7.6 

7.2 

7.2 

7.2 

6.7 

6.7 

7.2  7.2 

7.1 

7.1 

7.1 
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7.1 

7.1  7.0 

7.0 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0.0  0.0 

0.0 

39.0 

38.0  37.0 

38.0 

41.5 

40.3 

42.0  40.0 

39.9 

35.0 

38.1 

38.0  38.0 

37.9 

37.0 

37.1 

37.2  37.2 

37.2 

39.2 

38.6 

38.6  38.7 

32.5 

31.9 

30.6 

38.0  31.9 

32.5 

32.6 

38.3 

38.4  38.8 

38.0 

37.8 

43.9 

45.0  44.1 

45.0 

49.0 

50.0 

47.4  48.0 

46.0 

38.0 

41.0  40.0 

41.0 

41.2 

40.4 

38.0  41.1 

42.4 

45.0 

42.1 

42.0  42.3 

38.0 

41.1 

41.1 

41.1  40.0 

40.0 

39.9 

40.5 

35.5  35.1 

34.6 

39.0 

33.4 

33.7  38.5 

38.5 

38.3 

38.5 

38.6  38.9 

40.8 

40.9 

50.0 

48.7  48.7 

49.4 

49.7 

51.0 

51.0  50.0 

44.0 

FRICTION 

TABULAR 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.033 

0.033 

0.033 

0.033 

0.033  0.033  0.033 

0.033 

0.033 

0.033 

0.033  0.033  0.033 

0.033 

0.033 

0  033 

0.033  0.033  0.033 

0.033 

0.033 

0.033 

0.033  0.033  0.033 

0.033 

0.033 

0.033 

0.033  0.033  0.033 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045  0.045  0.045 

0.033 

0.045 

0.045 

0.045  0.045  0.033 

0.033 

0.033 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045 

0.045  0.045  0.045 

0.045 

0.045 

0.045r 

0.045  0.045  0.045 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

41.3 

40.0 

39.0 

39.4 

39.0 

39.0 

39.0 

39.0 

38.9 

38.2 

37.0 

37.0 

37.0 

37.1 

37.1 

38.2 

38.4 

38.4 

38.5 

38.5 

31.6 

31.5 

34.0 

31.5 

30.2 

37.2 

37.0 

37.6 

37.6 

37.9 

40.0 

39.9 

40.3 

43.7 

47.0 

48.0 

48.5 

47.7 

47.9 

50.0 

40.0 

38.9 

38.0 

38.0 

39.0 

41.8 

41.8 

41.7 

42.2 

42.1 

41.1 

41.1 

41.1 

41.1 

41.1 

39.9 

40.0 

40.6 

40.6 

40.6 

33.1 

32.9 

39.1 

33.0 

33.2 

38.3 

38.2 

38.3 

38.3 

38.0 

41.7 

43.8 

48.0 

46.6 

47.0 

51.5 

51.0 

49.0 

50.7 

51.0 

YX(10F7 . 3) 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.033 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.033 

0.033 

0.033 

0.033 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 
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0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

0.020 

0.020 

0.020  0.020 

0.020 

RECGAGE 

70 

27 

RECGAGE 

70 

28 

RECGAGE 

70 

29 

RECGAGE 

36 

49 

RECGAGE 

44 

48 

RECGAGE 

42 

48 

RECGAGE 

38 

49 

RECGAGE 

6 

5 

RECGAGE 

16 

45 

RECGAGE 

36 

57 

RECGAGE 

58 

81 

RECGAGE 

63 

32 

RECGAGE 

43 

48 

RECGAGE 

37 

50 

RECGAGE 

35 

50 

RECGAGE 

69 

27 

RECGAGE 

69 

28 

RECGAGE 

69 

29 

RECSNAPS INTERVAL 

RECRANG 

XRECRANG 

69 

26 

RECRANG 

XRECRANG 

45 

48 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

INLET  GAGE  3 
INLET  GAGE  2 
INLET  GAGE  1 
MSTAT4 
MSTAT1 
MSTAT2 
MSTAT3 
VINES  GAGE 
POTNETS  GAGE 
MASSEY  GAGE 
DEWEY  GAGE 
INLET  GAGE 
MSTAT1 
NSTAT3 
MSTAT4 
INLET  GAGE3 
INLET  GAGE2 
INLET  GAGE1 

720  7920  9360 

29 
53 
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